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1.  SUMMARY 


The  Advanced  NDE  Sensor  Modeling  program  has  developed  a  finite  element  method 
(FEM)  of  electromagnetic  NDE  for  optimizing  the  inspection  of  aging  aircraft  structure.  The 
significant  results  of  the  program  include  1)  excellent  correlation  of  electromagnetic  FEM 
predictions  with  MR  sensor  experimental  data  for  aircraft  fastened  structure  configurations,  2) 
evaluation  of  parameter  variation  that  affect  the  signal  to  noise  in  the  resulting  data,  and  3) 
optimized  MR  scanning  of  test  samples  showed  an  improvement  in  sensitivity  to  second  layer 
flaws  by  over  an  order  of  magnitude  compared  to  standard  eddy  current  inspection.  MR 
sensors  are  found  to  provide  robust  signals  of  the  magnetic  field  strengths  in  electromagnetic 
testing.  The  data,  however,  can  be  subject  to  interpretation  difficulties  due  to  the  field 
patterns  from  geometric  features. 

The  FEM  approach  accurately  represented  the  fields  and  response  of  MR  sensors  from 
AC  currents  induced  in  metal  geometries  for  cracks  near  fasteners  in  the  second  layer  under 
thick  (>0.16  inch)  skins.  Once  satisfactorily  calibrated,  the  model  was  used  to  conduct  a 
variety  of  parametric  studies  including  frequency,  sensor  liftoff,  conductivity  of  layers, 
fastener  to  edge  distance  and  crack  dimension.  Additionally,  meticulous  analysis  of  MR 
sensor  measurements  was  performed  to  development  automated  feature  recognition.  The  data 
was  further  processed  to  place  an  objective  value  to  cracks  as  they  appear  in  MR  sensor  data. 
As  a  result  of  quantitative  data  reduction,  it  has  been  possible  to  establish  a  signal-to-noise 
ratio  (SNR)  value  for  the  MR  sensing.  With  this  evaluation  it  was  possible  to  show  that  MR 
sensors  offer  considerable  value  over  conventional  eddy  current  inspections  with  the 
possibility  of  improvements  in  sensitivity  by  over  an  order  of  magnitude. 

The  results  of  this  program  validate  the  importance  of  MR  sensors  for  deep,  second  layer 
inspections  in  aging  aircraft  and  should  lead  to  cost  savings  through  reduced  teardown  and 
faster  inspection  at  improved  sensitivity  over  presently  applied  methods.  Further  efforts  are 
warranted  to  enhance  the  usefulness  of  MR  sensors  through  reverse  modeling  for  sizing, 
evaluation  of  a  greater  number  of  parametric  effects  that  could  not  be  covered  in  this  program, 
and  investigating  MR  sensor  orientation  and  geometry  variations.  In  addition  to  this  final 
report,  the  program  has  generated  several  publicly  released  reports  as  listed  in  Table  1. 
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Table  1  Publicly  Released  Reports  from  the  Advanced  NDE  Sensor  Modeling  for  Multi- 

Site  Inspection  Program 


Title/Authors 

Publication 

“Automated  Analysis  of  Eddy  Current  Giant 
Magnetoresistive  Data,  “  Y.  Deng,  X.  Liu,  Z.  Zeng,  B. 
Koltenbah,  R.  Bossi,  G.  Steffes,  and  L. 

Udpa, 

Review  of  Progress  in 
Quantitative  Nondestructive 
Evaluation,  (28)2008  (Inpress) 

“GMR  signal  analysis  for  Detecting  Cracks  Under 
Fastener  Heads  in  Multilayer  Aircraft  Structures,”  Z. 
Zeng,  Y.  Deng,  X.  Liu,  L.  Udpa,  S.  S.  Udpa,  B. 
Koltenbah,  R.  Bossi,  and  G.  Steffes 

IEEE  Transactions  on 

Magnetics  (under  review) 

“Feature  Based  Analysis  of  Magnetoresistance  Sensor 
Data  for  Nondestructive  Testing  Applications,”  Y. 

Deng,  X.  Liu,  Z.  Zeng,  B.  Koltenbah,  R.  Bossi,  G. 
Steffes  and  L.  Udpa 

13  th  International  Workshop  on 
Electromagnetic  Nondestructive 
Evaluation,  Korea,  2008 

“Electromagnetic  Modeling  of  Magnetoresistive 
Sensors,”  B.  Koltenbah,  R.  Bossi,  G.  Steffes,  L.  Udpa, 
Y.  Deng,  X.  Liu  and  Z.  Zeng 

ASNT  Fall  Conference  and 
Quality  Testing  Show,  Las 

Vegas,  Nov.  2007 

2.  INTRODUCTION 

2.1  AGING  AIRCRAFT  NDE  NEEDS 

The  development  of  NDE  methods  is  vital  in  determining  the  state  of  health  of  aging 
aircraft  materials.  With  continuing  reductions  in  new  weapons  systems  procurements,  there 
has  been  and  continues  to  be  a  need  to  increase  the  service  life  of  aircraft  currently  in  the  Air 
Force  inventory.  This  requirement  has  made  the  detection  and  characterization  of  corrosion 
and  cracking  a  major  Air  Force  logistics  needs.  Table  2  shows  a  listing  of  several  types  of  Air 
Force  aircraft  and  their  average  age[1,2].  Design  lifetimes  for  these  aircraft  are  typically  20  to 
30  years.  Some  of  the  aircraft  referenced  in  Table  2  have  more  than  exceeded  their  design 
lifetimes.  In  fact,  the  current  strategies  include  extending  the  service  lives  of  aircraft  such  as 
the  B-52  and  KC-135  up  to  80  years. 
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Table  2  Air  Force  Supported  Aircraft [2] 


Aircraft  Type 

Number  of 
Aircraft 

Average  Age 
(Years) 

Comments 

A-10 

196 

26 

Retain  25+  years  -  no  replacement 
identified[1] 

B-1B 

67 

20 

B-52 

85 

46 

No  replacement  identified 

C-5 

33 

21 

C17 

153 

7 

C130 

277 

33 

E-3 

32 

28 

E-4 

4 

33 

F-15  C-D 

292 

24 

Retire  -  replace  with  F-22 

F-15E 

223 

16 

F-16 

700 

16 

Retire  -  replace  with  Joint  Strike 
Fighter  (JSF) 

KC-10 

59 

23 

KC-135 

425 

46 

Replacement  over  time  with  some 
aircraft  reaching  80  years 

T-37 

87 

43 

T-38 

462 

40 

For  aging  aircraft  inspection  the  inspection  technique  needs  to  be  readily  fieldable,  one¬ 
sided,  applicable  to  variable  materials  and  geometries,  and  sensitive  to  defects  in  sublayers. 
Table  3  compares  the  advantages  and  limitation  of  the  general  categories  of  NDE 
methodologies.  Of  these  methods,  it  is  evident  that  the  emphasis  of  study  for  aging  aircraft 
NDE  should  be  the  eddy  current  method  because  it  is  the  only  approach  that  is  relatively  low 
cost,  one-sided,  and  has  the  ability  to  penetrate  multiple  layers.  Eddy  current  NDE  methods 
create  electromagnetic  fields  through  some  form  of  alternating  current  flow.  These  fields 
penetrate  into  the  aircraft  material  and  induce  eddy  currents,  which  themselves  generate  a 
responding  electromagnetic  field.  The  superposition  of  these  electromagnetic  fields  is  then 
measured.  Changes  in  the  electromagnetic  field  can  be  readily  detected  by  the  eddy  current 
device  provided  the  detector  has  sufficient  sensitivity.  Standard  eddy  current  testing  (ECT) 
methods  using  wound  coils,  however,  suffer  from  limitations  due  to  the  principles  of  coil 
technology  for  frequency  of  operation  and  physical  size.  This  results  in  significant  limitations 
in  sensitivity  as  a  function  of  depth  of  penetration  for  sublayer  defect  detection.  Lower 
frequencies  required  for  penetration  use  large  coils  which  decrease  resolution.  These 
limitations  can  be  overcome  by  using  advanced  electromagnetic  sensors  and  modeling  studies 
to  improve  the  detection  sensitivity  to  meet  aging  aircraft  inspection  needs. 
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Table  3  Advantages  and  Limitations  of  NDE  Methods  for  Aging  Aircraft 


NDE  Methodology 

Advantages 

Limitations 

VT  -  visual  and/or 
Enhanced  Visual 
Testing 

Low  cost, 

Easy  to  implement 

Surface  only 

For  subsurface  requires  invasive 
measures 

ET  -  eddy  current 
testing 

One  sided 

Low  cost 

Multi-layer  sensitivity 

Sensitivity  at  depth 

Interpretation  of  signals 

UT  -  ultrasonic 
testing 

One  sided 

Low  cost 

First  layer  (marginal  or  no 
sensitivity  to  second  layer) 

Couplant  required 

RT  -  radiographic 
testing 

Multi  layer  sensitivity 

High  resolution 

Two  sided 

Backscatter  technique  for  one¬ 
sided  is  too  slow  and  lacks 
sensitivity  to  small  defects 

IRT  -  thermographic 
testing 

Image  display 

Non-contact 

First  layer,  no  known  sensitivity 
to  second  layer  features 

PT  and  MT  - 
penetrant  and 
magnetic  particle 

One  sided 

Low  cost 

First  layer  only  for  surface  or 
near  surface  defects 

2.2  ELECTROMAGNETIC  NDE 

Electromagnetic  NDE,  especially  in  the  form  of  ECT,  is  the  most  widely  used  inspection 
method  for  metallic  aircraft  inspection.  Of  paramount  importance  is  the  sensitivity  of  the 
electromagnetic  methods  to  aging  aircraft  defects  in  multilayered  geometries.  As  noted 
above,  ECT  suffers  from  limitations  for  detecting  defects  in  multilayer  complex  structure 
configurations  due  to  the  coil  size  required  for  the  penetrating  frequencies  of  operation  in 
thick  structure.  The  standard  eddy  current  method  easily  becomes  limited  in  sensitivity  as  a 
function  of  the  depth  of  penetration.  This  is  fundamental  to  the  physics  of  the 
electromagnetic  field  problem  and  the  interactions  of  the  induced  currents  with  the  materials 
and  defects,  and  the  interaction  of  the  altered  fields  with  the  sensors. 

Advanced  electromagnetic  sensing  methods,  such  as  MR  sensors,  have  been  developed 
that  offer  significant  advantages  over  traditional  ET  sensors  for  sensitivity  at  depth.  The  Air 
Force  Research  Laboratory  NDE  Branch  funded  a  series  of  programs  to  develop  MR  sensors 
into  useful  tools  for  the  in-service  inspection  of  at  aircraft  that  would  provide  improved 
capability  over  currently  applied  methods [3,4,5].  However,  the  implementation  of  MR  sensors 
is  subject  to  a  great  number  of  variables  when  applied  to  multi-site  geometric  configurations 
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and  defect  types.  These  include  variables  due  to  the  method,  such  as  frequency,  beam  size, 
excitation,  sample  spacing,  etc.,  variables  due  to  the  object  such  as  materials,  sizes,  shapes, 
configuration,  edges,  gaps,  etc.  and  variables  due  to  the  defects  such  as  length,  width,  volume, 
shape,  distribution,  orientation,  etc.  All  three  sets  of  variables  will  be  combined  in  any 
inspection  causing  uncertainty  in  the  interpretation  of  the  signals.  The  result  is  that  the 
sensitivity  to  defect  types,  sizes  and  orientations  is  rarely  optimized.  Intuitive  or  empirical 
approaches  to  optimize  the  inspections  are  inefficient  and  expensive  because  of  the  large 
number  of  variables  that  affect  the  sensitivity.  By  employing  modeling  of  the  electromagnetic 
field  generation  and  the  interaction  of  the  fields  with  materials,  geometries,  features  and 
sensors  the  electromagnetic  NDE  sensitivity  can  be  assessed  and  inspection  parameters 
optimized  efficiently.  Modeling  of  the  physics  problem  provides  the  best  approach  for 
optimization  and  will  result  in  improved  probability  of  detection  (POD)  based  on  how  each 
kind  of  sensor  system  will  behave  for  a  given  geometry,  excitation  method,  crack  or  corrosion 
size,  and  orientation. 

2.3  MODELING 


Electromagnetic  modeling  has  also  been  under  development  and  studied  for  a  number  of 
years,  particularly  at  Michigan  State  University  (MSU)  for  NDE  applications  and  at  Boeing 
for  general  physics  problems.  Models  can  be  broadly  classified  as  either  analytical  or 
numerical.  The  early  work  of  Dodd  and  Deeds,  who  derived  closed  form  integral  expressions 
for  the  impedance  of  a  coil  placed  above  a  layered  half-space  conducting  medium,  set  the 
stage  for  the  development  of  a  series  of  models  characterizing  eddy  current  NDE 
phenomena^.  Although  these  models  continue  to  be  used  extensively  to  this  day,  a  number 
of  extensions  have  been  proposed  in  recent  years  using  volume  integral  and  boundary  integral 
methods.  However,  such  analytical  models  are,  in  general,  limited  in  their  application  to 
simple  specimen  and  defect  geometries.  Numerical  models  have  gained  popularity  in  recent 
years  due  to  their  ability  to  simulate  arbitrary-shaped  defects  in  complex  sample  geometries. 

Numerical  models  can  help  in  visualizing  the  electromagnetic  field  distribution  and  the 
manner  in  which  the  field  interacts  with  the  flaw  in  the  test  specimen.  Such  models  can  be 
used  to  conduct  systematic  parametric  studies  as  the  material,  defect  and  experimental 
parameters  are  varied.  The  knowledge  of  the  qualitative  and  quantitative  influence  of  these 
parameters  on  the  sensor  measurement  can  be  useful  in  determining  improvements  needed  in 
sensor  and  system  design  and  for  maximizing  the  probability  of  detection  (POD)  of  critical 
flaws  in  structures.  Numerical  models  can  also  serve  as  a  test  bed  for  generating  defect 
signatures  for  a  variety  of  defects.  These  defect  signatures  can  be  employed  for  developing 
inversion  algorithms  as  well  as  for  training  inspection  personnel.  Finally,  computational 
models  can  be  used  to  explore  the  feasibility  of  novel  inspection  methods  and  evaluate  new 
types  of  sensors. 
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2.4  ELECTROMAGNETIC  SENSORS 


2.4.1  Eddy  Current  Sensors 

Eddy  current  NDE  techniques  are  based  on  the  principles  of  electromagnetic  induction. 
These  techniques  find  widespread  use  for  inspecting  aircraft  structures  largely  due  to  their 
noncontact  nature,  simplicity  of  operation  and  sensitivity  to  anomalies  in  thin  multilayer 
conducting  structures.  The  basic  eddy  current  probe  consists  of  a  coil  excited  by  an 
alternating  current  source.  When  the  probe  is  brought  in  close  proximity  to  a  conducting 
sample,  eddy  currents  are  induced  in  the  test  specimen.  The  secondary  field  generated  by  the 
induced  currents  alters  the  net  flux  linkages  associated  with  the  probe  coil.  This  change  in 
flux  linkages  is  reflected  as  a  change  in  the  terminal  impedance  of  probe.  The  impedance 
change  of  the  probe  coil,  as  the  sample  is  scanned,  constitutes  the  eddy  current  signal.  Coil 
type  transducers  offer  the  advantages  of  being  very  accurate,  robust  and  simple  to  construct. 
Unfortunately,  the  measurements  tend  to  be  somewhat  noisy  since  the  output  voltage  is 
related  to  the  time-derivative  of  the  field.  However,  this  is  seldom  a  concern  in  the  case  of 
conventional  single  frequency  eddy  current  methods,  since  the  coherent  detection  scheme 
employed  in  most  eddy  current  instruments  offers  excellent  noise  rejection  characteristics. 
Major  issues  in  eddy  current  inspection  include  the  limited  depth  of  penetration  and  the 
challenges  involved  in  interpreting  the  complex  probe  impedance  signal. 

Eddy  current  inspection  systems  have  grown  in  sophistication  with  advancements  in  the 
fields  of  sensors,  microelectronics  and  computers.  Concurrently,  advancements  in  our  ability 
to  model  the  underlying  physical  processes  give  us  the  ability  to  develop  novel  probe 
designs1-71.  Sensors  that  measure  the  field  due  to  induced  currents  directly  as  opposed  to  coil 
sensors  that  measure  the  derivative  of  the  field  are  attracting  attention.  Giant 
magnetoresistive  (GMR)  probes18’91,  Hall  sensors,  fluxset  probe[10],  magneto-optic 
sensors11 1,121  and  meandering  wire  magnetometers1131  are  all  examples  of  sensors  being 
considered  seriously  as  alternatives  to  the  classical  eddy  current  coil.  Advantages  offered  by 
these  sensors  range  from  high  levels  of  sensitivity  to  small  size. 

2.4.2  GMR  Sensor 

Among  sensors  that  are  capable  of  measuring  the  field  due  to  induced  currents  directly, 
GMR  probes  appear  to  be  the  most  promising.  These  probes  offer  exceptional  levels  of 
sensitivity,  small  size  and  low  cost.  A  simple  GMR  probe  configuration  is  shown  in  Figure  1, 
where  the  excitation  coil  encircles  the  GMR  sensor^141.  The  magnetic  field  generated  by  the 
excitation  coil  induces  currents  in  the  specimen.  The  induced  currents  are  distorted  if  a  crack 
is  encountered.  The  GMR  sensor  picks  up  perturbations  in  the  fields  associated  with  the 
induced  eddy  currents.  Signal  processing  algorithms  may  be  required  to  demodulate  the  signal 
detected  by  the  GMR  sensor. 
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Figure  1  Experiment  Setup  with  Conventional  Eddy  Current  Testing 

Figure  2  shows  an  alternative  approach  to  using  the  GMR  sensor  where  the  sinusoidal 
excitation  current  is  applied  to  an  induction  foil.  An  array  of  GMR  sensors  can  be  used  to 
pick  up  the  normal  component  of  magnetic  flux  density  generated  by  the  induced  currents.  In 
the  absence  of  any  discontinuity,  the  magnetic  flux  is  tangential  to  the  specimen  surface 
resulting  in  a  null  signal.  Anomalies  in  the  specimen  result  in  a  normal  component  of  the 
magnetic  flux  density.  The  normal  component  of  the  magnetic  field  is  sensed  by  the  GMR 
sensor.  Alternately,  a  pulse  excitation  current  can  be  applied  to  the  induction  foil  and  an  array 
of  GMR  sensors  can  be  used  to  measure  the  transient  magnetic  field  at  different  locations 
simultaneously.  The  output  signals  from  the  GMR  sensor  array  contain  time  and  space 
information.  The  information  can  be  used  to  estimate  the  defect  depth  and  location. 


1  D  array  of  MR  sensors 
it  current 
sample 


Figure  2  MR  Sensor  Array  with  Sheet  Current  Excitation  Foil 

2.5  ADVANCED  NDE  MODELING  PROGRAM  BENEFIT 

The  Advanced  NDE  Sensor  Modeling  for  Multisite  Inspection  program  was  a  three  year 
effort  concentrated  on  the  development  of  appropriate  electromagnetic  models  suited  to  MR 
sensors  to  determine  the  optimum  sensitivity  capability.  The  investigation  modeled  standards 
that  represent  damage  situations  identified  in  the  Air  Force  and  at  Boeing  as  serious  concerns 
involving  cracking  in  sublayers  of  thick  (>0.25  inch)  metallic  structure.  Testing  and 
correlation  of  real  data  with  modeled  data  was  performed  to  validate  and  improve  the 
modeling.  The  verified  models  were  used  to  predict  the  optimum  parameters  for  MR 
inspection  of  simulated  aging  aircraft  structure  configurations.  The  influence  of  the 
complexities  of  aging  aircraft  inspection,  such  as  multiple  layers,  thickness  variations,  gaps, 
material  changes,  defect  size  and  defect  orientation  were  examined  and  analyzed  by  the 
modeling  studies.  The  improvement  in  defect  detection  and  sizing,  through  the  use  of 
advanced  sensors  with  optimized  model  predicted  inspection  parameters,  was  assessed 
compared  to  conventional  eddy  current  and  found  to  be  about  an  order  of  magnitude  better. 
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3.  METHODS  ASSUMPTIONS  PROCEDURES 


3.1  PROGRAM  PLAN 

The  program  plan  for  “Advanced  NDE  Sensor  Modeling  for  Multi-Site  Inspection,”  is 
shown  in  Figure  3.  The  program  schedule  is  shown  in  Figure  4.  The  program  involved  the 
cooperative  efforts  of  Boeing  Phantom  Works  and  MSU  Electrical  and  Computer  Engineering 
Department.  Boeing,  MSU  and  the  Air  Force  set  the  direction  of  the  program  beginning  with 
a  kickoff  meeting  in  February  of  2006.  At  that  time  the  electromagnetic  methods,  test  sample 
configuration  and  modeling  code  approaches  were  selected.  Task  1.1.1  Literature  Review  and 
Tasks  1.1.2  and  1.1.3,  Test  Samples  and  Model  Development,  are  discussed  in  sections  3.2, 

3.3  and  3.4  of  this  report.  Section  4  of  the  report  discusses  the  results  of  Tasks  1.2,  1.3  and 
1.4. 


Advanced  NDE  Sensor  Modeling  for  Multi-Site  Inspection  Program  Plan 


Figure  3  Program  Plan 
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Figure  4  Program  Schedule 


3.2  LITERATURE  REVIEW 

Computational  methods  for  modeling  electromagnetic  NDE  problems  fall  into  two  major 
classes,  namely,  analytical  and  numerical  approaches.  Analytical  models  are  capable  of 
providing  ideal,  closed  form  solutions  to  problems.  However,  such  analytical  models  are,  in 
general,  limited  in  their  application  to  simple  specimen  and  defect  geometries,  such  as  infinite 
half  plane  samples  and  regular,  rectangular  or  hemispherical  crack  geometry.  Numerical 
models  have  gained  popularity  in  recent  years  due  to  their  ability  to  simulate  arbitrary  shaped 
defects  in  complex  sample  geometries.  Numerical  models  can  simulate  problems  with  high 
dimensionality  as  well  as  nonlinear,  anisotropic  and  inhomogeneous  material  properties. 

With  the  development  of  faster  and  more  powerful  computers,  the  development  and 
implementation  of  the  numerical  models  have  become  a  major  research  focus  in  NDE115’161. 
Several  numerical  modeling  methods  have  been  developed  for  solving  the  governing 
equations  describing  various  NDE  phenomena,  such  as  finite  difference  method  (FDM), 

FEM,  boundary  element  methods  (BEM),  boundary  integral  methods  (BIM)  and  volume 
integral  methods  (VIM).  Each  method  has  its  own  advantages  and  limitations  and  is 
appropriate  for  different  kinds  of  numerical  problems  [10,17’18’19l 

The  finite  difference  method  is  the  simplest  numerical  modeling  method  for  solving 
partial  differential  equations  (PDE).  The  advantage  is  its  relative  ease  and  intuitive  method 
for  approximating  partial  derivatives  by  corresponding  finite  difference  formulation.  It  is 
widely  applied  for  direct  current,  quasi-static,  transient  fields,  and  linear  problems.  A 
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limitation  of  the  finite  difference  technique  is  its  poor  convergence  property,  particularly  in 
the  case  of  irregular  geometries.  Also,  finite  difference  formulations  do  not  lend  themselves 
easily  for  modeling  distributed  parameters  such  as  the  current  densities,  conductivities  and 
permeability'141. 

FEM  evolved  in  the  late  1950s  as  a  numerical  technique  in  structural  analysis  and  was 
quickly  adapted  as  a  major  numerical  modeling  method  in  various  engineering  fields  to  solve 
PDE.  The  finite  element  domain  discretization  scheme  allows  modeling  of  complex  shapes, 
making  it  efficient  and  relatively  more  accurate'101.  Compared  with  FDM,  FEM  has  many 
advantages  including  its  ease  of  imposing  essential  boundary  conditions  and  modeling 
complex  geometries.  FEM  formulations  can  be  easily  extended  to  handle  higher  order 
approximation  and  thereby  lead  to  faster  convergence  and  better  accuracy.  The  resulting 
matrix  equations  are  in  general  sparse,  banded  and  diagonally  dominant  making  the  solutions 
process  very  stable'121.  FE  models  have  been  developed  and  used  for  a  variety  of  two-  and 
three-dimensional  eddy  current  NDE  problems170,23]. 

The  drawbacks  of  FEM  include  large  computer  resources  especially  for  nonlinear  and 
time-dependent  problem.  Also,  the  finite  element  method  is  not  well  suited  for  open  region 
problems'121  encountered  in  wave  regimes.  In  the  area  of  antenna  and  electromagnetic  wave 
propagation,  integral  equation  based  approaches  are  more  commonly  used.  BEM  and  VIM 
methods  based  on  the  integral  equations'12,131  along  with  absorbing  boundary  conditions  are 
typically  used  in  wave  propagation  problems  where  the  Green’s  function  is  readily  available. 

Boundary  integral  methods  solve  Maxwell’s  equations  by  integrating  the  appropriate 
Green’s  function  over  the  boundary  using  the  given  boundary  conditions.  Hence  there  is  a 
need  to  discretize  only  the  boundary  surface  in  the  integral  equations,  rather  than  calculating 
all  values  in  solution  domain.  In  the  post-processing  phase,  the  integral  equation  is  used  again 
to  calculate  the  required  physical  quantity  at  any  point  in  the  solution  domain. 

In  the  case  of  inhomogeneous  problems,  VIM  is  commonly  applied  to  discretize  the 
volumetric  source  for  solution.  In  VIM,  the  field  is  determined  at  a  point  by  summing  the 
effects  of  the  sources  at  all  points.  The  volume  integral  computes  a  convolution  of  the  source 
points  with  the  appropriate  Green’s  function.  The  advantage  is  it  is  only  necessary  to 
construct  a  mesh  over  the  test  sample  and  solve  for  the  currents  in  the  test  sample.  But  this 
method  also  requires  computation  of  the  Green’s  function  which  in  itself  may  not  be  trivial. 

Integral  equation  based  formulation  is  in  general  more  efficient  than  differential  equation 
model  with  regard  to  computer  resources  but  the  resulting  global  stiffness  matrix  equation  is 
generally  fully  populated  and  ill  conditioned,  leading  to  instabilities'211.  More  importantly 
VIM  and  BEM  are  inherently  linear  and  can  be  used  to  solve  linear  homogeneous  problems. 
Lastly  these  methods  are  useful  only  when  the  Green’s  function  for  the  problem  is  easily 
available. 

In  this  project,  we  use  the  finite  element  method  for  modeling  the  multilayered  geometry 
largely  due  to  the  ease  of  modeling  ferritic  and  nonferritic  fasteners  and  complex  defect 
shapes,  as  well  as  the  inherent  stability  and  accuracy  of  the  solution  procedure. 
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3.3  TEST  SAMPLES 


In  this  program,  where  modeling  will  be  used  to  predict  and  optimize  performance,  test 
standards  play  a  vital  role.  The  test  samples  need  to  be  well  suited  for  model  development  yet 
representative  of  the  aging  aircraft  issues  facing  the  Air  Force.  Figure  4  lists  some  of  the  test 
sample  considered  for  this  program  and  Figure  5  through  10  are  photographs  of  the  various 
samples. 


Table  4  Representative  List  of  Boeing  Test  Samples  Options 


Sample 

Size 

Description 

Calibration 

Standard 

0.625-inch-thick  aluminum  plate  with  notch 

Leading  Edge 

Part 

Leading  edge  structure 

Leading  Edge 

Spar 

Leading  edge  spar 

Wing  Splice 
Standard 

20  x  13  x  1 
inch 

0.5  inch  aluminum  with  0.25  inch  double 
lap  splice 

Lower  Wing 
Structure 

Standards 

Four  samples:  S-l,  S-2,  S-3  and  S-4 
representative  of  707  type  center  wing  tank 
structure  inspection 

Lower  Structure 

12x4  inch 

W-67  configuration 
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Figure  6  Leading  Edge  Part 


Figure  7  Leading  Edge 


Figure  8  Wing  Splice  Standard 
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Figure  9  Lower  Wing  Structure  Standards 


Of  these  samples,  the  criteria  for  being  suitable  to  model  and  representative  of  a  critical 
Air  Force  need  focused  on  the  W-67  location  on  KC-135  (Boeing  707  type)  aircraft  center 
wing  box  structure  configuration  of  samples  S- 1 ,  S-2,  S-3  and  S-4  provided  by  Paul 
Rutherford121  \  Currently  standard  eddy  current  inspection  is  used  on  these  samples.  These 
standards  represent  a  splice  with  four  rows  of  fasteners  and  two  rows  on  each  side  of  the 
splice.  The  fasteners  are  titanium  on  one  side  of  the  splice  and  steel  on  the  other.  S-l  and  S-2 
have  0.16  inch  top  layers  while  samples  S-3  and  S-4  have  0.25  inch  top  layers.  Edge  views  of 
S-l  and  S-4  are  shown  in  Figure  10.  The  bottom  splice  plate  is  0.313  inch  thick.  Notches 
have  been  inserted  at  the  fastener  holes.  Figure  1 1  through  14  show  the  layout  and  section 
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cuts  for  the  four  samples.  Figure  15  shows  the  notch  defect  size  and  configuration  in  the 
samples.  The  notch  configurations  were  triangular  including  0. 10  inch  length  by  0. 13  inch 
depth,  0.13  inch  length  by  0.17  inch  depth,  0.15  inch  length  by  0.2  inch  depth,  0.17  inch 
length  by  0.22  inch  depth  for  S-l  and  S-3  and  0.2  inch  length  by  0.26  inch  depth,  0.22  inch 
length  by  0.29  inch  depth,  0.25  inch  length  by  0.313  inch  depth  and  0.30  inch  length  by  0.313 
inch  depth  for  S-2  and  S-4.  Table  5  shows  the  results  of  eddy  current  inspection  on  the  S-2 
and  S-4  samples.  The  signal  to  noise  is  interpreted  by  the  operator,  where  values  of  1  indicate 
that  the  defect  cannot  be  detected.  In  general  the  standard  eddy  current  inspection  is  not 
sensitive  in  the  second  layer  to  notches  less  than  0.25  inch  in  length. 


Figure  11  S-l  Drawing 


Figure  14  S-4  Drawing 
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Figure  15  Notch  Sizes  and  Configurations 


Table  5  Eddy  Current  Detection  of  Notches  in  Samples  S-2  and  S-4 


Signal-to-Hoise  (X:1) 

Sirin  Thickness 

F  ariaier  Material 

Fastener  Row 

Probe 

V-Gain 

H-Gain 

0.30-inch  Flaw 

0l2f- inch  Flaw 

0.22- inch  Flaw 

0.20-inchFlaw 

0.16 

steel 

outside 

SPC-181 

65 

50 

5 

2.5 

2 

2 

0.16 

steel 

outside 

SPC-223/1 00HZ 

65 

50 

5 

3 

2 

2 

QO 

0.16 

steel 

inside 

SPC-223/1 00H2 

65 

50 

3 

2 

1.5 

1.5 

oZ 

0.16 

titanium 

outside 

SPC-181 

75 

60 

3 

1.5 

1.5 

1.5 

0.16 

titanium 

outside 

SPC-223/1 00HZ 

75 

60 

3.5 

2 

2 

1.5 

0.16 

titanium 

inside 

SPC-223/1 00HZ 

75 

60 

2.5 

1.5 

1 

1 

.25 

steel 

outside 

SPC-181 

70 

55 

3 

1.5 

1 

1 

.25 

steel 

outside 

SPC-223/1 00HZ 

70 

55 

3.5 

2 

1 

1 

S4 

.25 

steel 

inside 

SPC-223/1 00HZ 

70 

55 

2.5 

1.5 

1 

1 

.25 

titanium 

outside 

SPC-181 

80 

65 

2 

1 

1 

1 

.25 

titanium 

outside 

SPC-223/1 00HZ 

80 

65 

3 

1.5 

1 

1 

.25 

titanium 

inside 

SPC-223/1 00HZ 

80 

65 

2 

1.5 

1 

1 

3.4  MODEL  DEVELOPMENT 
3.4.1  Introduction 

The  FEM  is  one  of  the  most  widely  used  tools  for  analysis.  The  method  finds  applications 
in  diverse  areas  of  engineering  and  science  largely  due  to  its  ability  to  model  intricate 
geometries  efficiently  and  accurately.  In  the  last  few  decades,  FEM  has  been  applied  with 
great  success  in  such  areas  as  the  study  of  DC  and  low  frequency  electromagnetic  fields  in 
electrical  machines,  design  of  structures,  modeling  thermal  fields  and  NDE.  The  FEM  was 
selected  for  this  project  to  provide  an  accurate  and  robust  simulation  of  the  electromagnetic 
fields. 

A  formulation  of  the  problem  from  first  principles  derived  from  Maxwell’s  equations 
governing  the  differential  equations  to  be  solved  in  the  model  is  provided  by  Jin[22]  and  shown 
in  Equations  (1)  to  (7).  Two-  and  three-dimensional  electromagnetic  finite  element  models 
for  electromagnetic  NDE  applications  have  been  developed  over  the  last  two  decades  and 
used  extensively  by  the  MSU  team.  Examples  of  successful  modeling  efforts  include  two- 
and  three-dimensional  eddy  current  models  for  simulating  steam  generator  tube 
inspections1-231,  three-dimensional  magneto-static  model  for  simulating  gas  transmission 
pipeline  inspection1741  and  three-dimensional  modeling  magneto-optic  imaging  inspection1 751 . 
More  recent  work  comprises  a  three-dimensional  model  for  microwave  inspection  of  re-bars 
in  concrete  structures1761.  The  major  steps  in  FEM  are  described  briefly  below. 

The  Maxwell’s  equation  in  differential  form  are  given  by 


SB 

(1) 

V  x  E  =  —  — — 

dt 

dD 

VxH  =  J+J  +  — 

(2) 

dt 

V  B  =  0 

(3) 
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V-D  =  p, 


(4) 


With  the  constitutive  relations: 

B  =  pH 

(5) 

tel 

ii 

(6) 

ii 

s 

(7) 

In  Equation  (1)  through  (7), 

E  =  electric  field  intensity  (volts/meter) 

D  =  electric  flux  density  (coulombs/square  meter) 

H  =  magnetic  filed  intensity  (amperes/meter) 

B  =  magnetic  flux  density  (webers/square  meter) 

Js  =  source  current  density  (amperes/square  meter) 

J  =  conduction  current  density  (amperes/square  meter) 
p  =  electric  charge  density  (coulombs/cubic  meter) 
p  =  permeability  (henrys/peter) 
s  =  permittivity  (farads/meter) 

<j  =  conductivity  (siemens/meter). 

The  vector  potential  formulation,  A-V Formulation  is  used  for  this  project.  From  (3),  the 
magnetic  flux  density  is  divergence  free  and  can  be  expressed  as  the  curl  of  a  vector,  i.e., 

B=VxA,  (8) 

where  A  is  the  magnetic  vector  potential.  The  magnetic  field  intensity  can  also  be  expressed 
in  terms  of  .4,  as 

H  =  —VxA  (9) 

Substituting  (8)  into  the  Faraday’s  Law  (1),  we  get 
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(10) 


Vx£  =  -V  x 


f 


Vx 


E  +  M 

V  &  ; 


8A 

~dt 


=  0 


(11) 


Then  the  electric  field  intensity  can  be  expressed  as 
E  =  -ja>A-VV 


(12) 


The  j<x>  term  is  introduced  since  we  are  dealing  with  the  time-harmonic  problems  for  this 
project.  Equation  (12)  tells  us  that  the  electric  field  intensity  E  can  be  expressed  in  terms  of 
the  magnetic  vector  potential  A  and  the  scalar  function  V. 


3.4.2  Geometry  and  Mesh  Generation 

The  test  geometry  truly  simulated  is  that  of  the  S2  sample  as  described  in  Section  3.3.  It 
consists  of  two  layers  of  aluminum.  The  top  layer  is  0. 16  inch  thick  and  the  bottom  layer  is 
0.313  inch  thick.  The  sample  contains  a  top  layer  and  second  layer  edges  close  to  the 
fastener.  Notches  are  introduced  at  the  fastener  holes  as  shown  in  Figure  15.  The  three- 
dimensional  domain  is  discretized  using  brick  or  hexahedral  elements  interconnected  at  eight 
nodes.  The  mesh  size  and  complexity  increases  particularly  when  a  small  air  gap  of  0.005 
inch  is  introduced  between  the  layers  and  also  between  the  fastener  and  plate.  The  density  of 
mesh  elements  were  optimized  and  validated  by  comparing  the  model  prediction  with 
experimental  measurements. 


3.4.3  Interpolation/Shape  functions  and  Global  Matrix  Assembly 

The  magnetic  vector  A  and  electric  scalar  potential  V  in  an  element  e  is  expressed  in 
terms  of  selected  linear  or  quadratic  interpolation  functions  as  follows 

8  24 


=  'ZN'iA'ii+NjA‘>, 


y  +  N-A-z 


=  2>. 


Ae 

k^k 


(13) 


j= i 


k= 1 


where 
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ve=ZN7; 

7=1 
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(14) 

'N‘x 

*  =  37-2 

Njy 

*  =  37-1  , 

(15) 

Njz 

k  =  3  j 
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and 


Ae 

*  =  37-2 

Ae 

Ayj 

*  =  37  - 1  , 

(16) 

Ae 

Azj 

*  =  37 

where  Axj ,  Ayj ,  Azj  and  V  .  are  the  three  Cartesian  components  of  the  vector  potential  and  the 
scalar  potential  at  node  j;  and  AC  is  the  shape  function  associated  with  node  j  that  assumes 
the  value  of  one  at  this  node  and  the  value  of  zero  at  any  other  node.  The  vectors  x,  y,  z  are 
the  Cartesian  unit  vectors. 

Using  the  Galerkin  formulation127-1,  and  imposing  the  Coulomb  gauge  on  (8)  and  (12),  we 
obtain 
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for  i  =  l,2,---,24 
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II  jaxffNl  ■  NjdV}Aj  +  XI  oVNek  -  VNjdVjVj  + 

j  i  e  2=1 

lao  ^  -W)-hdS  =  0 , 


(17a) 


(17b) 


for  k  =  1,2,  •  •  •  ,8  where  h  is  the  outward  unit  vector  and  dQe  is  the  boundary  of  the  sub- 
domain  Qe .  Combining  these  equations  together,  we  obtain  a  32x32  matrix  equation, 

[GY[UY=[QY.  (18) 

The  global  stiffness  matrix  is  obtained  by  assembling  each  element  matrix  [GY  together. 
The  value  at  each  entry  is  the  sum  of  values  contributed  by  all  the  connected  elements: 

GU  =  Q,  (19) 
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where  G  is  a  complex,  symmetric  sparse  matrix,  U  is  the  vector  of  unknowns  consisting  of  the 
electric  scalar  potential  and  the  three  components  of  the  magnetic  vector  potential  at  each 
node  and  Q  is  the  load  vector  incorporating  the  current  source. 

3.4.4  Boundary  and  Interface  Conditions 

To  obtain  a  unique  solution,  appropriate  boundary  conditions  need  to  be  imposed  before 
solving  the  system  of  equations.  Either  Dirichlet  boundary  conditions  (values  of  A  or  V  on 
the  boundary)  or  Neumann  boundary  conditions  (values  of  the  magnetic  flux  density  B  or  the 
magnetic  field  intensity  H  on  the  boundary)  need  to  be  specified.  Since  Neumann  boundary 
conditions  are  usually  included  implicitly  in  the  finite  element  formulation,  we  only  need  to 
impose  Dirichlet  boundary  conditions.  Current  continuity  conditions  will  be  explicitly 
imposed  at  the  interface  boundaries  to  avoid  spurious  solutions. 

3.4.5  Matrix  Solution 

The  finite  element  procedure  results  in  a  linear  algebraic  system  of  equations  that  must  be 
solved  to  determine  the  unknown  coefficients  of  the  shape  functions.  The  equations  can  be 
solved  using  either  direct  or  iterative  methods.  The  Gaussian  elimination  method  is  a  direct 
approach  that  can  be  used  to  solve  either  full  or  sparse  matrix  equations.  In  the  finite  element 
method,  the  stiffness  matrix  is  usually  sparse  and  banded.  In  such  cases,  it  is  more 
advantageous  to  employ  iterative  methods  since  they  are  more  efficient  with  respect  to 
computational  cost  and  data  storage.  A  number  of  iterative  methods  have  been  proposed  over 
the  years.  However,  only  a  small  number  of  them  can  be  used  for  solving  the  complex-valued 
matrix  equations  arising  in  quasi-static  electromagnetic  problems.  We  propose  to  use  the 
Transpose  Free  Quasi-Minimal-Residual  (TFQMR)  method  [28]  where  the  associated  Krylov 
subspace  is  Kn  =  span{Q,GQ,G2Q,  --,Gn~lQ} .  The  exact  solution  U  =  G  lQ  could  be 
approximated  by  Un  e  Kn  that  minimizes  the  norm  of  the  residual  rn=Q-  GUn .  The  TFQMR 
algorithm  relaxes  this  minimization  requirement.  It  aims  to  minimize  a  different,  data- 
dependent  norm  that  in  practice  is  similar  to  || rn  || . 

An  important  issue  in  iterative  methods  is  its  convergence  properties.  The  convergence 
of  iterative  process  depends  on  the  properties  of  the  matrix,  such  as  its  eigen  values,  the 
singular  value  or  condition  number.  Usually  a  large  condition  number  will  result  in  a  large 
solution  error  or  sometimes  even  failure  to  converge.  Thus,  the  process  of  preconditioning  is 
essential  for  the  success  of  iterative  methods.  The  overall  solution  procedure  requires  O(N) 
multiplications  and  additions1291. 

3.4.6  Post-Processing 

After  we  solve  for  the  potentials,  physical  and  measurable  quantities  of  interest,  such  as 
the  magnetic  flux  density  B,  coil  impedance  and  induced  current  density  J,  can  be  calculated 
using: 


B  =  VxA 

(20) 

E  =  -jcoA-VV 

(21) 
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./  =  (7  E  . 


(22) 


where  Equations  (20)  through  (22)  are  the  same  as  Equations  (8),  (7)  and  (12)  and 
presented  here  again  for  the  reader’s  convenience. 


GMR  sensors  measure  the  normal  component  of  the  magnetic  flux  density  B.  The 
different  components  of  the  magnetic  flux  density  can  be  obtained  from  the  nodal  values  of 
magnetic  vector  potential  as: 


ay  _<K 

dy  dz 


(23a) 


J  dz  dx 


(23b) 


dx  dy 


(23c) 


The  current  density  J  inside  element  e  is  determined  at  the  element  center  from  the 
relation 
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The  components  Jey  and  Jez  can  be  calculated  similarly. 

3.4.7  Modeling  GMR  Probes  with  Induction  Foil  Excitation 

The  excitation  source  in  this  project  is  applied  to  an  infinite  planar  induction  foil  located 
parallel  to  the  surface  of  the  specimen.  For  a  linear  sinusoidal  current  excitation  with  angular 
frequency  of  co ,  the  source  current  is 


j iin  =  Re[yoe-"  ] =  Jo  cos  cot ,  (25) 

where  Jo  is  the  magnitude  of  the  current.  However  a  probe  with  a  linear  excitation 
current  flow  pattern  is  sensitive  only  to  cracks  that  are  perpendicular  to  the  current  direction. 
In  this  project,  to  best  detect  a  crack  oriented  in  the  x  direction  the  current  flow  is  assumed  to 

be  along  the  y  direction,  i.e.,  ^,in  ~  where  ^  is  the  unit  vector  along  they-axis.  The 

D 

normal  component  of  magnetic  flux  density  z  will  be  measured  by  the  GMR  sensor  array. 
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4.  RESULTS  AND  DISCUSSIONS 


4.1  MR  SENSORS  MEASUREMENTS 

MR  sensor  measurement  and  data  analysis  have  been  of  major  focus  throughout  Boeing’s 
development  history  of  this  technology1-3,4,51.  Although  MR  sensing  techniques  are  similar  to 
and  are  based  upon  eddy  current  measurement  methods,  the  MR  sensor  measurement  and  data 
processing  methods  have  evolved  to  be  distinct  and  unique  from  traditional  eddy  current 
technology.  This  distinction  is  enough  to  classify  MR  sensors  as  a  new  modality  of 
electromagnetic  NDE  with  unique  advantages  over  eddy  current  for  particular  measurement 
needs.  This  uniqueness,  however,  requires  the  establishment  of  new  calibration  and 
measurement  techniques  as  well  as  new  data  interpretation  methodologies  and  advanced 
feature  extraction  algorithms. 

In  support  of  this  program,  measurements  were  taken  at  both  the  Boeing  Phantom  Works 
Physics  NDE  Laboratory  in  Seattle  using  the  MR  Sensor  Test  Bed  Scan  System  and  the 
Boeing  Phantom  Works  MR  Sensor/MAUS  Scan  System  in  St.  Louis.  The  Seattle-based 
measurements  were  used  initially  to  validate  and  calibrate  MSU  modeling  results.  Following 
this,  the  measurements  were  used  in  both  the  development  and  testing  of  advanced  feature 
recognition  and  signal  to  noise  (SNR)  extraction  algorithms  at  both  Boeing  and  MSU. 

Section  4.1  is  organized  as  follows:  Section  4.1.1  details  the  MR  sensor  test  systems 
developed  by  Boeing  and  used  in  this  program.  Section  4.1.2  presents  the  MR  test  sample 
data  in  support  of  the  modeling  and  feature  extraction  efforts.  Section  4.1.3  presents  the 
algorithms  developed  for  and  the  results  of  SNR  extraction  from  the  measurement  data. 

4.1.1  MR  Sensor  Test  System  and  Measurement  Process 

This  section  presents  the  MR  Sensor  test  systems  and  measurement  processes.  The 
existing  MR  systems  at  Boeing  Phantom  Works  include  a  prototype  Test  Bed  Scan  System 
unit  and  the  MAUS-based  system  presently  under  development.  Both  are  capable  of 
accomplishing  this  task  of  producing  measurement  data  in  support  of  the  modeling  and  data 
analysis  efforts  of  this  program.  The  laboratory-based  unit  was  used  for  taking  most  of  the 
experimental  test  data.  The  measurement  process  for  the  Test  Bed  Scan  System  is  detailed 
below. 

4. 1 . 1 . 1  Test  Bed  Scan  System 


The  Test  Bed  Scan  System  (shown  in  Figure  16)  consists  of  the  following  major 
components:  (1)  MR  Sensor  Head,  (2)  PXI  Chassis,  (3)  XY  Scan  Bridge,  and  (4)  Executive 
Computer. 
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Boeing  PW  Physics  Lab 
MR  Sensor  Scan  System 


MR  Sensors  Housed  in  a  “Can'1 
Attached  to  a  "Shoe” 
[Excitation  Current  Sheet) 


Magneto -Res  I  stive  Sensors  (GIVER)  Card 
32  Sensor  Array  - 1"  Width  -  “Z-Direction 

Sensor  array  under  black  epoxy 


Figure  16  Laboratory  MR  Sensor  System 


The  MR  Sensor  Head  consists  of  one  MR  sensor  card  with  a  linear  array  of  32  GMR 
sensors  arranged  in  a  row  1  inch  in  length  (see  bottom-right  picture  in  Figure  16).  The 
sensors  (and  card)  are  aligned  normal  to  the  scan/index  plane  (the  z-direction)  such  that  a  1 
inch  swath  is  swept  by  the  sensors  along  the  scan  direction  (the  x-direction).  The  card 
includes  cable  connectors  for  the  sensors.  The  MR  Sensor  Head  is  housed  in  a  metal  box  for 
shielding  (see  upper-right  picture  in  Figure  16).  At  the  bottom  of  the  head  is  a  “shoe”  to 
which  a  copper  foil  is  adhered.  This  foil  contains  strips  that  deliver  a  linear,  unidirectional 
excitation  current.  The  direction  of  current  is  along  the  line  of  sensors,  which  is  along  the 
index  direction  (the  y-direction).  The  Sensor  Head  also  contains  hardware  for  a  DC  strip 
current.  This  strip  current  is  provided  to  shift  the  operating  point  of  the  MR  sensors  such  that 
they  are  optimized  for  maximum  sensitivity  along  the  most  linear  portion  of  their  magnetic 
field  to  current  response  curves. 


The  PXI  Chassis  contains  both  Boeing-developed  hardware  for  excitation  signal  source, 
strap  current  and  encoder  signal  processing  and  off-the-shelf  hardware  components  for  both 
D/A  and  A/D  conversion  of  sensor  signals  and  digital  control  of  various  system  components. 
The  latter  tasks  are  accomplished  through  five  multifunction  cards  from  National  Instruments. 
A  voltage  output  channel  (D/A  converter)  of  one  card  is  used  to  “seed”  the  AC  excitation 
signal  for  desired  amplitude  and  frequency.  This  excitation  signal  is  pre-processed  in  the 
Boeing  hardware,  which  is  housed  in  one  slot  of  the  PXI  chassis  and  uses  the  PXI  power 
supply  bus.  The  excitation  signal  is  connected  through  cable  to  the  excitation  sheet  at  the 
bottom  of  the  MR  Sensor  Head.  Thirty-two  of  the  available  voltage  input  channels  (A/D 
converters)  are  used  to  input  the  return  signals  from  the  sensors,  which  are  first  pre-processed 
through  amplifiers  on  the  Boeing  hardware.  In  addition  to  these  32  signals,  two  input 
channels  are  devoted  to  reading  and  digitizing  the  excitation  signal  and  its  quadrature  signal 
in  the  same  manner  as  the  sensor  signals. 
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The  XY  Scan  Bridge  (upper  left  picture  of  Figure  16)  allows  for  accurate  position  and 
motion  control  within  0.01  inch  (and  theoretically  even  smaller)  in  both  in-plane  directions. 
The  position  and  scan  speed  of  the  MR  sensors  can  therefore  be  controlled  to  adequate 
precision  to  produce  high  resolution  c-scans. 

The  Executive  Computer  controls  the  PXI  Chassis  through  a  connected  PC-card  /  PXI- 
card  interface.  Control  software  is  written  in  LabVIEW  to  provide  extensive  graphical 
interface  and  control.  C-scan  data  is  saved  in  a  specialized  scan  data  file  format  (.SCD- 
extension).  A  converter  program  has  been  developed  to  read  and  parse  the  full  SCD  data 
format  into  the  Matlab  programming  environment.  From  Matlab,  the  data  can  be  manipulated 
or  stored  in  any  number  of  standard  formats. 

4. 1.1.2  MR/MAUS  System 


Some  data  was  also  taken  using  the  MR  /  MAUS  system  under  development  at  Boeing 
Phantom  Works  St.  Louis.  Figure  17  shows  an  early  proof-of-concept  arrangement  showing 
(a)  the  MR  Sensor  Head  mounted  on  a  MAUS  V  with  the  Test  Bed  System  in  the  background 
controlling  the  sensor  and  data  acquisition  portion  of  the  test  and  (b)  the  MAUS  V  control 
hardware  controlling  the  motion.  This  was  not  a  full  concept  integration.  The  new 
MR/MAUS  system  is  shown  in  Figure  18. 


(a) 


(b) 


Figure  17  Photographs  of  MR  Sensors  Attached  to  the  MAUS  V  System 

Note:  (a)  MR  sensors  being  swept  over  a  test  sample,  (b)  The  MAUS  V  data  acquisition 
system 
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Figure  18  Photographs  of  the  New  MR/MAUS  System  Under  Development  at  Boeing  St. 

Louis 


4. 1.1.3  Demodulation  Process 


It  is  instructive  to  examine  the  measurement  data  processing  method  before  examining 
and  interpreting  the  data.  The  following  method  describes  the  one  used  on  the  Test  Bed 
System,  which  is  the  system  from  which  most  of  the  MR  sensor  data  used  in  this  program  was 
obtained.  The  amplitude  and  frequency  of  the  AC  excitation  signal  is  selected  and  a  scan  is 
made.  The  32  sensors’  full  waveform  signals  are  digitized  along  with  the  excitation  signal 
and  its  quadrature.  Figure  19  shows  example  waveforms.  The  excitation  signal  and  its 
quadrature  have  the  following  time-dependent  waveforms: 

S{]{t)=  Aosm{lx ft) 

S9 0  (t)=  Aa  sin(2/r ft  +  ;r/2)  =  A0  cos(2;r ft) 

In  response,  the  nth  sensor  input  signal  has  the  following  form  (see  Figure  19): 

Sf>(/)=4,sin(2<r/r  +  A).  (27) 


Each  sensor  is  characterized  with  its  own  amplitude  An  and  phase  shift  <f>n  with  respect  to 

the  excitation  signal  SQ.  These  parameters  are  ultimately  functions  of  phase,  geometry, 
electronics,  etc. 

Through  experience  it  was  determined  that  the  amplitude  and  phase  themselves  were  not 
all  that  useful  in  producing  meaningful  c-scans.  Rather,  a  demodulation  process  was 
developed  to  express  the  results  in  two  related  orthogonal  components  that  could  more  readily 
be  manipulated  in  an  interactive  fashion.  The  in-phase  and  quadrature  demodulation 
functions  can  be  written  as  follows: 


vi  i 

\dt  =-A0Ancos{<j)n) 
)  2 

vr  i 

\dt  =-4,4,  sin  (<0 


0 

/NT 


(28) 
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Eh  Sty  WF  Graph 


Figure  19  Plots  from  Boeing  Seattle  Test  Bed  System 

Note:  (a)  in-phase  (white)  and  quadrature  (red)  excitation  signal,  (b)  32  sensor  signal 
responses  and  (c)  down-sampled  XY  Scan  Bridge  motor  pulses 


In  Equation  (28)  the  integration  is  done  over  a  multiple  of  the  oscillation  period  T  =  1//  . 
This  span  of  time  may  be  considered  a  moving  window  along  the  time  domain  of  the  signals. 
Down-sampled  motor  pulses  from  the  XY  Scan  Bridge  (see  Figure  19(c))  are  used  to 
determine  the  points  in  time  to  apply  the  demodulation  algorithm  of  Equation  (28).  These 
two  components  are  functions  of  the  amplitude  and  phase,  which  themselves  are  easily 
obtained  as  follows: 


■i=^-VY+E 

<t>n  =arctan  (T„/yJ 


(29) 


4. 1 . 1 .4  Detection  Angle 

The  demodulation  for  in-phase  (Xn)  and  quadrature  ( Yn)  were  described  in  Equation  (28). 
It  was  determined  that  demodulation  at  some  detection  angle  #  yielded  more  favorable  results 
in  which  the  response  from  grosser  geometrical  features,  such  as  fastener  holes,  was 
minimized  and  those  of  cracks  or  voids  were  enhanced.  In  earlier-developed  systems  this 
detection  angle  was  selected  through  hardware,  and  additional  scans  had  to  be  made  whenever 
that  angle  was  changed.  One  of  the  main  driving  forces  for  developing  the  Test  Bed  System 
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was  to  generate  the  two  component  demodulated  signals  through  software  in  order  to 
minimize  the  burden  on  the  electronics  and  allow  for  software  selection  of  detection  angle 
from  a  single  scan. 

The  nth  sensor  signal  is  considered  along  with  an  artificially-introduced  detection  angle  &. 

S'';)(t)=A.sm(2nft  +  t,-d).  (30) 

Equation  (27)  is  simply  Equation  (30)  at  detection  angle  6  =  0 .  If  the  same  in-phase 
demodulation  algorithm  of  Equation  (28)  is  applied  to  Equation  (30),  one  obtains  the 
following  expression: 


NT  /NT  , 

F,(6)=  \dts,(t)$\t)  jd,  =  -A,A,o  os(t,-0) 

=  \AoAn  cos(^)cos(6»)+|  An  sin(^)sin(6»).  (31) 

=  Xn  cos(#)+  Yn  sin(#) 

In  the  final  expression  of  Equation  (31),  one  can  see  that  this  demodulated  signal  with 
arbitrary  detection  angle  reduces  to  a  simple  algebraic  sum  of  the  in-phase  and  quadrature 
components.  Equation  (31)  demonstrates  that  it  is  sufficient  for  each  sensor  to  calculate  and 
store  only  the  two  orthogonal  demodulated  components.  From  these,  the  amplitude  and  phase 
of  the  original  time-dependent  waveform  can  be  reconstructed  from  Equation  (29),  and  these 
two  components  can  also  be  used  to  construct  a  demodulated  signal  at  any  desired  detection 
angle.  Indeed,  the  interactive  manipulation  of  this  detection  angle  and  the  resulting  change  in 
c-scan  has  become  a  vital  tool  for  examining  MR  sensor  results. 

4. 1.1.5  Calibration 


The  calibration  of  the  MR  sensors  can  now  be  better  described  following  this  explanation 
of  demodulation  function  and  detection  angle.  Each  sensor  has  its  own  unique  response  curve 
to  magnetic  field.  As  mentioned  above,  the  operating  point  of  the  sensors  can  be  shifted 
together  using  the  strap  current,  however  each  sensor  still  has  its  own  unique  dependency.  A 
common  zero  in  the  in-phase  and  quadrature  components  is  first  obtained  by  scanning  at  a 
stationary  position  on  a  thick  sample  devoid  of  any  features  and  far  enough  away  from  any 
edges.  This  is  referred  to  as  Calibration  1.  Figure  5  shows  a  full  picture  of  the  calibration 
standard,  and  Figure  16  shows  a  close-up  picture  of  the  standard  with  the  sensor  head 
positioned  for  taking  the  zero  calibration  measurements.  The  two  demodulated  components 
for  each  sensor  that  are  measured  at  this  zero  location  are  recorded  and  are  then  subtracted 
from  all  subsequent  scans. 

The  effects  of  applying  Calibration  1  can  be  seen  in  Figure  20(a)  and  Figure  20(b). 

These  figures  show  the  in-phase  (0=  0°)  demodulated  signal  from  scans  across  the  calibration 
standard  before  and  after  zero-calibration  (Calibration  1),  respectively.  It  is  clear  to  see  that 
the  low-signal  portion  of  the  curves,  which  correspond  to  areas  of  the  sample  farther  from  the 
crack  and  edges,  are  brought  tighter  together  following  this  calibration  step. 
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Figure  20  Plots  of  In-phase  Demodulated  Signal  across  the  Calibration  Standard  of  Figure 
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Note:  (a)  uncalibrated,  (b)  zero  calibration  and  (c)  gain  calibration 


A  gain  factor  (Calibration  2)  is  then  determined  from  the  measurements  shown  in  Figure 
20(b)  and  applied  to  each  sensor.  The  results  of  applying  both  Calibrations  1  and  2  are  shown 
in  the  tightly  grouped  plot  of  Figure  20(c).  Calibration  2  is  determined  as  follows. 
Examination  of  the  measurements  across  this  crack  show  an  “s”-shape  response,  as  seen  in 
Figure  20,  which  is  dependent  upon  the  detection  angle  as  described  in  Equation  (31).  The 
peak-to-peak  or  area  under  the  curves  is  then  used  to  determine  a  gain  factor  for  each  sensor. 
In  practice,  one  sensor,  typically  one  whose  response  lies  in  the  middle  of  the  group  of  32 
responses,  is  chosen  to  have  a  gain  factor  of  1,  and  the  others  are  assigned  a  gain  factor 
accordingly.  Again,  the  final  result  of  zero  subtraction  (Calibration  1)  and  gain  factor 
(Calibration  2)  are  shown  in  Figure  20(c),  where  the  32  curves  fit  much  more  tightly  together 
and  give  more  similar  responses  across  the  calibration  sample. 

Mathematically,  the  fully  calibrated  demodulation  signal  for  the  nth  sensor  has  the 
following  form: 
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f,'(0)=*>s(0)+r;sin(0) 

K=c.(X,-a ,) 

K=c.(r,-b,) 


(32) 


In  Equation  (32),  the  primed  quantities  denote  calibrated  values,  a„  is  the  in-phase  zero 
calibration,  bn  is  the  quadrature  zero  calibration,  and  cn  is  the  gain  calibration.  X„  and  Y„  are 
the  uncalibrated  in-phase  and  quadrature  demodulation  results,  respectively,  from  Equation 
(28).  Establishment  of  the  three  calibration  coefficients  is  done  using  the  procedures 
described  above.  The  resulting  calibrated  demodulation  results  can  be  displayed  effectively 
as  a  c-scan  with  arbitrary  adjustable  detection  angle  6.  Optimization  of  this  detection  angle 
has  to  date  been  accomplished  interactively  by  the  operator.  The  operator  observes  the 
changing  patterns  of  the  c-scan  as  he  adjusts  6.  At  optimum  6,  patterns  and  amplitudes  due  to 
the  presence  of  the  feature  of  interest  will  be  best  highlighted  in  comparison  with  other 
portions  of  the  c-scan  that  are  devoid  of  the  feature.  For  example,  scans  in  search  of  cracks 
around  fasteners  should  yield  differences  between  flawed  and  normal  fasteners  of  like  size, 
shape  and  material  composition  when  the  detection  angle  is  properly  tuned  and  the  cracks  are 
large  enough  to  highlight  such  differences. 

4.1.2  Test  Sample  MR  Data 

To  date  the  extraction  of  features  of  interest  from  c-scans  has  been  accomplished 
subjectively  by  an  operator  who  requires  much  experience  and  training  in  establishing  a 
method  of  adjusting  the  detection  angle,  c-scan  contrast,  etc.  This  section  highlights  a 
comprehensive  data  set  of  results  that  were  measured  in  April  2007  from  the  lower  wing 
structure  standards  of  Figure  9.  These  measurements  were  initially  taken  to  assist  in 
validation  of  the  electromagnetic  FEM  calculations,  presented  in  Section  4.2.1.  This  set  was 
also  used  as  in  the  development  of  a  more  objective,  quantitative  set  of  algorithms  for  feature 
extraction  and  establishment  of  signal  to  noise  values,  which  are  presented  in  Section  4.1.3. 
An  even  more  extensive  set  of  measurements  on  the  same  samples  were  taken  in  June  2008  in 
part  as  optimized  scans  following  the  FEM  calculation  results,  which  are  presented  in  Section 
4.3. 

4. 1.2.1  Subjective  Data  Analysis  Tools  and  File  Conversion  Procedures 

The  measurement  results  were  analyzed  using  a  set  of  post-processing  tools  developed  at 
Boeing  Phantom  Works  Seattle.  As  mentioned  in  Section  4. 1.1.1  data  was  recorded  in  a 
custom  scan  data  binary  format  (files  with  extension  “SCD”).  In  this  format,  the  data  and 
nearly  all  control  settings  of  the  data  acquisition  software  are  recorded.  It  was  convenient  to 
extract  a  reduced  set  of  data  from  this  format  into  a  structure  that  was  more  universal  for  post¬ 
processing  and  distribution  to  others.  Matlab  functions  were  developed  to  read  and  parse  the 
scan  data  from  SCD  files  into  the  Matlab  environment.  Details  of  these  functions  are 
presented  in  Appendix  A.  The  reduced  pertinent  data  could  then  be  processed,  visualized  and 
further  analyzed  in  the  Matlab  environment  as  desired.  In  addition,  extracted  data  could  be 
saved  in  MAT  files.  This  is  a  standard  data  format  that  is  recognized  by  many  programming 
environments,  not  just  Matlab. 
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It  was  desired  to  build  an  interactive  post-processing  tool  in  a  programming  environment 
that  had  a  wide  range  of  visualization  and  processing  tools.  Mathematica  was  chosen  because 
of  a  recent  major  upgrade  with  Version  6  that  contained  even  more  advanced  graphics  as  well 
as  new  interactive  manipulation  tools.  Within  this  environment,  a  panel  of  tools  was  created 
to  load  scan  data  from  MAT  fdes  and  manipulate  the  presentation  of  data  to  best  optimize 
(subjectively)  the  c-scans  to  highlight  features  of  interest.  This  Mathematica  notebook  was 
developed  as  a  platform  on  which  new  post-processing  and  visualization  algorithms  could  be 
applied  and  tested.  Most  of  the  c-scan  results  presented  in  this  report  were  generated  by  this 
tool.  Details  of  the  Mathematica  notebook  tool  are  given  in  Appendix  B. 

4. 1.2.2  C-Scans  of  S-2  Test  Standard:  Titanium  and  Steel  Sides 


Scans  were  made  of  the  S-l,  S-2,  S-3  and  S-4  samples  at  frequencies  of  100,  150,  200, 
250,  300,  350,  400,  450,  and  500  Hz.  Figure  21  shows  the  400  Hz  results.  Figure  21(a) 
shows  a  schematic  of  the  S-2  sample,  which  is  the  standard  used  through  most  of  the 
modeling  and  feature  extraction  work.  The  scan  direction  was  along  x  and  index  direction 
along  y  in  this  diagram  with  the  excitation  current  in  the  y  direction.  In  this  orientation,  the 
steel  side  is  on  the  left,  the  titanium  side  on  the  right.  The  various  fasteners  are  labeled  here 
as  SI  through  S21  for  the  steel  side,  odd-numbered  on  the  outside,  even-numbered  on  the 
inside.  T1  through  T21  are  similarly-labeled  for  the  titanium  side.  The  schematic  also  shows 
the  location  and  approximate  size  of  the  various  cracks,  which  are  detailed  in  Figure  15.  The 
crack  length  shown  here  is  the  top  dimension  of  the  notch  cracks  in  inches.  Only  the  largest 
two  crack  configurations  are  through  cracks.  The  geometry  and  fastener  numbering  scheme 
of  Figure  21  should  be  used  as  reference  through  the  rest  of  the  discussion  on  the 
measurements  on  the  S-2  sample. 
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Figure  21  Plots  of  S-2  Scan 

Note:  (a)  geometry,  crack  dimensions  in  inches  and  fastener  numbering,  (b)  full-body  in- 
phase  0=0°  and  (c)  quadrature  0=  90°  c-scans 


Figure  21(b)  and  (c)  show  the  full-body  in-phase  and  quadrature  c-scans,  respectively,  of 
the  S-2  standard  taken  at  a  frequency  of '/=  400  [Hz],  The  c-scans  will  be  analyzed  in  more 
detail  below.  One  can  note  the  larger  signal  on  the  steel  side  vs.  the  titanium  side. 
Furthermore,  it  so  happens  that  the  quadrature  c-scan,  or  7-component,  corresponding  to  a 
detection  angle  of  0=  90°,  is  close  to  optimum  angle  for  enhancing  the  crack  signatures.  This 
is  just  serendipitous  as  it  is  not  generally  the  case  for  other  frequencies  or  other  samples 
scanned.  The  largest  steel  cracks  at  S14  and  S 15  are  easily  seen  even  on  these  wide-area 
plots.  Note  that  the  features  at  S19  and  T19  are  surface  cracks,  which  were  not  of  interest  to 
this  work. 


As  stated  above,  the  optimum  detection  angle  is  discovered  somewhat  subjectively  by  the 
operator  through  interactive  adjustment  of  the  detection  angle,  data  contrast,  color  palette,  etc. 
In  practice,  a  particular  feature  of  the  data,  perhaps  the  magnitude  and/or  shape  of  a  particular 
lobe  surrounding  the  fasteners,  is  examined  comparatively  from  one  fastener  to  another  in 
search  of  differences.  When  these  differences  are  pronounced  enough,  the  operator  may 
designate  that  a  flaw  has  been  visualized. 

An  example  of  this  analysis  of  specific  c-scan  features  is  shown  in  Figure  22.  Here  only 
a  section  of  the  titanium  side  has  been  enhanced  to  illustrate  this  point.  In  Figure  22(b),  the 
results  at  0=  69.84°  are  shown  with  the  contrast  scale  squeezed  to  enhance  the  smaller  crack. 


32 


The  larger  0.3  in  cracks  at  T14  and  T15  are  quite  obvious  and  have  been  saturated  here. 
Observing  the  inside  row  (even-numbered  fasteners),  one  can  see  that  the  blue  lobe  to  the 
right  of  the  fastener,  the  crack  side,  has  progressively  larger  magnitude  and  size  as  seen  at  T4, 
T8  and  T12  than  the  non-crack  ones  at  T6  and  T10.  The  operator  observed  this  lobe  as  the 
detection  angle  was  adjusted  and  determined,  albeit  subjectively,  that  the  differences  in  lobes 
were  best  enhanced  at  this  particular  angle. 
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Figure  22  Analysis  of  the  S-2  Titanium  Side  Cracks  at /=  400  Hz  and  0=  69.84°. 

Note:  The  plots  are  of  the  (a)  area  of  interest,  (b)  c-scan  data,  (c)  closeup  of  T1 1  and  T12 
cracks  with  bottom  layer  edge  signature. 

For  the  outside  row  (odd-numbered  fasteners),  the  operation  was  a  bit  more  difficult  due 
to  the  proximity  to  the  edge  of  the  bottom  layer  (dashed  line  in  the  plots).  The  response  from 
this  edge  can  suppress  crack  signal  and  identification  if  the  fastener  is  too  close  to  the  edge. 

In  this  case,  it  is  apparent  in  Figure  22(b)  that  the  mirror  lobe  to  the  left  of  the  fasteners, 
opposite  side  from  the  cracks,  can  be  observed  and  optimized  for  best  detection.  The  red 
lobes  of  T3,  T7  and  Til  are  progressively  larger  in  magnitude  and  size  compared  with  the 
non-crack  ones  at  T5,  T9  and  T13.  Optimizing  on  the  crack-side  lobe  of  the  outside  row 
fasteners  is  not  as  successful.  One  can  see  that  the  right-side  lobe  of  T1 1  is  notably  larger,  but 
this  discernment  is  harder  to  make  for  T7  and  T3  when  compared  with  the  non-crack  lobes. 
Figure  22(c)  shows  an  enhancement  of  the  edge  effect  on  the  data.  The  signal  at  the  edge  is 
obviously  saturated  and  tends  to  suppress  the  right  lobe. 

On  the  steel  side,  the  discernment  of  cracks  against  the  large  response  of  steel  fasteners  is 
even  more  challenging.  Figure  23  shows  the  c-scan  of  the  steel  fastener  section  of  the  S-2 
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sample  at /=  400  [Hz]  and  0=  89.6.  This  c-scan  is  actually  of  the  absolute  value  of  the 
demodulated  signal,  which  mathematically  is  |F'|  in  Equation  (32),  and  the  plot  was  made 

with  a  different  color  palette  than  that  of  Figure  22.  The  0.3  inch  cracks  are  easily  detected 
and  are  saturated  on  this  scale.  For  the  steel  fasteners,  the  observer  must  analyze  the  shape 
and  magnitude  of  the  halo  on  this  plot.  For  the  outside  row,  one  can  easily  see  that  the 
magnitudes  of  SI  1  and  S7  are  notably  larger  than  the  others.  The  halo  of  S3,  although  a 
weaker  signal,  does  demonstrate  the  same  connection  with  the  small  right-side  lobes  of  the 
other  fasteners  with  crack.  However,  S3  is  admittedly  more  difficult  to  discern  from  the 
fasteners  with  no  crack.  For  the  inside  row,  the  halos  of  the  fasteners  with  crack  are  notably 
more  complete  and  connected  in  shape  than  those  without  cracks.  Interestingly,  S8  is  more 
difficult  to  see  than  S4,  despite  being  from  a  slightly  larger  notch  crack.  The  discovery  of 
cracks  on  the  steel  side  using  this  somewhat  subjective  methodology  is  indeed  problematic, 
which  is  of  course  the  motivation  for  developing  more  precise  quantitative  algorithms  to 
apply  to  these  kinds  of  cases. 
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Figure  23  Analysis  of  the  S-2  Steel  Side  Cracks  at  f  =  400  Hz  and  q  =  89.6°. 
Note:  The  plots  are  of  the  (a)  area  of  interest,  (b)  c-scan  data. 

4. 1.2.3  C-Scans  of  S-2  Test  Standard:  Frequency  Variation 


There  is  another  dimension  in  which  to  explore  the  data,  namely  the  frequency  of  scan. 

In  general,  lower  frequency  scans  highlight  deeper  features,  however  it  is  the  nature  of  all 
magnetic  detection  techniques  that  deeper  features  are  more  diffused  than  shallower  features 
that  are  more  sharply  detected  at  higher  frequencies.  This  property  is  related  to  the  skin  depth 
of  the  material. 
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It  stands  to  reason,  therefore,  that  the  larger  cracks  on  the  S-2  sample  would  appear 
prominently  over  a  wider  range  of  frequencies  since  they  are  through  cracks  and  would  be 
enhanced  at  all  scan  depths  through  the  second  layer.  For  the  comer  notch  cracks,  however, 
this  should  not  be  the  case  as  these  cracks  lie  within  a  narrow  range  of  depth  close  to  the 
interface  between  the  top  and  bottom  layers.  However,  it  is  important  to  understand  that  the 
fields  measured  by  the  sensors  are  themselves  detected  right  at  the  immediate  environment  of 
the  sensors,  which  is  at  some  elevation  above  the  top  of  the  sample,  indeed  just  above  the 
excitation  sheet.  These  fields  are  themselves  partly  induced  by  eddy  currents  present  at  all 
depths  through  the  material.  In  other  words,  the  final  magnetic  field  measured  by  the  sensors 
could  be  thought  of  as  some  complex  integration  of  all  currents  through  the  whole  the  depth 
of  material.  The  weighting  factor  of  one  depth  to  another  is  considered  then  a  strong  function 
of  the  frequency. 

Having  gone  through  this  somewhat  heuristic  explanation,  it  is  illustrative  to  compare  the 
c-scans  of  the  S-2  sample  at  the  various  frequencies.  Figure  24  shows  the  inside-row 
titanium-side  results  for  frequencies  from /=  100  to  500  Hz  in  steps  of  50  Hz.  The  outside- 
row  has  been  masked  out  to  aid  in  viewing  the  inside-row  only.  Again,  the  optimum 
detection  angle  for  each  frequency  was  determined  subjectively  by  observing  the  relative 
contrast,  scale  and  shape  of  the  lobe  to  the  right  (crack-side)  of  the  fasteners.  This  optimized 
detection  angle  generally  increases  with  frequency,  and  the  changes  are  smaller  between 
higher  frequencies.  The  contrast  was  further  reduced  to  enhance  the  smaller  cracks  since  the 
larger  ones  were  obvious  at  all  frequencies.  Comparison  between  the  smaller  crack  fasteners 
(T04  and  T08)  and  the  no  crack  fasteners  (T06  and  T10)  shows  that  there  really  is  some 
difference  at  all  frequencies,  however  this  difference  is  best  enhanced  above,  say  200  Hz.  400 
Hz  has  been  stated  as  the  best  frequency  for  this  row  on  the  S-2  sample;  however,  this 
determination  is  open  to  interpretation  due  to  its  subjective  nature.  As  stated  above,  the  lower 
frequencies  are  more  diffused,  whereas  the  higher  ones  are  sharper. 


36 


lil-h  **  15  II  MFM  -ILL  »  I '  I  ?  Rri*AI4.  n  I  ?  it  H?«  U-L  «  3«*T-*fc-2ei 1  ?  it  »3?«  U4,  *  !^-**-;*  1 7  J7  »rr;wi  ALL  *6  I*  *t  DTJM  ALL  n  I*  «l  MT4fc*  M  l  •*  3MT-*H*  1 5  *f  flTrt  M  l-  *  !^-**-3*  H  N*p  w  ALL  #B  IW-W-S*  I*  I  <  *r 


li|  1  L -  ft  ■’ 

0=.T1 

« 

tmr*m  M-  Smii.ikim- 

% 

— r  J  -  4  v 

u 

■!  .  1  "I  -  T  '  teBflllll 

II 

,JL,  t— »  iVl  i  — p  J  -  t  -  V  ^ik.^ 

1  i  J  •  J  • 

.  |«J. . ; 

1  i-wmy.:,  m  ~  ftfe— ill K$  ■  <? $1  1  +  P«H  -  T  -  (?  5. 

•  t 

1 _ _ i 

o 

/=  ioo : 

Hi 

150 

200 

250 

300  350  400 

450  500 

0=  Tla 

ft 

|  # 

• 

i  •  ft 

*  ft  1 

k  k 

O 

II 

C\ 

o 

o 

O 

TP 

42.5° 

56° 

59.5°  63.5°  67.5° 

--4 

o 

o 

O  T10 

| 

m 

L  3  t  k 

i  r  — 

• 

o 

O  T08 

■ 

T 

■Ll  k 

El  f 

•  >  1 

CM 

.  O  T06 

i 

f 

j  ^  L _ 

r 

ft 

l 

o 

O  T04 

m  T  f 

r  _  ^ 

• 

mm  4 

O 

TiTi  »  >*  j's  i'i  "i'i  ti  Ij  w  jj  iV  Ti  ii'u'rt  iTiT  ‘i.<  ti’S'w  »  *7  'n  i’»  ii  **'  n“j7  mT*  jj  «  i>  Jv+  "s.j  i*"-i  >,*  jTi"/  TnV‘ri'i*Tj  w'  Yi7i  11  j*  n  i'i  ji'Ti  31  jf*  j'j  i 


Figure  24  C-Scans  of  S-2  Test  Standard  at  Various  Frequencies  Showing 
Titanium-side  Inside-row 


Figure  25  shows  c-scans  with  frequency  that  best  highlight  the  titanium-side  outside-row. 
Comparison  with  Figure  24  shows  that  the  optimum  frequencies  were  determined  to  be  close 
to,  but  not  generally  equal  to  their  counterparts  for  the  inside-row.  The  proximity  to  the 
bottom  layer  edge  meant  that  the  operator  had  to  optimize  on  the  secondary  peak  on  the  other 
side  from  the  cracks.  It  can  also  be  noted  that  in  this  case  the  discernment  between  cracks  and 
no  cracks  was  based  more  on  magnitude  of  the  lobes  than  in  shape  and  size. 
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Figure  25  C-Scans  of  S-2  Test  Standard  at  Various  Frequencies  Showing  Titanium-side 

Outside-row 


Figure  26  shows  c-scans  with  frequency  that  best  highlight  the  steel-side  inside-row. 
Once  again  it  is  apparent  that  interpretation  of  data  about  steel  fasteners  is  more  complex. 
There  is  a  difference  in  halo  shape  of  the  lobes  about  the  cracks  as  compared  with  those 
without  cracks;  however,  this  difference  is  subtle  in  some  cases.  Interestingly  the  intensity  of 
S04  is  greater  than  that  of  S08  despite  having  a  smaller  notch  crack.  It  is  both  the 
completeness  and  intensity  of  these  halos  that  is  used  to  determine  the  presence  of  the  cracks. 
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Figure  26  C-Scans  of  S-2  Test  Standard  at  Various  Frequencies  Showing  Steel-side 

Inside-row 

Figure  27  shows  c-scans  with  frequency  that  best  highlight  the  steel-side  outside-row. 

Here  again  it  is  a  combination  of  the  intensity,  shape  and  completeness  of  the  lobe  that 
indicates  the  presence  of  a  crack.  The  smallest  crack  fastener  S02  is,  as  expected,  the  most 
difficult  to  discern. 
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Figure  27  C-Scans  of  S-2  Test  Standard  at  Various  Frequencies  Showing  Steel-side 

Outside-row 


The  previous  collection  of  plots  over  frequency  for  the  different  fastener  types  and  rows 
further  from  and  closer  to  the  bottom  edge  of  the  S-2  test  standard  showcases  the  complexity 
of  interpretation  of  MR  Sensor  data.  This  interpretation  has  been  done  in  a  mostly  subjective 
manner,  where  the  operator  gains  experience  in  determining  what  kind  of  signature  to 
anticipate  for  given  scenarios. 

One  of  the  tasks  of  this  program  was  to  build  more  automated  algorithms  and  procedures 
for  quantifying  the  results  and  ultimately  producing  SNRs  as  measureable  quantities  to 
describe  a  particular  feature.  The  results  of  this  task  are  described  in  the  next  section. 

4.1.3  Automated  Defect  Detection 

This  section  describes  discriminatory  features  extracted  from  the  MR  measurements  that 
are  then  used  for  automatic  interpretation  of  individual,  segmented  fastener  image  data.  The 
detection  of  crack  is  primarily  described  using  a  skewness  function  that  best  quantifies  the 
asymmetry  in  the  left  and  right  lobes  of  the  fastener  image.  Classification  can  then  be 
performed  by  thresholding  the  one-  or  multi-dimensional  feature  values.  The  overall 
performance  is  then  quantified  using  SNR  values  calculated  in  the  feature  space. 
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4. 1.3.1  Amplitude -Based  Features 


An  initial  skewness  function,  calculated  using  the  peak  values  of  the  fastener  image,  is 
defined  in  Equation  (33). 


A 

B , 


r  i  ^2 1 , 


(33) 


where  Bt  and  B2  are  peak  values  of  the  two  lobes  of  magnitude  image  calculated  from 
the  in-phase  and  quadrature  MR  measurements  as  depicted  in  Figure  28.  The  values  of  the  Si 
function  calculated  for  each  fastener  image  are  tabulated  in  Table  6.  The  rows  highlighted  in 
red  correspond  to  fasteners  with  crack  and  consequently  have  higher  values  of  Si  whereas  the 
non-highlighted  rows  are  defect- free  fasteners  which  are  characterized  by  a  lower  value  of  Si. 
This  feature  results  in  a  single  false  call  at  fastener  3  which  is  in  truth,  defect- free  but  the 
asymmetry  in  the  data  results  in  a  vale  of  Si  that  is  above  the  threshold.  The  performance 
using  these  features  was  not  optimal  since  the  skewness  function  Si  only  uses  information  in 
the  peak  values  and  ignores  any  information  in  the  overall  shape  of  fastener  images.  The  next 
section  defines  shape  based  features. 


Figure  28  Peak  Amplitude  Definitions  for  Skewness  Function  1 
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Table  6  Classification  Results  for  S-2  Sample  at  400  Hz  and  DA  =  220 


Rivet  # 

Crack 

Length 

Dn 01 

0" 

Dn_02 

0.20” 

I)n_03 

O'1 

I)n_04 

0  22” 

Dn 05 

O'’ 

Dn_06 

0.25” 

1— 1 

■  1 

O'’ 

I)n_08 

0.30" 

Dn 09 

O'’ 

Dn_10 

0.35” 

Dnll 

0” 

4. 1.3.2  Shape-based  features 

Typical  shape  differences  between  the  normal  (defect-free)  fastener  image  and  the 
abnormal  (with  defect)  fastener  image  are  illustrated  in  Figure  29.  Three-dimensional  shape 
information  can  be  represented  in  the  form  of  contour  plots  generated  using  a  set  of 
appropriately  chosen  threshold  values  on  the  image  data.  The  contours  are  represented  by  a 
set  of  vectors  S3  =  [iS)Ai ,  Sth2,...SthN  ] . 


(a)  (b) 

Figure  29  Different  Shape  Information  for  Crack-free  S-2  Sample  Fastener  Images 
Note:  (a)  S2_Ti_dn_09  (no  crack)  and  (b)  S2_Ti_dn_10  (0.35  inch  crack) 

The  contour  plots  of  selected  Titanium  fasteners  in  the  S-2  sample  data  at  400  Hz  are 
shown  in  Figure  30  (a)  -  (f).  These  plots  clearly  indicate  the  need  for  features  based  on  the 
shape  of  local  magnetic  fields  at  the  fastener  site. 
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Figure  30  Contour  Pots  for  S-2  Smple  Fastener  Images  at  400  Hz  and  DA  =  82 

Note:  (a)  crack-free  (b)  0.20  inch  subsurface  (c)  0.22  inch  subsurface  (d)  0.25  inch 
subsurface  (e)  0.30  inch  subsurface  and  (f)  0.25  inch  surface 

4.1. 3.2.1  IMAGE  PREPROCESSING  -  MEAN  SUBTRACTION 

The  rationale  for  the  preprocessing  algorithm,  described  in  this  section,  is  that  a  crack 
free  fastener  image  should  be  symmetric.  The  local  region  of  the  data  around  each  fastener 
was  first  extracted  using  an  automated  image  segmentation  algorithm  and  preprocessed  so 
that  in  the  case  of  a  crack  free  fastener,  the  resulting  image  and  its  contours  are  symmetric. 
The  simple  processing  steps  to  accomplish  shape-based  feature  extraction  for  S-2,  S-4 
samples  are  outlined  in  the  following  steps. 
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1 .  Derive  fastener  images  from  raw  MR  signals  with  ODA=0  (Re)  and  ODA=90(Im) 

2.  Extract  1-D  signal  through  peak  of  fastener  image  (Note  that  is  not  symmetric  about  0  - 
see  Figure  31) 

3.  Generate  corresponding  zero-mean  signal  by  subtracting  the  mean  from  the  signal. 

4.  Generate  symmetric  Magnitude  signal  image  using M  =  (Re2  +  Im2  )12 


The  one-dimensional  line  scans  before  and  after  (magnitude,  real,  and  imaginary 
components)  mean  subtraction  are  shown  in  Figure  32  for  a  defect-free  fastener. 


Figure  3 1  Illustration  of  1-D  Signal  Extraction. 


NP-Rivel  Outside  Row  #5  AP-RIvet  Outside  Row  #3 


(a) 


(b) 


Figure  32  One-Dimensional  Signals  of  Magnitude(blue),  real(green),  and  Imaginary(red) 

Components 

Note:  (a)  before  preprocessing  and  (b)  after  preprocessing 


The  fastener  images  of  fasteners  with  different  cracks,  before  and  after  preprocessing,  are 
shown  in  Figure  33  (a)  and  (b). 
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Figure  33  Images  of  Fasteners  with  Different  Crack  Dimensions-  (from  top  to  bottom) 
0.00  inch,  0.20  inch,  0.22  inch,  0.25  inch  and  0.30  inch 

Note:  (a)  before  and  (b)  after  preprocessing 

Five  shape-based  features  are  defined  and  expressed  in  Equations  (34)  to  (38). 
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(34) 


6=n 


Ok 

wLi 


i= 1 

F  Area  (X,- X,)-(Y,-Y2) 
3  | Slope]  |r3-y4|/pr3-A'4| 

^4  ~  Hmhw  (37) 


F5=HTail=YA-Y2 


(35) 

(36) 


(38) 


The  parameters  used  in  these  features  are  illustrated  graphically  in  Figure  34(a)  to  (d). 
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The  features  calculated  using  400  Hz  data  from  S-2  samples  were  used  to  devise  an 
appropriate  classification  rule  for  crack  detection.  Not  all  features  were  found  to  be  useful. 
The  optimal  features,  selected  manually,  were  different  in  the  case  of  signals  from  inside  and 
outside  fasteners.  Optimal  features  for  outside  fastener  data  were  Feature  2  and  Feature  3 
and,  for  inside  fastener  data,  were  Feature  2  and  Feature  4.  The  selected  features  and 
corresponding  one-dimensional  signals,  for  inside  and  outside  fastener  signals,  are  plotted  in  a 
scatter  plot  in  Figure  35  and  Figure  36,  respectively.  The  green  dots  in  the  plots  are  the  no¬ 
defect  fasteners.  The  red  flags  are  the  defective  fasteners.  Similar  analysis  was  performed  on 
data  from  S-4  sample  obtained  at  100  Hz  frequency. 


2-Dim  Feature  Space  Classification  Results  for  Inside  S-2 


Inside  Row  rivets 


Mag  Line  Plots  for  Inside  S-2 


2-Dim  Feature  Space  Classification  Results  for  Outside  S-2 


Outside  Row  rivets 


Mag  Line  Plots  for  Outside  S-2 


Figure  35  Two  Dimensional  Classification  Results  for  S-2  Sample  (0.16  inch  Top  Layer 
thickness)  and  Corresponding  1-D  Signal  Plots 
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2-Dim  Feature  Space  Classification  Results  for  Inside  S-4 


2-Dim  Feature  Space  Classification  Results  for  Outside  S-4 


Mag  Line  Plots  for  Outside  S-4 


Figure  36  Two  Dimensional  Classification  Results  for  S-4  Sample  (0.25  inch  Top  Layer 
thickness)  and  Corresponding  1-D  Signal  Plots 

4.1. 3.2.2  IMAGE  PREPROCESSING  -  OPTIMUM  DETECTION  ANGLE  (ODA) 

In  this  approach  the  in-phase  and  quadrature  MR  sensor  signals  are  mixed  to  generate  a 
real  valued  signal  using  a  detection  angle  a  according  to  equation  (39), 


S  ( a )  =  S0  cos(a)  +  Sn/  sin(a) , 

/2 


(39) 


where  Sa, 


So  and  SA 


are  the  mixed,  in-phase  and  quadrature  signals,  respectively. 


A  method  to  automate  the  estimation  of  optimum  value  of  detection  angle  a,  was 
developed  to  suppress  the  dominant  fastener  image  and  highlight  the  data  corresponding  to 
the  crack.  The  optimum  detection  angle  was  obtained  by  minimizing  the  peak-to-peak 
amplitude  of  the  mixed  signal.  The  detection  angle  preprocessing  algorithm  consists  of  the 
following  steps: 

1)  Initialize  detection  angle  to  0. 

2)  Increment  detection  angle  (0  to  180  with  step  size  of  10  degrees). 

3)  For  each  detection  angle  - 

1)  Generate  the  mixed  signal  using  equation  (4. 1 .3-7). 

2)  Segment  Fastener  image  in  the  mixed  signal  automatically. 

3)  Calculate  S  k  (a)  =  peak-to-peak  value  between  the  two  lobes. 
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4)  Plot  S  k  (a)  versus  a  -  the  angle  corresponding  to  minimum  of  this  curve  is 
treated  as  the  optimum  detection  angle  (ODA) 

4)  Generate  mixed  image  using  ODA  determined  in  step  3. 


The  plots  generated  in  step  3)  to  determine  ODA  and  corresponding  mixed  images  after 
applying  ODA  are  shown  in  Figure  37  for  a  set  of  S-2  Titanium  OUTSIDE  fasteners  at  400 
Hz  frequency.  In  these  images,  fasteners  2,  4,  6,  and  8  correspond  to  0.20  inch,  0.22  inch,  0.25 
inch,  and  0.30  inch  subsurface  defects,  respectively.  Also,  it  should  be  noted  that  the  ODA  is 
different  for  data  at  different  frequencies.  From  the  minima  of  the  Sph  (a)  curves  it  is  seen 

that  at  400  Hz  the  ODA  is  70  degrees,  while  at  100  Hz,  the  ODA  is  30  degrees. 
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Figure  37  Estimation  of  Optimum  Detection  Angle  Algorithm  and  Resulting  Images  for 
S-2  Titanium  Outside  Fasteners  at  400  Hz 

Although  the  ODA  processed  images  do  not  appear  to  be  informative,  the  asymmetry  in 
the  lobes  is  enhanced  which  is  more  visible  in  the  line  plot  along  the  fastener  center  as 
illustrated  in  Figure  38. 


Figure  38  Extraction  of  Line  Plot  from  ODA  Processed  Images 

The  extracted  line  plots  from  different  fastener  images  are  shown  in  Figure  39.  The 
enhanced  asymmetry  in  the  lobes  is  reflected  in  these  line  plots.  More  importantly  it  is  seen 
that  the  curves  from  defective  fasteners  (red  lines)  are  clearly  separated  from  line  plots 
obtained  from  defect-free  fasteners  (blue  lines).  The  corresponding  line  plots  for  outside  and 
inside  Steel  fastener  data  from  S-2  sample  are  shown  in  Figure  40  at  a  frequency  of  400  Hz 
and  ODA  of  90  degrees. 
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S-2  (Ti-Outside-400Hz)Signal-Out-Row:ODA  selected70 


S-2  (Ti-ln  side-400  Hz  )Signal-ln-Row:ODA  selected70 


(b) 


Figure  39  Line  Plots  for  S-2  Ti  Fasteners  at  400  Hz 
Note:  (a)  outside  and  (b)  inside  after  applying  ODA  =  70  where  blue  lines  are  defect  free 
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S-2  Freq-400  Outside  Stl:ODA  selected 90 


Figure  40  Line  plots  for  S-2  Steel  Fasteners  at  400  Hz 
Note:  (a)  outside  and  (b)  inside  after  applying  ODA  =  90  where  blue  lines  are  defect  free 

Results  of  processing  data  from  S-2  sample  with  Ti  and  Steel  fasteners  at  other 
frequencies  (200  Hz,  300  Hz,  and  400  Hz)  are  shown  in  Figure  41  and  Figure  42, 
respectively.  The  line  plots  for  defective  fasteners  (red  lines)  are  clearly  separated  from  the 
line  plots  of  defect-free  fasteners  (blue  lines). 
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Figure  41  Line  plots  for  S-2  Ti  Fasteners  at  Frequencies  from  200  Hz  to  400  Hz  where 

Blue  Lines  are  Defect  Free 


54 


S-2  Freq-200  Outside  Stl:ODA  selected80 


S-2  Freq-200  Inside  Stl:ODA  selected 80 


04  06  08 

X  Distanc«(inch) 


1.2  14  16 

X  OstanceOnch) 


0  15 


S-2  Freq-300  Outside  Stl:ODA  selected 270 


S-2  Freq-300  Inside  Stl:ODA  selected270 


0.2  04  06  08  1 

X  Distance<inch) 


1.2  14  1.6 

X  Distance<incti) 


S-2  Freq-400  Outside  Stl:ODA  selected 90 


S-2  Freq-400  Inside  Stl:ODA  selected90 


02  04  06  08  1 

X  Distarvc«(incfi) 


1  1.2  14  1.6 

X  DistanceOncft) 


Figure  42  Line  Plots  for  S-2  Steel  Fasteners  at  Frequencies  from  200  Hz  to  400  Hz  where 

Blue  Lines  are  Defect  Free 


New  features  based  on  these  line  plots  were  introduced  and  investigated.  The  new 
features  are  defined  in  Equations  (40)  and  (41). 


1  TV  1  M 

FMAE  =  —  X  LD  PI  “  —  X  LND_j  PI 


(40) 


1  TV  1  iVl 

FmSE  (I'D  P  ]  “  ~  X  LND  _  j  PI) 


(41) 
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Fmae  is  the  feature  calculated  from  mean  absolute  error  and  the  FMse  is  calculated  from 
mean  square  error.  In  Equations  4. 1.3-8  and  4. 1.3-9,  N  is  the  length  of  the  mixed  1-D  data 
vector,  M  is  the  total  number  of  nondefective  fasteners.  LD[i]  is  the  /th  pixel  of  mixed  data 
from  defective  fastener  and  Lndj[i]  is  the  /th  pixel  of  the  /th  mixed  data  from  a  non-defective 
fastener.  A  classification  rule  was  devised  based  on  partitioning  the  two-dimensional  feature 
vector  in  the  feature  space.  The  classification  results  using  data  at  different  frequencies  for  Ti 
and  Steel  fasteners,  are  shown  in  Figure  43  and  Figure  44,  respectively. 


Feature-MAE  Feature  MAE 


Figure  43  Classification  Results  for  S-2  Ti  Fasteners  at  Frequencies  from  200  Hz  to  400 
Hz  where  Green  Dots  are  Defect-Free  Fasteners. 
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Figure  44  Classification  Results  for  S-2  Steel  Fasteners  at  Frequencies  from  200  Hz  to 
400  Hz  where  Green  Dots  are  Defect  Free 


4, 1,3.3  Classification  Results 


The  results  presented  in  this  section  show  that  the  automated  signal  analysis  algorithm 
developed  have  considerable  promise  for  classifying  MR  data  from  Ti  and  steel  fasteners.  In 
the  case  of  steel  fasteners  the  data  are  classified  with  100  percent  POD  using  the  200  Hz  and 
400  Hz  data.  The  classification  performance  can  be  assessed  quantitatively  using  SNR 
values. 

It  is  seen  in  the  scatter  plots  in  feature  space,  the  points  corresponding  to  the  good 
fasteners  are  well  separated  from  those  due  to  defective  fasteners.  The  classification 
performance  can  be  quantified  using  signal  to  noise  ratio  (SNR),  as  defined  in  Equation  (42): 


57 


(42) 


snrm 


/  2 
/<Toi 


where  M  is  the  dimensionality  of  the  feature  vector.  Here  we  use  global  features  FMae  and 
Fmse,  so  M  =  2  ,  m0j  and  <r0;2  are  the  mean  and  variance  of  the  defect-free  fastener  data  (the 

green  circle  cluster  shown  in  Figure  45).  The  parameters  used  in  SNR  definition  are 
illustrated  in  Figure  45. 


rr#2 in  Fj  direction  Feature  2 

erwin  Fj  direction 


Figure  45  Illustration  of  the  SNR  Definition 


The  SNR  values  are  computed  using  scatter  plots  obtained  using  amplitude-based 
(magnitude  -  MAG)  features  defined  on  image  preprocessed  data  as  well  as  shape-based 
features  defined  on  the  ODA  processed  data.  These  values  are  compared  with  typical  SNR 
values  obtained  using  eddy  current  data  in  Table  7. 
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Table  7  SNR  Results  Comparison  for  ECT,  MAG,  and  ODA-Based  Methods 


S-2  Test  Standard 

(Ti,  400  Hz) 

S-4  Test  Standard 

(Ti,  100  Hz) 

SNR 

SNR*=jz,(F‘-" '")2 

V  «=1  ^0  i 

Thickness 

Row 

Methods 

0.30  inch 

0.25  inch 

0.22  inch 

0.20  inch 

0.16  inch 

Outside 

ECT 

3.00 

1.50 

1.50 

1.50 

MAG 

45.0 

18.0 

14.0 

9.60 

ODA 

317 

49.8 

16.4 

17.1 

0.16  inch 

Inside 

ECT 

2.50 

1.50 

1.00 

1.00 

MAG 

20.6 

8.20 

4.40 

4.00 

ODA 

292 

64.5 

17.1 

18.2 

0.25  inch 

Outside 

ECT 

2.00 

1.00 

1.00 

1.00 

MAG 

30.1 

11.4 

7.70 

5.50 

ODA 

114 

13.3 

4.60 

11.4 

0.25  inch 

Inside 

ECT 

2.00 

1.50 

1.00 

1.00 

MAG 

7.20 

4.60 

2.50 

3.30 

ODA 

44.1 

9.70 

4.10 

5.10 

The  SNR  values  are  directly  proportional  to  the  POD  and  results  presented  demonstrate 
that  the  POD  is  increased  significantly  by  processing  the  MR  signals  using  ODA.  The  table 
also  shows  the  advantage  of  MR  sensors  over  ECT  in  subsurface  crack  detection. 

The  SNR  versus  crack  area  plots  for  all  three  methods  are  shown  in  Figure  46  for  the  S-2 
sample  inside  Ti  fastener  data.  The  crack  (notch)  areas  are  derived  from  the  trapezoidal  and 
triangular  cross-section  shapes  as  shown  in  Figure  15.  The  figure  clearly  shows  the  advantage 
of  automated  processing  of  image  data  over  the  manually  interpreted  eddy  current  data. 
Further  it  is  seen  that  preprocessing  the  MR  data  using  ODA  is  particularly  advantageous 
since  this  method  enhances  the  defect  indication  relative  to  the  signal  generated  by  the 
fastener. 
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Figure  46  SNR  Improvement  after  Applying  ODA-based  Method.  SNR  is  Plotted  Versus 
Crack  Cross-sectional  Area.  The  Points  are  also  Labeled  with  Top  Crack  Dimension. 

Implementing  the  steps  1)  ODA  estimation,  2)  Mixing  with  ODA,  3)  Feature  extraction, 
4)  generation  of  scatter  plots,  and  5)  SNR  calculation,  on  the  MR  data  obtained  at  different 
frequencies,  one  can  generate  a  plot  of  SNR  versus  frequency  to  determine  the  optimum 
frequency  that  maximizes  the  POD.  The  SNR  versus  frequency  curves  for  inside  and  outside 
fastener  data  using  magnitude-based  features  are  shown  in  Figure  47  and  corresponding  plots 
obtained  using  ODA-based  results  are  shown  in  Figure  48. 
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(a) 


(b) 

Figure  47  SNR  Versus  Frequency  for  S-2  Titanium  Sample  using  MAG-based  Method 
Note:  (a)  Inside  row  fasteners  and  (b)  outside  row  fasteners 
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SNR  vs.  Frequency(S-2  Sample  Inside  Row  Rivets) 
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Figure  48  SNR  Versus  Frequency  for  S-2  Titanium  Sample  using  ODA-based  Method 
Note:  (a)  INSIDE  fasteners  (b)  OUTSIDE  fasteners 


The  optimum  frequencies  occur  at  450  Hz  and  500  Hz  for  inside  and  outside  S-2 
Titanium  fasteners.  Figure  49  shows  the  mixed  data  line  plots  and  scatter  plots  in  two- 
dimensional  feature  space  for  the  S-2  Ti  fastener  data  at  these  optimum  frequencies. 
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Figure  49  Line  Plots  and  Scatter  Plots  in  Two-Dimensional  Feature  Space  for  the  Peak 
Occurring  Position  for  S-2  Sample  Titanium  Fasteners. 
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4. 1.3.4  Extension  to  multiple  line  plots 


Results  presented  so  far  are  based  on  features  calculated  from  single  line  scans  selected 
carefully  through  the  center  of  the  rivet.  In  a  practical  situation,  this  would  correspond  to  a 
specific  sensor  in  the  GMR  array  sensor  to  study  the  sensitivity  of  the  algorithm  to  location  of 
the  line  scan  with  respect  to  the  fastener,  the  performance  was  tried  on  multiple  data  scans 
selected  as  illustrated  in  Figure  50.  The  central  line  L3  is  the  single  central  line  scan  used  in 
the  previous  section.  In  this  study,  five  consecutive  scan  lines  centered  about  the  L3  are  used 
in  feature  extraction  and  classification. 


ODA-O-Row  #9  ODA=70 


2.8  3  3.2  3.4  3.6  3.8 

XDistsnce(lnch) 


+2 

+1 

0 

-1 
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Figure  50  Multiple  Line  Plots  across  a  MR  Sensor  Fastener  Image 


The  data  acquired  on  the  S-2  standard  sample  with  Titanium  fasteners  at  400  Hz 
frequency  was  analyzed  using  multiple  line  scans  and  results  are  presented  below.  The  line 
plots  along  with  scatter  plots  of  the  features  are  shown  in  Figure  5 1  and  Figure  52  for  inside 
and  outside  fasteners,  respectively. 
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S-2  Freq-400  Inside  Ti:ODA-70  LinePos:  5  total 


X  Distance(inch) 


Figure  5 1  Collection  of  All  Line  Scans  for  Defect- free  and  Defective  Inside  Fasteners  and 

Scatter  Plot  of  Features 


S-2  Freq-400  Outside  Ti:ODA-70  LinePos:5  total 


X  Distance(inch) 


Sx30Freq-400  Outside  Ti:ODA-70  LinePos:  5  total 


Figure  52  Collection  of  All  Line  Scans  for  Defect- free  and  Defective  Outside  Fasteners 

and  Scatter  Plot  of  Features. 

The  SNR  values  for  the  cases  of  single  and  multiple  line  scans  are  shown  in  Table  8.  The 
SNR  values  obtained  earlier  for  single  line  scan  are  also  included  for  comparison. 
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Table  8  SNR  Values  for  Different  Methods 


S-2  (Ti,  0.16” 

,  400Hz) 

SNR(X:1) 

Skin 

Thickness 

Fastener 

Row 

0.30  inch 

0.25  inch 

0.22  inch 

0.20  inch 

0.16  inch 

Outside 

3(eddy  current) 

45.0(Mag) 

316.7(Single) 

891.4(Multi) 

1.5(eddy  current) 

18.0(Mag) 

49.8(Single) 

144.8(Multi) 

1.5(eddy  current) 

14.0(Mag) 

16.4(Single) 

49.3(Multi) 

1.5(eddy  current) 

9.6(Mag) 

17.1(Single) 

58.1  (Multi) 

0.16  inch 

Inside 

2.5(eddy 

current) 

20.6(Mag) 

292.1(Single) 

204.3(Multi) 

l.5(eddy  current) 

8.2(Mag) 

64.5(Single) 

45.1  (Multi) 

l(eddy  current) 

4.4(Mag) 

17.1  (Single) 

12.8(Multi) 

l(eddy  current) 

4.0(Mag) 

18.2(Single) 

13.3(Multi) 

In  this  study,  3,5,7  and  1 1  consecutive  scan  lines  about  the  central  line  L3  are  used  in 
feature  extraction  and  classification. 

The  data  acquired  on  the  S-2  standard  sample  with  Titanium  fasteners  at  400  Hz 
frequency  was  analyzed  and  results  are  presented  below.  Figure  53  and  Figure  54  present  the 
results  in  the  form  of  line  plots  and  scatter  plots,  respectively.  The  multiple  line  positions 
were  marked  on  the  scatter  plots,  so  the  SNR  versus  the  line  positions  can  be  read  directly 
from  the  plots.  From  Figure  54,  it  is  seen  that  the  central  line  always  gives  a  better  SNR  value 
for  outside  fasteners  while  for  inside  row  fasteners,  the  line  plots  around  the  central  one  give  a 
better  SNR. 
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3-2  Freq-400  Outside  Ti:ODA-70  LNum:1 


S-2  Freq-400  Inside  Ti:ODA-70  LNum:1 
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3-2  Freq-400  Outside  Ti:ODA-7Q  LNum:5 


S-2  Freq-400  Inside  Ti:ODA-70  LNum:5 


32  3.4  3.5 

X  CSSt0rtC*[irwtTl 


26  28  3  32 

x  Distance  inert  > 


(c) 


S-2  Freq-400  Outside  Ti:ODA-70  LNum;7 


S-2  Freq-400  Inside  Ti:ODA-70  LNum;7 
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(d) 


66 


(e) 

Figure  53  Collection  of  All  Line  Scans  for  Defect-free  and  Defective  Fasteners 

Note:  Outside  Row  (Left)  and  Inside  Row  (Right):  (a)  single  line  (central),  (b)  3  lines,  (c)  5 
lines,  (d)  7  lines  and  (e)  1 1  lines 
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Feature-MSE  Feature^MSE  Feature-MSE  Feature-MSE 


(a) 


(b) 


(c) 


(d) 
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(e) 

Figure  54  Scatter  Plot  of  Features  for  Multiple  Lines 

Note:  Outside  Row  (Left)  and  Inside  Row  (Right):  (a)  single  line  (central),  (b)  3  lines,  (c)  5 
lines,  (d)  7  lines  and  (e)  1 1  lines 


The  SNR  values  for  single  and  multiple  line  plots  are  shown  in  Table  9.  The  SNR  values 
obtained  using  amplitude  based  features  are  also  included  for  comparison. 

Table  9  SNR  Values  for  Different  Methods 


S-2  (Ti,  0. 
400 

16  inch, 
Hz) 

Signal-to-Noise  (X:1) 

Skin 

Thickness 

Fastener 

Row 

0.30” 

0.25” 

0.22” 

0.20” 

0.16” 

Outside 

3  (eddy  current) 

1.5  (eddy 
current) 

1.5  (eddy 
current) 

1.5  (eddy 
current) 

45.0(Mag) 

18.0(Mag) 

14.0(Mag) 

9.6(Mag) 

316.7(Single) 

49.8(Single) 

16.4(Single) 

17.1  (Single) 

455.9(3  lines) 

74.3(3  lines) 

23.6(3  lines) 

27.3(3  lines) 

355.5(5  lines) 

56.5(5  lines) 

18.7(5  lines) 

22.3(5  lines) 

290.9(7  lines) 

46.4(7  lines) 

15.8(7  lines) 

18.7(7  lines) 

78.2  (11  lines) 

13.1  (11  lines) 

4.8  (11  lines) 

5.6  (11  lines) 

0.16” 

Inside 

2. 5  (eddy 
current) 

1.5  (eddy 
current) 

1  (eddy 
current) 

1  (eddy 
current) 

20.6(Mag) 

8.2(Mag) 

4.4(Mag) 

4.0(Mag) 

292.1  (Single) 

64.5(Single) 

17.1  (Single) 

18.2(Single) 

252.5(3  lines) 

54.5(3  lines) 

15.2(3  lines) 

15.8(3  lines) 

189.9(5  lines) 

41.7(5  lines) 

11.8(5  lines) 

12.2(5  lines) 

112.9(7  lines) 

25.5(7  lines) 

7.5(7  lines) 

7.8(7  lines) 

43.8  (11  lines) 

10.8(11  lines) 

3.4  (11  lines) 

3.5  (11  lines) 
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4.2  OPTIMIZATION  OF  PARAMETERS  THROUGH  MODELING 


4.2.1  Model  validation 

The  finite  element  model  developed  for  MR  inspection  of  multilayer  test  samples  was 
first  validated  using  experimental  measurements  from  both  the  calibration  sample  with  notch 
and  a  S-2  standard  sample.  The  validation  results  are  presented  in  this  section. 

4.2.1. 1  Model  validation  using  notch  calibration  standard  sample. 

The  calibration  sample  is  made  of  standard  aluminum.  The  sample  is  a  plate  of 
dimensions  16  inches  X  10  inches  and  of  thickness  0.625  inch.  A  notch  of  dimensions  0.05 
inch  width  and  0.25  inch  depth  is  machined  at  the  center  of  the  sample.  A  photograph  of  the 
sample  is  shown  in  Figure  55. 


(a)  (b) 

Figure  55  Photographs  of  (a)  the  Sensor  Head  at  the  Zero-Balancing  Calibration  Position 
and  (b)  an  Overhead  View  of  the  Calibration  Sample  with  Slot  Down  the  Middle.  The  Sensors 
and  Excitation  Current  are  Aligned  Parallel  to  the  Slot. 

The  three-dimensional  domain  of  the  FE  model  and  the  corresponding  mesh  are  shown  in 
Figure  56  and  Figure  57,  respectively. 
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Figure  56  Sample  Geometry  Used  in  the  Model 


Figure  57  (a)  SidView  and  (b)Top  View  of  the  Finite  Element  Mesh 

The  FEM  solutions  produce  the  electric  scalar  potential  V  and  magnetic  vector  potential 
A  at  each  node  in  the  solution  region.  The  magnetic  flux  density  B=  VXA  is  computed  and 
its  normal  component,  is  compared  with  the  measured  MR  signal.  Figure  58  shows  three- 
dimensional  plot  of  the  real  and  imaginary  part  of  the  model  computed  magnetic  flux  density, 
The  two-dimensional  scan  across  the  slotted  sample  is  applied  to  obtain  the  sensor  signal. 
Figure  59  shows  geometry  of  the  scan. 
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Figure  58  Modeling  Results  of  the  Normal  Component  of  the  Magnetic  Flux  Density  (a) 

Real  Part  (b)  Imaginary  Part 


Excitation 

Current 


Figure  59  Geometry  of  Scan  Plan 


Experimental  signals  for  the  sample  were  available  for  various  demodulation  angles  as 
described  in  section  4. 1.1. 3.  The  in-phase  and  quadrature  components  were  demodulated  to 
generate  a  mixed  signal  using  varying  demodulation  angle  5.  The  mathematical  form  for  the 
calibrated  demodulated  signal  for  sensor  n  at  arbitrary  detection  angle  8  is: 

{So  ■  S„  )s  =  a ,  ((S„  ■  S„  )0  -  ( S •  s;)o)cos(5)  +  a,  ((S,  ■  S„)I/2  -  {s„  ■  S„°)j sin(<5)  (43) 


where  (S„  -S,)s  =^^cos(<fM -S) . 

In  finite  element  modeling,  the  zero  balancing  term  (sq  •  Sn°^  and  the  scaling  factor  an 
are  0  and  1,  respectively  since  these  two  terms  are  only  useful  in  experimental  data  analysis. 
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Since  the  initial  phase  of  excitation  current  in  the  FE  models  is  different  from  what  is 

"  "  71 

used  in  experiments,  the  detection  phase  5  is  calibrated  to  5  ,  where  8  =  —  -  8 .  A 

comparison  of  model  and  experimental  signals  are  shown  in  Figure  60,  where  the  calibrated 
detection  phase  8  spans  half  a  cycle  from  0  to  180  degrees  in  steps  of  30  degrees.  We  can 
conclude  that  the  model  result  accurately  predicts  the  experimental  measurements. 
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(a)  model  prediction  (b)  experimental  signal 

Figure  60  Comparison  of  Demodulated  and  Calibrated  Signal  Across  the  Calibration  Slot 


4.2. 1.2  Model  validation  using  S-2  Standard  Sample 


A  schematic  of  the  geometry  of  S-2  standard  sample  with  H  section  crack,  sensor 
configuration  and  scan  direction  are  shown  in  Figure  61. 
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Figure  61  Geometry  for  S-2  Sample. 

The  fastener  material  in  this  sample  is  Titanium.  The  side  view  shows  top  layer  truncated 
in  accordance  with  the  S-2  sample  geometry.  The  model  also  simulates  the  gap  between  the 
two  top  skins  of  S-2  sample.  The  current  is  parallel  to  the  notch  length  and  the  scan  direction 
is  perpendicular  to  it.  The  H  section  crack  is  not  applied  in  the  FEM  so  that  the  comparison  is 
between  simulated  and  experimental  signals  for  the  crack  free  rivet.  The  optimized  three- 
dimensional  FE  mesh  for  this  sample  with  edge  is  shown  in  Figure  62. 
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(a) 


Figure  62  FE  Mesh  for  S-2  Sample:  (a)  Side  View  (b)  Top  View 

The  modeling  results  at  400  Hz  excitation  frequency  are  shown  in  Figure  63.  The  real 
and  imaginary  of  the  normal  component  of  magnetic  flux  density  are  shown  in  Figure  63(a) 
and  (b).  The  magnitude  of  the  magnetic  flux  density  is  presented  in  Figure  63(c). 
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(a) 


(b) 


(c) 

Figure  63  Magnetic  Flux  Density  for  S-2  Sample  Geometry  in  Figure  62  at  400  Hz  (a) 
Magnitude  (b)  Real  Component  (c)  Imaginary  Component 


The  corresponding  experimental  data  from  the  rivet  is  shown  in  Figure  64. 
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Note:  The  line  scan  is  marked  by  the  red  dashed  line. 

A  comparison  between  simulation  and  experimental  results  is  presented.  Line  scans  in  the 
x-direction  across  the  center  of  the  fastener  is  extracted  from  the  real  and  imaginary  image 
data  (both  the  simulation  and  experiment)  which  are  then  mixed  using  different  demodulation 
or  detection  angles.  The  detection  angle  a  and  the  mixed  MR  sensors  signal  is  given  by 

Sa  =  S()  cos(a)  +  .S’,  sin(a) ,  (44) 

72 

where  Sa  ,  S0  and  Sn/  are  the  mixed,  in-phase  and  quadrature  signals,  respectively. 
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The  mixed  model  and  experiment  signals  are  compared  for  a  number  of  detection  angle 
values  in  Figure  65  through  Figure  67.  The  detection  angle  for  experimental  signals  is 
changed  from  80  to  110  degrees,  120  to  150  degrees  and  160  to  190  degrees  in  the  three 
figures,  respectively.  Both  the  magnitude  and  shape  of  the  signals  are  remarkably  similar. 
These  validation  results  for  the  model  allow  the  use  of  the  model  in  conducting  the  series  of 
simulations  for  the  parametric  study  described  in  the  following  sections. 
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Figure  65  Comparison  between  Experimental  and  Modeling  Signal  with  Various 

Detection  Phases. 
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Figure  66  Comparison  between  Experimental  and  Modeling  Signal  with  Various 

Detection  Phases 
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Figure  67  Comparison  between  Experimental  and  Modeling  Signal  with  Various 

Detection  Phases 


4.2.2  Parametric  Studies 

4.2.2. 1  Introduction 


Once  the  FEMs  were  validated,  then  studies  of  parameters  could  be  pursued  to  assess 
optimum  scanning  conditions  for  MR  sensors.  Table  10  lists  parameters  of  interest  that  could 
be  evaluated  within  the  program  constraints.  The  table  contains  the  range  of  the  parameters 
that  would  be  examined  in  the  modeling  studies. 
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Table  10  Parameter  Study  Table 


Parameter 

Range  of  Value  Tested 

Frequency  (kHz) 

0.1,  0.2,  0.4,  0.5,  0.7,  1.0,  2.0,  3.0, 

5.0,  7.0 

Sensor  liftoff  (inch) 

0.0050,  0.0095,  0.015,  0.1,  0.15, 

0.2 

Conductivity  of  top  layer 
(%IACS) 

Between  28  and  33 

Conductivity  of  bottom  layer 
(%IACS) 

Between  30  and  36 

Conductivity  of  fastener 
(%IACS) 

Titanium  (6A1,  4Y):  1.0,  2.2,  3.1 

Fastener  to  edge  distance  (inchO 

0.40,  0.50,  0.60,  0.70,  0.80 

Crack  dimensions  (inch) 

0.20,  0.22,  0.25,  0.30 

4. 2.2. 2  Frequency  Effect 

The  parameter  that  has  a  significant  impact  on  the  measured  MR  signal  is  the  excitation 
frequency.  Since  most  aircraft  geometry  contains  multiple  layers  with  cracks  in  second  and 
third  layer  the  operating  frequencies  is  largely  at  the  low  end.  The  range  of  frequencies 
considered  in  this  study  are:  0.1,  0.2,  0.4,  0.5,  0.7,  1.0,  2.0,  3.0,  5.0,  and  7.0  kHz.  The 
geometry  of  the  S2  sample  and  sensor  configuration  and  finite  element  mesh  are  shown  in 
Figure  61  and  Figure  62.  A  0.3  inch  second  layer  defect  is  introduced  under  the  fastener 
through  the  second  layer  on  the  side  of  the  stringer  edge.  The  normal  component  of  magnetic 
flux  density  Bz  is  calculated  for  each  frequency,  for  both,  the  crack  free  case  and  0.3-inch 
crack  case. 

Figure  68  shows  the  simulation  results  at  a  frequency  of  0.1  kHz.  The  images  in  3(a), 
3(c)  and  3(e)  are  the  magnitude,  real  and  imaginary  parts  of  the  normal  component  of  flux 
density  for  crack  free  case  at  a  liftoff  of  0.0095  inch.  Images  in  3(b),  3(d)  and  3(f)  are  the 
corresponding  magnitude,  real  and  imaginary  parts  of  the  normal  component  of  flux  density 
for  0.3-inch  crack  geometry.  Figure  69  through  Figure  73  present  similar  model  predicted 
results  for  frequencies  of  400,  500,  700,  2000,  and  7000  Hz.  The  real  imaginary  and  absolute 
magnitude  of  the  normal  component  magnetic  flux  density  for  both  defect  free  and  0.3-inch 
crack  are  presented. 
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(e)  (f) 

Figure  68  Simulation  Results  for  100  Hz  Frequency 
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(e) 


(f) 


Figure  69  Simulation  Results  for  400  Hz  Frequency 
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(e)  (f) 

Figure  70  Simulation  Results  for  0.5  kHz  Frequency 


85 
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Figure  71  Simulation  Results  for  0.7  kHz  Frequency 
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Figure  72  Simulation  Results  for  2.0  kHz  Frequency. 
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Figure  73  Simulation  Results  for  7.0  kHz  Frequency 


4.2.2.2.1  SKEWNESS  CALCULATIONS 


Although  the  defect  indication  is  visible  in  most  of  the  images,  a  simple  method  for 
determining  the  optimum  value  of  frequency  is  presented  in  this  section  via  a  quantitative 
measure  of  the  asymmetry  in  the  two  lobes  of  the  image.  A  simple  skewness  function  for 
quantifying  the  asymmetry  is  calculated  based  on  the  peak  values  of  the  two  lobes  of  the 
fastener  image,  and  is  defined  in  Equation  (45). 
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max{ 

B1 

•> 

b2 
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min{ 

•> 

b2 

} 

(45) 


The  parameters  of  this  function  are  shown  pictorially  in  Figure  74.  The  value  of  Si  is 
calculated  for  the  fastener  image  obtained  at  each  frequency  and  is  plotted  in  Figure  75. 


Figure  74  Surface  Plot  of  the  Image  Data  Showing  the  Asymmetry  in  Two  Lobes  of 

Fastener  Image 


Frequency(kFlz) 


Figure  75  Skewness  Versus.  Frequency  for  Results  with  0.3-inch  crack. 


From  Figure  75  it  is  seen  that  at  lower  frequencies  the  skewness  values  are  much  higher 
than  1  and  the  crack  is  easily  detected.  At  higher  frequency,  the  crack  is  hard  to  see  due  to  the 
skin  depth  effect.  The  skewness  plots  show  that  the  optimal  operating  frequency  is  0.4  kHz. 
Further,  the  modeling  results  indicate  that  as  the  frequency  increases,  the  images  are  sharper 
and  the  edge  effect  vanishes,  but  the  crack  detection  ability  decreases  because  of  the  skin 
depth  effect. 
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4.2.2.22  SNR  CALCULATIONS 


An  alternate  method  for  optimizing  frequency  of  operation  is  to  select  a  frequency  that 
maximizes  the  POD  or  equivalently  the  SNR  of  the  signal.  In  this  study,  SNR  values  are 
calculated  with  simulation  results  obtained  at  frequencies:  100  Hz,  400  Hz,  500  Hz,  700  Hz 
and  1000  Hz.  Both  crack- free  fastener  and  the  case  of  a  fastener  with  0.3  inch  subsurface 
crack  were  simulated.  The  fastener  image  at  each  frequency  is  analyzed  using  the  procedure 
described  in  section  4. 1.3. 2.2.  The  image  data  was  preprocessed  using  the  optimum  detection 
angle  at  each  frequency  and  the  line  scans  across  the  center  of  the  fastener  was  extracted.  The 
same  features  developed  in  section  4. 1.3. 2.2  using  experimental  data  from  S-2  sample  are 
applied  for  simulated  signal.  The  features  are  as  shown  in  Equations  (45)  and  (47),  which  are 
the  same  as  (4. 1.3-8)  and  (4. 1.3-9). 


1  N 

FuAE  =  Nt; 


LD\i\ 


1  M 

— Yl 


ND_j 


P] 


7=1 


(46) 


^  i=\ 


M 


7=1 


(47) 


FMae  is  the  feature  calculated  from  mean  absolute  error  and  the  FMse  is  calculated  from  mean 
square  error.  In  Equations  46  and  47,  N  is  the  length  of  the  mixed  1-D  data  vector,  M  is  the 
total  number  of  non-defective  rivets.  LD[i]  is  the  /'th  pixel  of  mixed  data  for  different 
frequencies  and  Lndj[}]  is  set  to  zero  vector  as  the  baseline  signal.  The  line  scans  at  different 
frequencies  for  the  0.3  inch  subsurface  crack  signals  are  shown  in  Figure  76. 
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x  10-9  Pa  rameter-FrequencyfO  DA  subtract  ion) 


Figure  76  Simulated  Line  Scans  at  Different  Frequencies  for  0.3-inch  Subsurface  Cracks 

after  ODA  Processing 

The  SNR  definition  used  earlier  is  modified.  The  definition  used  in  analyzing 
experimental  data  is  provided  in  Equation  (48).  The  modified  equation  for  analysis  of 
simulated  data  is  given  in  (49),  where  the  mean  and  variance  are  set  to  0  and  1 ,  respectively, 


SNR 


JM 

•  (49) 


where  M=  2  in  this  case. 

A  scatter  plot  of  two-dimensional  feature  vector  in  the  feature  space  is  shown  in  Figure 
77.  A  classification  rule  can  be  devised  based  on  partitioning  the  defect  free  fastener  (green 
dot)  from  the  feature  vector  corresponding  to  the  0.3  inch  crack  data  (red  triangle).  The  SNR 
value  is  proportional  to  Euclidean  distance  between  (0,0)  and  the  feature  at  each  frequency. 
The  SNR  versus  frequency  plots  are  shown  in  Figure  78  and  it  is  seen  that  the  SNR  maximum 
is  at  400  Hz  frequency. 


M 


(Pi 


M 


i= 1 


(48) 
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182-D  feature  space  for  SIM  data  at  diff.  freq. 
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Figure  77  Two-Dimensional  Feature  Space  for  Simulated  Data  at  Different  Frequencies 
x10-9SNR  vs.  Frequency  for  Simulated  Data 
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Figure  78  SNR  versus  Frequency  for  Simulated  Signal  for  S-2  Sample  with  0.3  inch 

Subsurface  Crack 

From  Figure  77,  the  model  predicted  optimum  frequency  is  400  Hz,  which  is  close  to  the 
optimum  frequency  of  450  Hz  obtained  from  the  experimental  signals  for  S-2  sample. 

4.2.2. 3  Sensor  Liftoff  Effect 


The  parametric  effect  of  sensor  liftoff  on  the  signals  was  studied.  Liftoff  values 
considered  were  0.0050,  0.0095,  0.0150,  0.1,  0.15,  and  0.2  inches.  The  geometry  and  finite 


92 


element  mesh  of  the  S2  sample  and  sensor  configuration  are  shown  in  Figure  61  and  Figure 
62.  A  0.3-inch  second  layer  defect  is  introduced  under  the  fastener.  The  normal  component 
of  magnetic  flux  density  Bz  is  calculated  at  excitation  frequency  of  400  Hz.  The  real, 
imaginary  and  magnitude  of  the  normal  component  of  magnetic  flux  density  images  are 
plotted  for  each  individual  liftoff  value  in  Figure  79  through  Figure  84. 


(a)  (b)  (c) 

Figure  79  Bz  Plots  for  0.0050-inch  Liftoff  (a)  Real,  (b)  Imaginary,  (c)  Magnitude 


(a)  (b)  (c) 

Figure  80  Bz  Plots  for  0.0095-inch  Liftoff  (a)  Real,  (b)  Imaginary,  (c)  Magnitude 


(a)  (b)  (c) 

Figure  81  Bz  Plots  for  0.0150-inch  Liftoff  (a)  Real,  (b)  Imaginary,  (c)  Magnitude 
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(a)  (b)  (c) 

Figure  82  Bz  Plots  for  0.1 -inch  Liftoff  (a)  Real,  (b)  Imaginary,  (c)  Magnitude 


(a)  (b)  (c) 

Figure  83  Bz  Plots  for  0.15-inch  Liftoff  (a)  Real,  (b)  Imaginary,  (c)  Magnitude 


(a)  (b)  (c) 

Figure  84  Bz  Plots  for  0.2-inch  Liftoff  (a)  Real,  (b)  Imaginary,  (c)  Magnitude 


To  see  the  variation  in  the  signals  with  sensor  liftoff,  the  peak  to  peak  values  in  line  scans 
across  the  center  of  the  fastener  are  considered.  The  real  and  imaginary  parts  of  the  line  scans 
are  plotted  in  Figure  85  (Real)  and  Figure  86  (Imaginary).  The  real  and  imaginary  data  are 
demodulated  using  the  optimum  detection  angle  of  70  degrees  at  400  Hz  and  the  resulting 
mixed  signals  are  shown  in  Figure  87.  The  peak  to  peak  values  of  Bz  (real,  imaginary, 
magnitude  and  mixed)  at  each  liftoff  are  listed  in  Table  1 1  and  plotted  in  Figure  88. 
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Sensor  Liftoff  Real  Line  Plot 


x  10'8 


xIO'4 


Figure  85  Line  Scans  across  the  Center  of  Fastener  of  Real  Part 


x  io’9  Sensor  Liftoff  Imaginary  Line  Plot 


Figure  86  Line  Scans  across  the  Center  of  Fastener  of  Imaginary  Part 
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x  10 9  Sensor  Liftoff  Mixed  Line  Plot  using  ODA 


Figure  87  Mixed  Line  Scans  Using  ODA 


Table  1 1  Peak-to-Peak  Values  of  Bz  for  Various  Sensor  Liftoffs 


Sensor  liftoff 
(inch) 

Real 

Imaginary 

Mixed 

0.005 

3.3484e-008 

1.3586e-008 

6.31 19e-009 

0.0095 

3.1232e-008 

1.2620e-008 

5.9523e-009 

0.015 

2.8929e-008 

1.1657e-008 

5.5312e-009 

0.053 

1.9075e-008 

7.6667e-009 

3.2867e-009 

0.1 

1.222  le-008 

5.3627e-009 

1.5625e-009 

0.15 

7.7481e-009 

4.5275e-009 

7.3745e-010 

0.2 

4.7604e-009 

4.2386e-009 

1.7312e-009 
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Peak  to  Peak  Value  us  Sensor  Liftoff 


Real 

— ■ — Imaginary 
Mixed 


Figure  88  Peak  Values  of  Real,  Imaginary,  and  Mixed  MR  Signals  Versus  Liftoff 


As  expected,  when  the  sensor  liftoff  increases,  the  measured  signal  decreases.  Since  the 
optimum  detection  angle  is  selected  based  on  minimum  peak-to-peak  value  in  the  image  data 
of  the  fastener,  it  is  seen  that  in  the  mixed  signals  the  overall  spread  in  the  peak  value  with 
liftoff  is  also  lower;  hence,  the  change  in  signal  values  with  liftoff  is  relatively  flatter  than  that 
of  raw  data  reported. 


4.2.2.4  Conductivity  of  Top  Laver  Effect 

The  conductivity  values  for  the  top  layer  were  selected  according  to  the  Parametric  Table 
presented  earlier.  The  operation  frequency  is  chosen  as  400  Hz  and  liftoff  is  chosen  as  0.0050 
inch.  The  geometry  of  the  S-2  standard  including  a  fastener  with  0.3  inch  defect  on  the 
bottom  layer  is  modeled.  The  nominal  conductivity  for  top  layer  plate  for  S-2  sample  is  29.6 
percent  IACS,  and  the  range  of  conductivity  of  the  top  layer  plate  is  assigned  values  from  28 
percent  to  33  percent  IACS.  Figure  89  through  Figure  94  present  simulation  results  of  real, 
imaginary,  and  magnitude  value  of  the  normal  component  of  magnetic  flux  density  for 
different  values  of  top  layer  conductivity.  Table  12  shows  the  peak  value  of  the  signal 
magnitude  at  different  values  of  conductivity. 
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(a) 


(b) 


(c) 


Figure  89  Bz  with  28%  lacs  Top  Layer  Conductivity  (a)  Real,  (b)  Imaginary,  (c) 

Magnitude 


(a) 


(b) 


(c) 


Figure  91  Bz  with  30%  lacs  Top  Layer  Conductivity  (a)  Real,  (b)  Imaginary,  (c) 

Magnitude 
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(a) 


(b) 


(c) 


Figure  92  Bz  with  31%  lacs  Top  Layer  Conductivity  (a)  Real,  (b)  Imaginary,  (c) 

Magnitude 


(a)  (b)  (c) 

Figure  93  Bz  with  32%  lacs  Top  Layer  Conductivity  (a)  Real,  (b)  Imaginary,  (c) 

Magnitude 


(a)  (b)  (c) 

Figure  94  Bz  with  33%  lacs  Top  Layer  Conductivity  (a)  Real,  (b)  Imaginary,  (c) 

Magnitude 
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Table  12  Peak  Values  of  Magnitude  of  Bz  for  Various  Top  Layer  Conductivities 


Conductivity 
(%  IACS) 

Conductivity 

(Siemens/m) 

Peak  value  of  the  signal  magnitude  (Tesla) 

28 

1.6240e7 

2.0141e-8 

29 

1.6820e7 

2.0150e-8 

30 

1.7391e7 

2.0206e-8 

31 

1.7980e7 

2.0312e-8 

32 

1.8560e7 

2.0358e-8 

33 

1.9140e7 

2.04 19e-8(  1.4%  difference  vs.  28%  IACS) 

Figure  95  shows  a  plot  of  peak  value  of  magnitude  signal  magnitude  versus  top  layer 
plate  conductivity. 


x  10  T'eak  Value  of  Signal  Magnitude  vs.  Tap  Layer  Plate  Conductivity 


Figure  95  Peak  Value  of  the  Signal  Magnitude  Versus  Top  Layer  Plate  Conductivity 

From  the  results  in  Table  12  and  the  Figure  95,  it  is  seen  that  the  peak  value  of  the  signal 
increases  slightly  when  top  layer  plate  conductivity  increases.  This  is  as  expected  since 
increase  in  the  top  layer  plate  conductivity  increases  the  induced  eddy  currents  thereby 
increasing  the  amplitude  of  the  measured  signal. 
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4.2.2. 5  Conductivity  of  Bottom  Laver  Effect 


The  conductivity  values  for  the  bottom  layer  were  selected  according  to  the  Parametric 
Table  presented  earlier.  For  bottom  layer  plate,  the  conductivity  value  ranges  from  30  percent 
IACS  to  36  percent  IACS.  The  operation  frequency  is  chosen  as  400  Hz  and  liftoff  is  chosen 
as  0.0050  inch.  The  geometry  of  the  S-2  standard  including  a  fastener  with  0.3  inch  defect  on 
the  bottom  layer  is  modeled.  The  nominal  conductivity  for  bottom  layer  in  the  S-2  standard  is 
33.5%  IACS,  and  the  range  of  the  bottom  layer  conductivity  is  chosen  to  vary  from  30% 

IACS  to  36%  IACS.  Figure  96  through  Figure  102  present  simulation  results  of  real, 
imaginary,  and  magnitude  value  of  the  normal  component  of  magnetic  flux  density  for 
different  values  of  bottom  layer  conductivity.  Table  13  shows  the  peak  value  of  the  signal 
magnitude  at  different  values  of  conductivity.  Figure  103  shows  the  plot  of  peak  value  of  the 
signal  magnitude  versus  bottom  layer  plate  conductivity. 


(a) 


(b) 


(c) 


Figure  97  Bz  with  31%  lacs  Bottom  Layer  Conductivity  (a)  Real,  (b)  Imaginary,  (c) 

Magnitude 
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Figure  100  Bz  with  34%  lacs  Bottom  Layer  Conductivity  (a)  Real,  (b)  Imaginary,  (c) 

Magnitude 
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(a) 


(b) 


(c) 


Figure  101  Bz  with  35%  lacs  Bottom  Layer  Conductivity  (a)  Real,  (b)  Imaginary,  (c) 

Magnitude 


Table  13  Peak  Values  of  Bz  Magnitude  for  Various  Bottom  Layer  Conductivities 


Conductivity  (% 
IACS) 

Conductivity 

(Siemens/m) 

Peak  value  of  the  signal  magnitude  (Tesla) 

30 

1.7391e7 

2.0390e-8 

31 

1.7980e7 

2.0346e-8 

32 

1.8560e7 

2.0271e-8 

33 

1.9140e7 

2.0216e-8 

34 

1.97 10e7 

2.0153e-8 

35 

2.0290e7 

2.0091e-8 

36 

2.0869e7 

1.9867e-8  (2.5%  difference  versus  30% 
IACS) 
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Figure  103  Peak  Value  of  the  Signal  Magnitude  Versus  Bottom  Layer  Plate  Conductivity 

From  the  results  presented  above,  it  is  seen  that  the  peak  value  of  the  signal  decreases 
when  bottom  layer  conductivity  increases.  When  the  bottom  layer  plate  conductivity 
increases,  the  induced  eddy  current  becomes  more  concentrated  around  the  bottom  surface 
and  the  observed  signal  on  the  top  surface  is  reduced. 

4.2.2. 6  Conductivity  of  Fastener  Effect 


The  conductivity  of  fastener  was  varied  as  described  in  the  parameter  table.  The 
geometry  of  the  S2  sample  with  a  0.3  inch  crack  and  finite  element  mesh  are  as  shown  in 
Figure  56  and  Figure  57.  The  operation  frequency  is  chosen  as  400  Hz  and  liftoff  is  chosen  as 
0.0050  inch.  The  different  conductivity  values  for  the  fastener  were  chosen  to  be  1.0  percent 
IACS,  2.2  percent  IACS,  and  3.1  percent  IACS.  The  real,  imaginary,  and  magnitude  values 
of  magnetic  flux  density  at  the  MR  sensor  for  different  conductivity  values  are  plotted  in 
Figure  104. 
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Figure  104  Real,  Imaginary,  and  Magnitude  of  Magnetic  Flux  Density  for  Conductivity 
Values  of  Ti  Fastener  (1.0  percent  lacs,  2.2  percent  lacs,  and  3.1  percent  lacs)  as  Indicated  on 

the  Left 

The  variation  of  the  peak  values  are  summarized  in  Table  14.  The  fastener  conductivities 
considered  in  this  study  does  not  affect  the  signals  with  any  significance. 


Table  14  Effect  of  Fastener  Conductivity  on  Peak  Value  of  Magnitude  Signal 


Fastener 

Conductivity 

1.0 

2.2 

3.1 

Peak  values 
of  magnitude 

2.0185e-008 

2.0212e-008 

2.0187e-008 

4.2.2. 7  Fastener  to  Edge  Effect 

A  commonly  encountered  problem  in  airframe  geometry  is  the  influence  of  surface  and 
subsurface  edge  on  fastener  and  crack  signals.  The  edge  discontinuity  behaves  as  a  large 
defect  and  generates  its  own  signature  that  can  affect  the  defect  signal  and  thereby  lead  to 
misinterpretation  of  MR  data.  This  section  describes  a  systematic  study  of  the  effect  of 
fastener-to-edge  distance  on  the  defect  signal.  The  S-2  sample  fastener  with  a  0.3  inch, 
second  layer  crack  geometry  presented  in  section  4.2.1  is  considered.  The  edge  is  in  the 
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second  layer  and  the  nominal  fastener-to-edge  distance  for  S-2  sample  is  0.6  inch.  The  range 
of  values  for  parametric  variations  is  0.4  to  0.8  inch. 

Simulations  were  performed  at  400  Hz  excitation  frequency  and  0.0095  inch  liftoff.  The 
normal  component  of  the  magnetic  flux  density  Bz  is  calculated  at  the  prescribed  liftoff  and 
the  resulting  image  data  are  plotted  for  each  simulation.  The  real,  imaginary,  magnitude  and 
demodulated  images  are  presented.  The  demodulated  image  is  derived  using  optimum 
detection  angle  (ODA)  which  is  70  degrees.  Figure  105  to  Figure  109  show  the  real, 
imaginary,  magnitude  and  demodulated  images  for  varying  fastener-to-edge  distance  in  the 
range  0.4  to  0.8  inch. 


Figure  105  Real,  Imaginary,  Magnitude,  Demodulated  Bz  for  Fastener  Edge  Distance  = 

0.4  inch 
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Figure  106  Real,  Imaginary,  Magnitude,  Demodulated  Bz  for  Fastener  Edge  Distance 

0.5  inch 


Figure  107  Real,  Imaginary,  Magnitude,  Demodulated  Bz  for  Fastener  Edge  Distance 

0.6  inch 
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Figure  108  Real,  Imaginary,  Magnitude,  Demodulated  Bz  for  Fastener  Edge  Distance  = 

0.7  inch 


Real 

4  * 

Imaginary 

< 

IHI 1 

Magnitude 

Demodulated 

Figure  109  Real,  Imaginary,  Magnitude,  Demodulated  Bz  for  Fastener  Edge  Distance  = 

0.8  inch 

The  line  plots  are  extracted  from  all  the  image  data  and  plotted.  Figure  110  and  Figure 
111  show  the  line  scans  across  the  center  of  the  rivet  image  of  real  and  imaginary  parts. 
Figure  112,  shows  the  line  scan  of  the  demodulated  signal  using  ODA.  When  the  fastener-to- 
edge  distance  increases,  the  edge  effect  on  the  right  side  of  the  fastener  is  decreased,  so  the 
magnitude  of  the  signal  on  the  right  lobe  will  decrease  slightly. 
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X  Distance(m) 


Figure  110  Line  Scans  Across  the  Center  of  the  Fastener  for  Real  Part  of  Bz 


X  Distance(m) 


Figure  111  Line  Scans  Across  the  Center  of  the  Fastener  for  Imaginary  Part  of  B 
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Parameter-Fastener-edge-dist  (ODA, bottom  0.3"  crack) 


-0.02  -0.015  -0.01  -0.005  0  0.005  0.01  0.015  0.02 

X  Distance(m) 


Figure  1 12  Line  Scans  Across  the  Center  of  the  Fastener  for  Mixed  Signal  Using  ODA 

To  further  examine  the  effect  of  edge  on  defect  signal  amplitude,  the  signal  from  a 
defective  fastener  for  the  largest  fastener-to-edge  distance  is  calculated  and  subtracted  from 
the  defect  signal,  in  each  case,  at  the  defect  location.  Assuming  that  there  is  minimal  effect 
on  defect  signal  when  the  edge  is  the  farthest,  this  value  reflects  the  effect  of  an  edge  in  the 
proximity  of  the  fastener.  The  defect  signal  amplitude  calculated  in  this  manner  is 
summarized  in  Table  15  and  plotted  in  Figure  1 13. 


Figure  1 13  Edge  Effect  on  Defect  Signal  Amplitude:  Real,  Imaginary,  and  Mixed  Signals 

Versus  Fastener-to-Edge  Distance 
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Table  15  The  Peak-to-Peak  Values  Versus  Fastener-to-Edge  Distance 


Fastener-edge  Distance  (inch) 

0.4 

0.5 

0.6 

0.7 

OO 

O 

Real  signal 

-0.2593  e-7 

0.1501  e-7 

0.2031  e-7 

0.1908  e-7 

0 

Imag  signal 

0.3132  e-6 

0.2369  e-6 

0.1898  e-6 

0.0810  e-6 

0 

Mixed  signal 

0.2430  e-6 

0.2326  e-6 

0.2352  e-6 

0.0936  e-6 

0 

From  Figure  1 13,  it  is  observed  that  the  closer  the  fastener  is  to  an  edge  the  more  effect 
on  defect  signal.  The  signal  due  to  the  edge  adds  to  the  defect  signal  which  increases  with 
proximity  of  fastener  to  the  edge. 

4.2.2.8  Crack  Dimension  Effect 


Based  on  the  S-2  flaw  schematics  diagram,  four  different  trapezoidal  crack  dimensions 
were  modeled  including  the  case  of  a  crack- free  fastener.  The  radial  dimensions  of  the  top  of 
crack  top  were  chosen  to  be  0.2,  0.22,  0.25  and  0.3  inch.  The  crack  geometry  was  tapered 
towards  the  bottom  surface.  The  sensor  liftoff  was  kept  at  0.0095  inch  and  the  frequency  was 
chosen  to  be  500  Hz.  The  normal  component  of  the  magnetic  flux  density  ( Bz )  is  calculated 
and  plotted.  The  real,  imaginary,  and  magnitude  of  the  complex  flux  density  are  plotted  in 
Figure  114  through  Figure  118. 


(a)  (b)  (c) 

Figure  1 14  Bz  for  0.2  inch  Crack  (a)  Real,  (b)  Imaginary,  (c)  Magnitude 
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(a)  (b)  (c) 

Figure  115  Bz  for  0.22  inch  Crack  (a)  Real,  (b)  Imaginary,  (c)  Magnitude 


(a)  (b)  (c) 

Figure  116  Bz  for  0.25  inch  Crack  (a)  Real,  (b)  Imaginary,  (c)  Magnitude 


(a)  (b)  (c) 

Figure  1 17  Bz  for  0.3  inch  Crack  (a)  Real,  (b)  Imaginary,  (c)  Magnitude 


112 


(a)  (b)  (c) 

Figure  118  Bz  for  no  Crack  (a)  Real,  (b)  Imaginary,  (c)  Magnitude 


In  order  to  see  the  variation  in  the  signals  with  increasing  crack  dimensions,  the  line 
scans  across  the  center  of  the  fastener  image  data  are  plotted.  The  real  part  of  the  magnetic 
flux  density  is  shown  in  Figure  119  and  Figure  120  is  the  image  of  the  imaginary  part.  The 
demodulated  signal  using  an  optimum  detection  angle  of  70  degree  at  500  Hz  is  shown  in 
Figure  121. 


1.5 

1 

0.5 

0 

-0.5 
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Figure  119  Line  Scans  across  the  Center  of  the  Fastener  for  Real  Image  Data 


Real  Line  Plot 


113 


x -lO'9  Imaginary  Line  Plot 


x  10'4 


Figure  120  Line  Scans  Across  the  Center  of  the  Fastener  For  Imaginary  Image  Data 
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Mixed  Line  Plot  using  ODA 


-4 

x  10 


Figure  121  Mixed  Line  Scans  Using  ODA  of  70  Degrees 

The  defect  signal  amplitude  is  calculated  using  the  peak-to-peak  value  of  all  3  signals, 
real,  imaginary,  and  mixed.  The  values  calculated  are  listed  in  Table  16  and  plotted  in  Figure 
122.  The  defect  signal  values  are  also  recalculated  after  subtracting  the  value  obtained  for  the 
crack-free  fastener.  These  values  are  summarized  in  and  plotted  in  Figure  123. 

Table  16  The  Peak-to-Peak  Values  of  Bz  for  Various  Crack  Dimensions 


Crack  width  on  top 
(inch) 

Real 

Imaginary 

Mixed 

0 

2.8433e-008 

1.4821e-008 

2.5146e-009 

0.2 

2.9422e-008 

1.4256e-008 

2.9636e-009 

0.22 

2.9683e-008 

1.3963e-008 

3.2886e-009 

0.25 

2.9941e-008 

1.3622e-008 

3.6405e-009 

0.3 

3.1045e-008 

1.2658e-008 

5.5201e-009 
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Peak-to-peak  Values  vs.  Crack  Dimension 


-9 

xIO 


Figure  122  Peak-to-peak  values  of  Real,  Imaginary,  and  Mixed  MR  Signals  Versus  Crack 

Area 

Table  17  The  Peak-to-Peak  Values  of  Bz  for  Various  Crack  Dimensions  (after  Subtracting 
the  Value  Obtained  for  the  Crack-Free  Fastener) 


Crack  width  on  top 
(inch) 

Real 

Imaginary 

Mixed 

0.2 

9.89e-010 

-5.65e-010 

4.49e-010 

0.22 

1.25e-009 

-8.58e-010 

7.74e-010 

0.25 

1.508e-009 

-1.199e-009 

1.1259e-009 

0.3 

2.612e-009 

-2.163e-009 

3.0055e-009 
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Peak  to  Peak  Value  vs  Crack  Dimension 


-  Real 

-  Imaginary 
Mixed 


Figure  123  Peak-to-Peak  Values  of  Real,  Imaginary,  and  Mixed  MR  Signals  Versus 

Crack  Area 


4.3  OPTIMIZED  SCANNING 


The  first  set  of  scans  of  the  S-2  sample,  shown  through  Sections  4. 1 . 1 .2  and  4. 1 . 1 .3,  were 
made  as  a  benchmark  from  which  the  modeling  calculations  could  be  validated,  as  was  shown 
in  Section  4.2.1.  The  frequency  parametric  study,  shown  in  Section  4.2.2.2,  indicated  at  one 
point  that  the  frequency  ought  to  be  increased  beyond  the  maximum  500  Hz  used  in  this  first 
set  of  scans.  Indeed,  as  indicated  in  the  above  description  of  these  scans  and  shown  in  Figs. 

4. 1.2-4  through  4. 1.2-7,  the  differentiation  between  fasteners  with  and  without  cracks  appears 
better  with  higher  frequencies,  but  a  visual  inspection  of  the  various  c-scans  did  not  easily 
yield  a  quantitative  justification  of  that  statement. 

This  modeling  result  of  higher  frequency  prompted  a  new  set  of  scans  to  be  made.  All 
four  of  the  test  standards:  S-l,  S-2,  S-3  and  S-4,  were  rescanned  over  the  same  frequency 
range  as  before, /=  100  to  500  Hz  at  steps  of  50  Hz.  Beyond  this,  the  S-2  standard  was 
scanned  from  550  to  800  Hz  in  steps  of  50  Hz.  In  addition,  each  of  the  scans  was  made  five 
times.  This  was  done  in  support  of  the  advanced  feature  recognition  algorithm  development 
to  provide  multiple  sets  of  data  for  statistical  analysis. 

4.3.1  Performance  on  New  data  from  S2  sample 

The  second  dataset  provided  by  Boeing  includes  data  at  15  frequencies  from  100  Hz  to 
800  Hz  in  steps  of  50  Hz.  For  each  frequency,  the  MR  data  was  acquired  five  times  using  five 
independent  scans.  For  instance,  the  Ti  fastener  images  of  S-2  sample  at  400  Hz  and  650  Hz 
for  five  scans  are  shown  in  Figure  124  and  Figure  125  with  optimum  detection  angles  of  70 
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and  80  degrees,  respectively.  As  we  can  see,  the  outcomes  of  five  different  scans  are  almost 
identical  to  each  other  in  two-dimensional  images. 


(a)  (b)  (c)  (d)  (e) 

Figure  124  Titanium  Fastener  Images  for  S-2  Test  Standard  at  400  Hz  with  ODA  =70 
degrees  in  Five  Experimental  Scans  (a)  to  (e) 


118 


Figure  125  Titanium  Fastener  Images  for  S-2  Test  Standard  at  650  Hz  with  ODA  =  80 
degrees  in  Five  Experimental  Scans  (a)  to  (e) 


The  data  acquired  on  the  S-2  standard  sample  with  Titanium  fasteners  at  different 
frequencies  was  analyzed.  Figure  126  and  Figure  127  present  the  results  in  the  form  of  line 
plots  and  scatter  plots  in  feature  space,  respectively. 
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x§=2  Freq-400  Inside  Ti:ODA-70  Scan:1  2  3  4  5 
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S-2  Freq-400  Outside  Ti:ODA-70  Scan:  1  2  3  4  5 
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(e) 


S-2  Freq-600  Inside  Ti:ODA-80  Scan:1  2  3  4  5 


S-2  Freq-600  Outside  Ti:ODA-70  Scan:  1  2  3  4  5 


(f) 
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Sigal(Line  plots)  Slgal(Llne  plots) 


(h) 

Figure  126  Collection  of  all  Line  Scans  for  Defect-Free  (blue)  and  Defective  (red) 

Fasteners 


Note:  Outside  row  (right)  and  inside  row  (left):  (a)  to  (h)  for  frequencies  100  Hz  to  800  Hz 
in  steps  of  100  Hz 
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(a) 


JStJ2  Freq-200  Outside  Ti:ODA-50  Scan:  1  2  3  4  5 


(b) 
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(d) 


Freq-400  Outside  Ti:ODA-70  Scan:  1  2  3  4  5 
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Freq-700  Outside  Ti:ODA-80  Scan:1  2  3  4  5 


(g) 


Freq-SOO  Outside  Ti:ODA-SO  Scan:  1  2  3  4  5 


(h) 

Figure  127  Scatter  Plot  of  Features-Defect-Free  (green),  Defective  (red) 

Note:  Outside  row  (right)  and  Inside  row  (left):  (a)  to  (h)  represent  frequencies  100  Hz  to 
800  Hz  in  steps  of  100  Hz 


The  SNR  values  calculated  for  the  S-2  sample  data  from  inspection  of  both,  outside  and 
inside  fasteners,  at  various  frequencies  are  summarized  in  Table  18  and  Table  19. 
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Table  18  SNR  Values  for  S-2  Sample  -  Outside  Fasteners  at  Different  Frequencies 

OUTSIDE 


Frequency(Hz) 

0.20  inch 

0.22  inch 

0.25  inch 

0.30  inch 

100 

4.5467 

0.8736 

6.4469 

69.3307 

150 

6.7753 

3.8642 

9.9404 

85.8019 

200 

7.5736 

4.2142 

9.6473 

76.8161 

250 

8.0441 

4.9579 

13.4174 

104.7986 

300 

15.6455 

9.822 

22.0077 

161.8626 

350 

8.2825 

6.293 

18.1838 

107.7558 

400 

15.5946 

9.8186 

26.3239 

152.1197 

450 

22.808 

15.3956 

34.3392 

175.6716 

500 

12.9118 

9.4311 

18.7284 

80.1985 

550 

1 1 .645 

10.3024 

26.6577 

145.4339 

600 

7.134 

7.499 

20.3077 

111.1863 

650 

6.746 

4.3718 

9.9863 

54.4714 

700 

4.0922 

2.382 

6.0518 

35.148 

750 

5.2685 

4.012 

8.1053 

38.6236 

800 

5.2444 

4.0292 

5.9675 

30.0911 
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Table  19  SNR  Values  for  S-2  Sample  -  Inside  Fasteners  at  Various  Frequencies 

INSIDE 


Frequency(Hz) 

0.20  inch 

0.22  inch 

0.25  inch 

0.30  inch 

100 

0.7722 

0.5452 

3.2376 

16.4065 

150 

1.5302 

1.7715 

5.5316 

28.6856 

200 

1.3708 

1.7229 

6.1701 

32.6099 

250 

1.3591 

1 .6483 

6.7992 

36.843 

300 

2.3389 

2.5246 

8.3488 

44.8557 

350 

1.5757 

1.5825 

6.915 

40.3348 

400 

2.1199 

2.1575 

7.6973 

41.719 

450 

3.4818 

3.1604 

9.8738 

46.4146 

500 

3.4771 

3.3807 

8.9754 

45.0073 

550 

6.7281 

7.4048 

25.413 

116.7166 

600 

6.4248 

3.798 

12.2827 

52.7032 

650 

5.8295 

5.6663 

12.4962 

54.6064 

700 

8.7571 

4.4901 

12.1609 

57.5303 

750 

5.8697 

3.73 

9.124 

41.7375 

800 

5.8722 

4.9402 

9.7067 

45.3913 
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Figure  128  Scatter  Plot  of  Features  Defect-Free  (green),  Defective  (red) 

Note:  Outside  row  (right)  and  Inside  row  (left)  at  the  optimum  frequencies:  450  Hz  and  550 
Hz,  respectively 

From  the  analysis  results  of  experimental  data  at  different  frequencies,  one  can  generate  a 
plot  of  SNR  versus  frequency  as  presented  in  Figure  129  and  Figure  130,  for  outside  and 
inside  fasteners.  These  results  show  that  the  optimum  frequency  that  maximizes  the  SNR 
criteria  is  450  Hz  for  outside  fasteners  and  that  for  inside  fasteners  is  550  Hz. 


Figure  129  SNR  (Raw  and  Curve-Fitted)  Plots  Versus  Frequency  for  S-2  Outside 

Fasteners 
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SNR  vs  Frequency  for  S-2  INSIDE  Rivets 


SNR  vs  Frequency  for  S-2  INSIDE  Rivets 


Curve-Fitted  SNR  vs  Frequency  for  S-2  INSIDE  Rivets 


Curve- Fitted  SNR  vs  Frequency  for  S-2  INSIDE  Rivets 


Figure  130  SNR  (Raw  and  Curve-Fitted)  Plots  Versus  Frequency  for  S-2  Inside  Fasteners 


4.3.2  Optimized  Scan  Conclusions 

As  was  mentioned  in  the  previous  section,  a  full,  comprehensive  set  of  scans  was  made  of 
the  S-2  test  standard  that  spanned  frequencies  from  100  to  800  Hz.  This  span  was  larger  than 
the  original  set,  shown  in  Section  4. 1.3.2.  Based  on  earlier  work  it  was  suggested  that 
optimum  frequency  would  be  found  higher  than  the  500  Hz  that  was  scanned  at  that  time. 
Subjectively,  it  was  determined  that  400-500  Hz  was  the  best  scan  frequency  range,  but  that 
was  based  more  upon  the  visual  interpretation  of  the  c-scans  rather  than  any  quantitative 
determination. 


It  was  shown  in  Figure  129  and  130  that  this  new  set  was  indeed  optimized  close  to  that 
range,  now  based  upon  contour  features  quantitatively  determined  through  equations  rather 
than  the  subjective  look  of  scan  plots.  Both  the  modeling  with  frequency  parameter  and  the 
application  of  advanced  feature  recognition  algorithms  agree  with  the  more  anecdotal 
conclusion  of  this  best  frequency  range. 


4.4  TRANSITION  TO  MAUS 


As  was  mentioned  briefly  in  Section  4. 1.1.2,  an  MR  Sensor  module  is  currently  being 
added  to  the  MAUS  system  in  Boeing  St.  Louis.  The  findings  of  this  study  have  been  made 
available  to  that  program,  and  in  particular  the  feature  recognition  algorithms  presented 
throughout  Section  4.1.3  will  be  taken  into  consideration  for  incorporating  automatic  feature 
detection  with  the  MAUS  system. 
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As  this  new  field  deployable  system  comes  online,  there  are  interesting  and  pertinent 
findings  and  conclusions  that  have  been  reported.  One  of  the  more  impactful  insights  to  this 
kind  of  modeling  and  laboratory  testing  effort  involves  the  comparison  between  fabricated 
test  standards  and  actual  airframe  structures  that  are  being  scanned.  As  was  shown  in  Figure 
18,  the  MR/MAUS  system  has  been  used  to  scan  an  actual  wing.  There  are  two  points  to 
note: 

First,  the  arrangement  of  steel  and  titanium  fasteners,  such  as  seen  in  the  S-2  test  standard 
that  was  so  prominently  studied  in  this  program,  is  not  a  realistic  example  of  what  is  actually 
encountered  in  the  field  on  real  structures.  Rather,  there  is  typically  a  non-ordered  collection 
of  many  kinds  of  fasteners,  particularly  on  older  aircraft.  This  fact  may  itself  not  be  of  any 
immediate  concern  when  making  more  fundamental  developments  on  standards  such  as  the  S- 
2;  however,  it  is  important  to  note  that  the  different  materials  are  not  usually  so  isolated. 

Second,  on  typical  airframe  structures  such  as  the  wings  that  have  been  scanned  using  the 
MR/MAUS  system,  there  is  no  bottom  layer  edge  encountered  that  is  so  close  to  the  fasteners 
that  it  causes  the  kind  of  interference  seen  in  the  test  measurements  and  considered  in  the 
modeling  work.  The  absence  of  such  an  edge  is  welcome,  of  course,  in  that  it  tends  to 
overwhelm  the  smaller  MR  sensor  response  from  cracks  parallel  to  and  near  the  edge.  The 
report  from  the  MR/MAUS  effort  is  that  the  kind  of  bottom  layer  geometry  that  comprises  the 
S-2  test  standard  is  simply  not  a  realistic  representation  of  the  kind  of  structures  being 
encountered  in  the  field.  This  implies  that  the  test  standard  arrangements  themselves  may  be 
driving  unnecessary  efforts  in  discerning  cracks  near  edges  when  this  situation  may  indeed 
not  be  common,  at  least  in  the  MR/MAUS  scans  presently  under  consideration. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  results  and  findings  of  Advanced  Nondestructive  Evaluation  Sensor  Modeling  for 
Multi-Site  Inspection  have  been  presented.  A  FEM  scheme  was  developed  to  accurately 
represent  the  fields  and  response  of  MR  sensors  from  AC  currents  induced  in  complex  metal 
geometries  that  are  of  interest  to  the  Air  Force,  namely,  cracks  of  various  shapes  and  sizes  in 
the  vicinity  of  fasteners  that  splice  metal  layers  together  as  on  an  airframe.  This  model  was 
developed  in  calibration  with  measurements  made  in  a  laboratory  environment.  Once 
satisfactorily  calibrated,  the  model  was  used  to  conduct  a  variety  of  parametric  studies 
including  frequency,  sensor  liftoff,  conductivity  of  layers,  fastener  to  edge  distance  and  crack 
dimension. 

In  addition  to  the  modeling  efforts  presented,  the  results  of  meticulous  analysis  of  the 
laboratory  measurements  were  presented  as  part  of  the  development  of  automated  feature 
recognition.  The  data  was  further  processed  to  reduce  the  results  to  figures  of  merit,  various 
feature  values  that  incorporated  both  magnitude  and  shape  of  signal  and  ultimately  SNRs.  It 
has  been  noted  that  applying  meaningful  quantitative  conclusions  to  c-scan  results  through 
subjective  observation  is  challenging  and  not  necessarily  repeatable,  particularly  in  borderline 
cases  of  deep  and/or  small  cracks.  These  automated  algorithms  have  served  to  place  an 
objective  value  to  cracks  as  they  appear  in  MR  sensor  data.  As  a  result  of  quantitative  data 
reduction  it  has  been  possible  to  establish  a  signal  to  noise  ratio  value  for  the  MR  sensing. 
With  this  evaluation  it  was  possible  to  show  that  MR  sensors  offer  considerable  value  over 
conventional  eddy  current  inspections  with  the  possibility  of  improvements  in  sensitivity  by 
over  an  order  of  magnitude. 

The  results  of  this  study  lead  to  the  following  recommendations: 

1 .  MR  Sensor  Implementation  -  The  improved  detection  capability  of  MR  sensors 
measured  in  this  program  indicated  that  implementation  under  the  MR  MAUS  program  is 
well  justified.  However,  additional  operational  schemes  should  also  be  considered  to 
obtain  the  potential  value  of  this  inspection  method. 

2.  Reverse  Modeling  -  This  program  has  consisted  primarily  of  forward  modeling,  in  which 
a  known  geometry,  arrangement  of  cracks,  etc.  is  modeled  in  order  to  obtain  sensor 
response  to  some  parametric  variation  such  as  frequency.  A  more  challenging,  yet  more 
rewarding  effort  would  be  to  follow  with  reverse  modeling  that  would  ultimately  take  the 
results  and  backward  model  into  the  more  likely  scenario  that  is  being  encountered.  For 
example,  given  a  result,  what  are  the  most  likely  size  and  direction  of  the  underlying 
crack? 

3.  Expanded  Parametric  Studies  -  There  are  additional  parametric  studies  that  need  be 
conducted,  partly  from  the  list  generated  through  this  program,  but  more  importantly 
through  feedback  from  the  MR/MAUS  program  or  in  response  to  current  Air  Force 
concerns. 

4.  Advanced  MR  Sensor  Development  -  This  work  has  utilized  straightforward  AC,  time- 
independent  excitation  signals.  There  is  interest  in  examining  pulsed  MR  Sensor  scanning 
in  which  time  response  is  incorporated  into  the  scan. 
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5.  Sensor  Orientation  -  This  work  has  concentrated  on  MR  sensing  using  Bz  measurements 
with  the  available  orientation  of  the  MR  devices.  The  modeling  results  indicate  that  Bx 
and  By  may  contain  additional  information  to  discriminate  important  characteristics  and 
should  be  investigate. 
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APPENDIX  A:  MATLAB  FUNCTION  ROUTINES  FOR  MR  SENSOR  DATA 
MANIPULATION 


A.l  INTRODUCTION 

Appendix  A  contains  a  description  of  Matlab  functions  created  by  Dr.  Ben  Koltenbah 
(Boeing)  for  reading  and  parsing  data  from  the  native  scan  data  format  (SCD  file  format)  of 
the  Test  Bed  Data  Acquisition  System  to  the  Matlab  environment.  These  routines  also  serve 
to  save  pertinent  data  in  the  more  standard  MAT  file  format,  which  can  be  loaded  into  a 
number  of  common  programming  environments  including  Mathematica.  Appendix  B 
contains  a  description  of  a  Mathematica  notebook  that  Dr.  Koltenbah  developed  for 
interactive  post-processing  of  MR  Sensor  data,  and  this  notebook  requires  data  files  to  first  be 
converted  from  SCD  to  MAT  using  the  Matlab  routines  that  are  described  here. 

These  functions  are  provided  “as  is,”  and  it  is  noted  that  these  are  not  required 
deliverables  on  this  or  any  other  program.  Dr.  Koltenbah  included  them  along  with  the  MR 
sensor  data  that  he  acquired  in  order  to  allow  all  collaborators  the  ability  to  quickly  read  the 
data  into  the  Matlab  environment  and  post-process  the  results.  The  routines  are  provided  in 
this  section  for  completeness. 

The  following  Matlab  functions  are  described  in  Appendix  A,  where  details  are  provided 
below: 


PostProcExtraction003  -  This  is  a  combined  function  that  reads  the  original 
SCD  file  format,  parses  the  data,  adds  some  common  post-processing  steps, 
and  saves  a  subset  of  the  data  into  the  MAT  file  format.  This  function  calls  the 
next  two  in  the  list. 

MRS_ReadSCDFile_01  -  This  function  reads  all  the  data  contents  from  an  SCD 
file  and  collects  them  into  data  structure  variables  within  the  Matlab 
environment. 

MRS_ExtractData_02  -  This  function  applies  common  post-processing  steps  to 
the  raw  data,  which  prepares  the  data  for  post-processing. 

MatlabExampleOOl  -  This  is  a  standalone  function  that  reads  a  MAT  data  file  and 
plots  the  results.  It  is  provided  as  an  example  to  show  basic  manipulation  of 
MR  sensor  data  in  the  Matlab  environment. 

A.2  POSTPROCEXTRACTION003 

As  stated  above,  the  MR  sensor  data  acquired  on  the  Test  Bed  System  is  recorded  in  a 
native  scan  data  format  as  files  with  SCD  extension.  The  function  PostProcExtraction003 
was  written  to  do  the  following  steps: 

1 .  Read  an  SCD  file 
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2.  Parse  the  data  into  structure  variables  within  the  Matlab  environment 

3.  Apply  common  post-processing  steps  to  convert  raw  data  into  calibrated 
component  data 

4.  Save  the  processed  data  into  a  MAT  file. 

Figure  A-l  shows  the  Matlab  code  for  this  function.  There  are  annotations  and 
comments  within  the  text,  displayed  as  green,  which  aid  in  following  the  procedural  steps  of 
the  code.  The  inputs  require  a  directory  name  and  file  name  (both  strings).  The  directory 
name  must  end  with  a  “\”  character.  The  file  name  is  an  SCD  file,  and  the  string  must  contain 
the  full  extension.  The  following  is  an  example  of  its  use  from  the  Matlab  command  prompt: 

>>  PostProcExtraction003 ( ' . \ ' , ' S-2  F400  ALL  06.scd'); 

In  this  example,  the  file  “S- 2  F400  ALL  06  .  scd”,  which  is  scan  data  of  the  S-2  test 
standard  acquired  at  400  Hz,  is  read,  parsed  and  processed,  and  then  a  subset  of  the  processed 
data  is  saved  as  “S-2  F400  ALL  06.  mat”.  There  is  no  text  output  generated. 

This  function  calls  the  two  aforementioned  subroutines,  which  are  highlighted  in  red 
within  the  code  listing  of  Figure  A-l.  These  will  be  detailed  next. 
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%  PostProcExtraction003  -  Convert  SCD  File  to  MAT  Format 

%  Dr.  Benjamin  E.C.  Koltenbah 
%  January  6,  2008 

%  This  routine  will  read  and  parse  a  SCD  file  and  create  a  reduce  data  set 
%  MAT  file  that  can  be  loaded  later  and  manipulated  from  within  Matlab. 

%  Inputs 

%  DirName  Directory  name  -  include  trailing  1 \ ’ 

%  FileName  File  name  with  '.scd'  extension 

%  Outputs 

%  none 

%  Note:  The  created  MAT  file  will  have  the  same  root  name  as  the  SCD  file, 
function  PostProcExtraction003 (DirName, FileName) 

%  Read  data  contents  of  SCD  file. 


SCDFile= [DirName  FileName]; 

If ix=l ; 

DAT A=MRS_Re a dSCDFile_ 01  (SCDFile,  lfix)  ;  %  Read  and  parse  SCD  file 

%  NOTE:  At  this  point,  DATA  contains  all  the  original  data  and  information 
%  from  the  SCD  file.  This  could  be  further  manipulated  from  within  Matlab. 

DAT  A=MRS_Ex  tract  Da  ta_  02  ( DATA)  ;  %  Create  extra  fields  of  data 

%  NOTE:  Now  DATA  includes  some  additional  enhanced  data,  namely  calibrated 
%  ZX2  and  ZY2  component  arrays.  See  MRS_ExtractData_02  for  more  details. 

%  NOTE:  One  could  now  set  a  detection  angle  delta  and  define  the  following 
%  to  be  plotted: 

%  Z  =  DATA. ZX2 . *cos (delta)  +  DATA. ZY2 . *sin (delta) 


%  Prepare  reduced  output  to  be  saved  in  MAT  file. 
k=strf ind (FileName, ' . scd’ ) ; 

RootName=FileName ( 1 : k-1 ) ;  %  Root  file  name 


X (1) =DATA. ScanStp; 

X (2) =DATA. IndexStp; 

X (3) =DATA . Excitation_Signal_Amplitude (1) ; 
X (4) =DATA . Excitation_Signal_Frequency (1) ; 


Scan  step  size  (in) 
Index  step  size  (in) 
Amplitude  (V) 
Frequency  (Hz) 


ZX=DATA . ZX2 ' ; 


Calibrated  ZX-component 


ZY=DATA . ZY2 '  ; 


Calibrated  ZY-component 


DateTime=DATA . DateTime ; 


Date  and  Time  string 


%  Save  the  MAT  file  -  (note  the  -v6  parameter) 

FileName3= [DirName  RootName  ' .mat' ] ; 

save (FileName 3 , ' RootName ’ , 1 DateTime ' , ' X ’ , 1 ZX ' , ’ ZY ' , ’ -v6 ’ ) ; 


end 


Figure  A-l  Matlab  Code  of  Function  PostProcExtraction003 


A.3  MRS  READSCDFILE  01 

This  function,  as  its  name  suggests,  will  read  an  SCD  file  and  parse  the  data  into  the 
Matlab  programming  environment.  The  code  is  listed  in  the  following  sequence  of  figures, 
Figure  A-2  through  A-5. 
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MRS  ReadSCDFile  01  -  Read  SCD  File  -  Version  01 


Author : 
Date : 


Dr.  Benjamin  E.C.  Koltenbah 
July  11,  2006 


This  routine  will  read  and  parse  a  SCD  file  generated  by  the  MRS  Testbed 
Software.  The  data  is  parsed  and  placed  in  a  DATA  structure. 


Inputs 

SCDFile 
If  ix 


SCD  File  -  Full  Path  Name 
Optional  flag  to  fix  missing  "2" 
Default:  lfix=0  (don't  fix) 


from  "Dy2"  variable 


Outputs 


DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 


DATA 

DATA 


DATA 

DATA 


Data  structure  -  field  names  match  internal  SCD  file  names 
DataVer  Data  File  Version 

SoftName  DAQ  Software  Version 

DateTime  DAQ  Date  and  Time 

Caption  Scan  caption  (if  included) 

AcqVx  Acquisition  speed  in  scan  direction  (X)  (in/s) 

AcqVy  Acquisition  speed  in  index  direction  (Y)  (in/s) 

AcqVx  Acquisition  speed  in  scan  direction  (X)  (in/s) 

AcqVy  Acquisition  speed  in  index  direction  (Y)  (in/s) 

Nx  Number  of  data  points  in  scan  direction  (X) 

Ny  Number  of  data  points  in  index  direction  (Y) 

ScanStp  Scan  Step  Size  (in) 

IndexStp  Index  Step  Size  (in) 

Dx  Requested  Scan  Area  Width  (X)  (in) 

Dx2  Actual  Scan  Area  Width  (X)  (in) 

Dy  Requested  Scan  Area  Height  (Y)  (in) 

Dy2  Actual  Scan  Area  Height  (Y)  (in) 

NumSens  Number  of  Sensors  in  Array 

NumExSig  Number  of  Excitations 

Sensor_Phase_Shif t 

Phase  Shift  (deg) 

[NumSens  NumExSig] 

Sensor_Phase_Shif t 

Quadrature  Phase  Shift  (deg) 

[NumSens  NumExSig] 

CScan_Detection_Phase 

Detection  Phase  Used  in  C-Scans  (deg) 

[NumSens  NumExSig] 

Boolean  Array  of  Active  Sensors  [NumSens] 

Sensor  Offsets  in  Array  Layout  (in  X  and  Y) 

[2  NumSens] ,  1 :X,  2 : Y 
Excitation_Signal_Active 

Booelan  Array  Active  Excitation  Signals 
[NumExSig] 

Strap  Current  (Saved  in  V  units)  [2] 

Z-Range  (VA2)  for  C-Scan  Visualization 
[2  NumExSig],  l:Z-Min,  2:Z-Max 
Sensor  Number  for  Each  Scan  Data  Value 
Numbered  from  0  to  NumSens-1,  -1  is  Empty  Data 
[Nx  Ny  NumExSig] 

Raw_Demodulated_Value 

Raw  Demodulated  Values  (VA2)  -  Uncalibrated 
[Nx  Ny  NumExSig] 

Raw_Demodulated_Quad_Value 

Raw  Demodulated  Quad  Values  (VA2)  -  Uncalibrated 
[Nx  Ny  NumExSig] 

Calibration  1  Values  [2  NumSens  NumExSig] 

1 : In-Phase  Value,  2: Quadrature  Value 
Calibration  2  Values  [NumSens  NumExSig] 


. Active_Sensors 
.Sensor  Offsets 


Strap_Current 

Z_Range 


DATA. Sensor  Index 


DATA 

DATA . Calibration_l 
DATA . Calibration  2 


Figure  A-2  Matlab  Code  of  Function  MRSReadSCDFileOl:  Header 

Figure  A-2  shows  the  header  of  the  function.  This  information  is  important  if  one  wishes 
to  understand  the  data  structures  into  which  the  scan  data  is  parsed.  The  inputs  are  a  string  of 
the  full  path  to  the  SCD  fde  and  a  logical  that  should  be  applied  to  older  data  to  fix  an  early 
minor  error  in  the  data.  The  output  is  the  DATA  structure,  the  full  contents  of  which  are 
detailed  in  the  header  in  Figure  A-2. 
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Figure  A-3  shows  the  opening  steps  of  the  file,  where  the  various  expected  variables  from 
the  SCD  file  are  listed  along  with  their  respective  data  types.  Several  variables  are  initialized 
for  the  file  read  that  is  to  follow. 


function  DATA=MRS_ReadSCDFile_01 (SCDFile, lfix) 

%  Define  lfix=0  if  not  input 
if  (nargin<2) 

If ix=0 ; 

end 


% -  Define  Variable  Names  (SCD  File  Syntax) 

VAR= { ' DataVer ' ,  ' SoftName ' ,  ' DateTime ' ,  '  Caption ' ,  .  .  . 

' AcqVx ' , ' AcqVy ' , ' TrVx ' , ' TrVy '  ,  .  .  . 

' Nx ' , ' Ny ' ,  ' ScanStp ' ,  ' IndexStp ' ,  ' StanzaStp ’ , . . . 

'Dx', ' Dx2 ' , 1 Dyf , ' Dy2 ' , . . . 

' NumSens ' ,  ' NumExSig ' ,  . . . 

'Sensor  Phase  Shif t Sensor  Quad  Phase  Shift',... 

'CScan  Detection  Phase',... 

'Active  Sensors ',' Sensor  Of f sets ',' Excitation  Signal  Active',... 
'Excitation  Signal  Amplitude ',' Excitation  Signal  Frequency',... 

' Strap  Current ' , ' Z-Range ' , ' Sensor  Index ' , . . . 

'Raw  Demodulated  Value', 'Raw  Demodulated  Quad  Value',... 
'Calibration  1 ',' Calibration  2'}; 


% -  Define  Corresponding  Data  Types 

%  F  Single  float  value  on  single  data  line 

%  S  String  on  single  data  line 

%  1  ID  Array  of  floats  (float64) 

%  2  2D  Array  of  floats  (float64) 

%  3  3D  Array  of  floats  (float64) 

%  -1  ID  Array  of  booleans  (int8) 

%  -3  3D  Array  of  integers  (int32) 

TYP= { 'F', 'S', 'S', 'S',... 

'F' , 'F' , 'F' , 'F' , . . . 

'F' , 'F' , 'F' , 'F' , 'F' , . . . 

'p  »p»  »p»  »p» 

'  F '  ^  '  F  '  ^ 

'2',  '2',  .  .  . 

'1',  .  .  . 

'-1', '2', '-1', . . . 

'1',  '1',  .  .  . 

'1',  '2',  '-3',  .  .  . 

f  O  f  !  O  ! 

'  3  '  '  '  2  '  } ;  ” 

num=length (VAR) ;  %  Number  of  variable  names 

fid=f open (SCDFile, ' r ') ;  %  Open  file  -  read  only 

DAT= [ ] ;  %  Initialize  DAT 


Figure  A-3  Matlab  Code  of  Function  MRS_ReadSCDFile_01:  List  of  Variables, 
Names,  and  Corresponding  Data  Type,  Initialization  for  File  Read 

Figure  A-4  shows  the  section  where  the  file  is  read  and  parsed.  Note  the  use  of  the 
“ieee-be”  option  in  the  ffead  functions  due  to  the  fact  that  the  machine  format  of  the  SCD  data 
is  IEEE  floating  point  with  big-endian  byte  ordering.  This  is  the  native  format  from  the 
Lab  VIEW  data  acquisition  program. 
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% -  Read  File  Until  End 

while  (~feof (fid) ) 

str=fgetl (fid) ; 
st.T=  TrimS tr{  str )  ; 
lstr=length (str) ; 

If lag=0 ; 

%  Loop  through  variable  names 
for  n=l:num 

var=VAR{n} ; 
typ=TYP { n } ; 
lvar=length (var) ; 

%  var  is  longer  than  str  ■ 
if  (lvar>lstr) ;  continue; 

If lag= (str (1 : lvar) ==var) ; 

%  Catch  redundancy  of  Dx2 
if  (lflag  &  length (DAT) ==] 
if  (lfix  &  n==16) 

str= [str (1:2)  '2' 

lstr=lstr+l ; 

end 

continued- 

end 


%  Get  next  line 

%  Trim  beginning  and  exclude  comments 
%  Length  of  line 

%  Initialize  lflag 

-  search  for  match 

%  nth  VAR 
%  nth  TYP 

%  Length  of  variable  name 

skip 

end; 

%  Check  match 

and  Dy2 
d 

str  (3 : lstr) ] ; 


%  If  matched,  read  data 
if  (lflag) 

switch  (typ) 

case  {  '  F  '  , ' S ' }  %  Single  line  data 

DAT { n } =GetVal ( str ,  typ) ;  %  Read  value  from  str 
case  {'1'}  %  ID  float  (float64) 

nn=f read ( f id, 1 , ' int32 ' , ' ieee-be ' ) ; 

DAT { n } =f read ( f id, nn, ' float 6 4 ' , ' ieee-be ' ) ; 
case  {'2'}  %  2D  float  (float64) 

nnl=f read ( f id, 1 , ' int32 ' , ' ieee-be ' ) ; 
nn2=f read ( f id, 1 , ' int32 ' , ' ieee-be ' ) ; 
xtmp=f read (fid, nnl*nn2 , ' float 64 ' , ' ieee-be ' ) ; 
DAT{n}=reshape (xtmp, [nn2  nnl]); 
case  {’-1'}  %  ID  boolean  (int8) 

nn=f read ( f id, 1 , ' int32 ' , 1 ieee-be ' ) ; 

DAT { n } =f read ( f id, nn, ' int8 ' , ' ieee-be ' ) ; 
case  {’-3'}  %  3D  integer  (int32) 

nnl=f read ( f id, 1 , ' int32 ’ , ’ ieee-be ' ) ; 
nn2=f read ( f id, 1 , ’ int32 ' , ’ ieee-be ' ) ; 
nn3=f read ( f id, 1 , ’ int32 ’ , ’ ieee-be ' ) ; 
xtmp=fread (fid, nnl*nn2*nn3, ' int32 ' , ' ieee-be ' ) ; 

DAT { n } =reshape (xtmp, [nn3  nn2  nnl]); 
case  {'3'}  %  3D  float  (float64) 

nnl =f read ( f id, 1 , ' int32 ' , ' ieee-be ' ) ; 
nn2=f read ( f id, 1 , ' int32 ’ , ' ieee-be ' ) ; 
nn3=f read ( f id, 1 , ' int32 ’ , ' ieee-be ' ) ; 
xtmp=fread (fid, nnl*nn2*nn3, ’ float 6 4 ' , ’ ieee-be ' ) ; 

DAT {n}=reshape (xtmp, [nn3  nn2  nnl]); 
otherwise  %  Unrecognized  -  set  to  null 

DAT { n } = [ ] ; 

end 

break  %  Break  for-loop 

end 

end 

end 


Figure  A-4  Matlab  Code  of  Function  MRSReadSCDFileOl:  File  Read  Loop  and 

Data  Parsing  Section 


Figure  A-  5  shows  the  last  section  of  code  for  this  function,  where  the  DATA  output 
structure  variable  is  assembled  with  all  the  SCD  variables  that  have  been  read.  There  are  two 
support  functions,  to  trim  strings  of  unnecessary  white  space  and  to  extract  numerical  values 
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from  string  sections.  The  code  has  been  highlighted  in  Figure  A-4  to  show  where  these 
support  functions  are  called  from. 


f close (fid) ; 

DATA= [ ] ; 
for  n=l : num 

str=VAR{ n } ; 

str=regexprep (str, '  ' , ) ; 

str=regexprep (str, 

DATA=set field (DATA, str, DAT {n} ) ; 

end 

end 


function  str=TrimStr ( strO ) 

str=strO ; 

l=length (str) ; 

if  (1<2);  return;  end 

if  (str  (1:2)=='  !  ') 

str=str (2:1) ; 
str=strtrim (str ) ; 

end 

idx=findstr (str, ' ! ’ ) ; 
num=length (idx) ; 
if  (num==0) ;  return;  end 
n2=idx (num) -1; 
str=str (1 :n2) ; 
str=strtrim (str ) ; 

end 


function  val=GetVal (str , typ) ; 
val= [ ] ; 

idx=f indstr (str , ' =  '); 
num=length (idx) ; 
if  (num==0) ;  return;  end 

nl=idx (num) +2 ; 
n2=length (str) ; 
str2=str (nl : n2 ) ; 

switch  (typ) 

case  { ’ F 1 } 

val=str2num (str2 ) ; 
case  {'S’} 

val=str2 ; 

end 

end 


Figure  A-  5  Matlab  Code  of  Function  MRS_ReadSCDFile_01:  Final  Assembly  of 
DATA  Output  Structure  and  Auxiliary  Functions 

A.4  MRS_EXTRACTDATA_02 

This  function  is  used  to  take  the  raw-formed  data  that  was  stored  in  the  SCD  fde,  and  is 
now  parsed  in  the  Matlab  environment.  Raw  data  is  saved  by  the  data  acquisition  software 
along  with  calibration  numbers,  as  described  in  Section  4. 1 . 1 .5.  It  is  useful  to  refer  again  to 
Equation  (32),  which  is  reproduced  here: 
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f;(0)=*;cos(0)+j>n(0) 

K=c.(X.-a.)  ■  (A-l) 

r.=c.(y,-b,) 


Recall  that  the  primed  quantities  are  calibrated  values,  and  the  unprimed  are  uncalibrated 
or  raw  data.  Given  recorded  raw  demodulated  component  values  of  {Xn ,  Yn )  and  a  set  of 

calibration  values  ( an,bn,cn )  for  each  sensor  n,  it  is  desired  to  produce  c-scans  of  quantity 
F'n{d),  which  represents  the  calibrated  demodulated  value  at  some  detection  angle.  As 
Equation  (A-l)  shows,  this  final  value  is  a  mixing  of  calibrated  components  (X',7„').  This 

function  will  add  the  calibrated  components  to  the  data  structure  DATA  as  it  is  typically 
computationally  intensive  to  perform  the  above  calibration  operation  of  Equation  (A-l) 
continuously  in  an  interactive  environment. 

The  code  for  this  function  is  listed  in  the  next  two  figures.  Figure  A-6  shows  the  header 
and  initial  steps  of  data  extraction  and  manipulation.  The  three  calibration  variables  of 
Equation  (A-l)  are  populated  on  a  two-dimensional  array  that  is  the  size  of  the  raw 
component  data  arrays.  These  must  be  correlated  properly  with  their  respective  sensor  as 
mapped  on  the  data  arrays.  There  is  an  additional  step  that  is  not  listed  here  as  part  of 
Equation  (A-l).  Due  to  the  nature  of  the  data  acquisition  hardware,  namely  the  multichannel 
A/D  converters  described  in  Section  4. 1.1.1,  the  digital  data  is  sampled  in  a  multiplexed 
fashion.  This  actually  introduces  a  small  but  constant  time  (or  phase)  lag  for  each  sensor  for 
which  a  correction  must  be  applied  in  this  step.  Therefore,  the  true  definition  of  each 
calibrated  in-phase  and  quadrature  component  requires  not  only  the  zero  subtractions  and  gain 
factors  but  also  a  small  mixing  of  the  uncalibrated  raw  components  by  this  phase  correction. 
Figure  A-6  shows  the  initialization  of  the  calibration  and  phasing  corrections  and  the 
application  of  the  calibration  part  only. 
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%  MRS_ExtractData_01  -  Extract  MRS  Data  from  DATA  Structure 

%  Author:  Dr.  Benjamin  E.C.  Koltenbah 

%  Date:  November  15,  2006 


%  This  routine  will  extract  raw  uncalibrated  data  from  the  DATA  structure 
%  and  compute  2D  arrays  of  coordinates  and  calibrated  data  components.  The 
%  results  are  stored  in  the  DATA  structure. 


function  DATA=MRS_ExtractData_02 (DATA) 


Nx=DATA . Nx; 
Ny=DATA . Ny; 


Number  of  points  in  scan  direction 
Number  of  points  in  index  direction 


dx=DATA. ScanStp; 
dy=DATA . IndexStp; 


Scan  direction  step  size  (in) 
Index  direction  step  size  (in) 


x= ( 0 : Nx-1 ) . *dx; 
y= ( 0 : Ny-1 ) .  *dy; 

[ Y, X] =meshgrid (y, x) ; 


Scan  coordinates  (in) 

Index  coordinates  (in) 

2D  arrays  of  index  and  scan  coordinates  (in) 


clx=DATA. Calibration^! (1, : , 1) ;  %  Call 

cly=DATA. Calibration^! (2 , : , 1 ) ;  %  Call 

c2=DATA. Calibration_2 ( : , 1 ) ;  %  Cal2 

phx=DATA. Sensor_Phase_Shift ( : , 1 ) ; 
phy=DATA. Sensor_Quad_Phase_Shif t ( : , 1)  ; 
phx=phx . *pi . / 1 8  0 . 0 ; 
phy=phy . *pi . /180 . 0; 


values  for  in-phase  demod  values  (VA2) 
values  for  quadrature  demod  values  (VA2) 
values 

%  Phase  shift  for  in-phase  (deg) 

%  Phase  shift  for  quadrature  (deg) 

%  Convert  from  (deg)  to  (rad) 


ZX=DATA . Raw_Demodulated_Value ( : , : , 1) ;  %  Raw  in-phase  demod  values  (VA2) 

ZY=DATA. Raw_Demodulated_Quad_Value ( : , : , 1 ) ;  %  Raw  quadrature  demod  values  (VA2) 


IDX=DATA. Sensor_Index ( : , : , 1 ) +1 ;  %  Sensor  index  2D  map 


sz=size ( ZX) ; 

ZX=ZX  (  : )  ; 

ZY=ZY ( : ) ; 

IDX=IDX ( : ) ; 

ClX=zeros (size (IDX) ) ; 
C1Y=C1X; 

C2=C1X; 

PHX=C1X; 

PHY=C1X; 

idx=f ind ( IDX>0 )  ; 

C1X (idx) =clx (IDX (idx)  )  ; 
C1Y (idx) =cly (IDX (idx) )  ; 
C2 (idx) =c2 ( IDX (idx) ) ; 
PHX (idx) =phx (IDX (idx) ) ; 
PHY (idx) =phy (IDX (idx) ) ; 

ZX=C2 . * (ZX-C1X) ; 

ZY=C2 . * (ZY-C1Y) ; 


2D  array  size 

Reshape  to  single  column  arrays 

Initialize  full  cal  and  phase  arrays 


Define  non-zero  indexed  points 
Define  full  cal  and  phase  arrays 


Calibrated  in-phase  demod  values  (VA2) 
Calibrated  quadrature  demod  values  (VA2) 


Figure  A-6  Matlab  Code  of  Function  MRS_ExtactData_01:  Header,  Data  Parsing, 
and  Initialization  of  Calibration  and  Phase  Corrections 


Figure  A-7  shows  the  Matlab  code  for  the  application  of  the  phasing  correction  and  final 
assembly  of  the  updated  DATA  structure. 
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%  Reshape  2D  arrays 
X=reshape  (X, sz ) ; 

Y=reshape (Y, sz) ; 

ZX=reshape (ZX, sz) ; 

ZY=re shape ( ZY, sz ) ; 

PHX=reshape (PHX, sz ) ; 

PHY=reshape (PHY, sz) ; 

%  Here  modified  ZX2  and  ZY2  are  defined  to  eliminate  need  to  use 
%  PHX  and  PHY  in  evaluations 
ZX2=ZX. *cos (PHX)  -  ZY. *sin (PHY) ; 

ZY2=ZX. *sin (PHX)  +  ZY . *cos ( PHY) ; 


Store  2D  arrays  in  DATA  structure 


DATA . X=X ; 

DATA . Y=Y ; 

DATA . ZX=ZX; 
DATA . ZY=ZY; 
DATA. PHX= PHX; 
DATA. PHY=PHY; 


x-  or  scan-coordinates  (in) 
y-  or  index-coordinates  (in) 
Calibrated  in-phase  component  (VA2) 
Calibrated  quadrature  component  (VA2) 
In-phase  comp  phase  adjust  (rad) 

Quad  comp  phase  adjust  (rad) 


DATA . ZX2=ZX2 ; 
DATA. ZY2=ZY2; 


Phase  adjusted  ZX 
Phase  adjusted  ZY 


end 


Figure  A-7  Matlab  Code  of  Function  MRS_ExtactData_01:  Final  Assembly  of 

Updated  DATA  Structure 


A.5  MATLABEXAMPLE001 


This  function  is  provided  as  an  example  of  how  the  saved  results  from  the  previous  three 
functions  can  easily  be  read  and  manipulated  in  the  Matlab  environment.  Figure  A-  8  shows 
the  code  for  this  example.  Here,  the  SCD  file  “S- 2  F4  00  ALL  06  .mat”  is  loaded  into 
Matlab,  and  in  just  a  few  lines  of  code  a  c-scan  figure  of  calibrated  values  is  plotted.  The 
MAT  file  was  created  in  the  step  listed  in  Section  A.2  using  PostProcExtraction003, 
which  internally  used  both  MRS_ReadSCDFile_01  and  MRS_ExtractData_02. 

The  MAT  file  is  loaded,  and  the  subset  of  data  contained  therein  is  used  to  construct  a 
magnitude  plot,  where 


X?+Y'2 


(A-2) 
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%  MatlabExampleOOl  -  Example  Using  Matlab  Routines  for  MRS  Data  Processing 


%  This  routine  shows  the  simplest  means  to  load  MRS  scan  data  from  reduced 
%  data  set  MAT  files.  The  original  SCD  files  can  be  loaded  via  other 
%  routines.  Refer  to  PostProcExtraction003 . m  for  more  information  on  this. 


function  MatlabExampleOOl 
%  Load  data 

MATFile= '  S-2  F400  ALL  06. mat'; 
load (MATFile) ; 


Data  now  loaded: 

RootName  File  name  (without  extension) 

DateTime  Acquisition  date  and  time  string 

X  Array  of  scalar  parameters  (see  below) 

ZX  X-  or  in  phase  component  array 

ZY  Y-  or  quadrature  component  array 


dx=X ( 1 ) ; 
dy=X (2 ) ; 
Amp=X (3) ; 
Freq=X ( 4 ) 


%  Scan  step  size  (in) 

%  Index  step  size  (in) 

%  Driver  amplitude  (V) 

%  Dirver  frequency  (Hz) 


%  In  this  example,  the  magnitude  is  plotted 


ZZ=sqrt (ZX. A2  +  ZY.A2);  %  Magnitude 


%  Define  XX  and  YY  coordinate  2D  arrays  for  surf  plot 
Nx=size  (ZZ, 2) ; 

Ny=size ( ZZ , 1 ) ; 
x= ( 0 : Nx-1 ) . *dx; 
y= (0 :Ny-l) . *dy; 

[XX, YY] =meshgrid (x, y) ; 


%  Make  the  surface  plot 

figure (1) 

h=surf (XX, YY, ZZ) ; 

set (h, ’ EdgeColor ' , ' none ' ) ; 

view (0,90); 

set (gca, ' DataAspectRatio ' , [ 1  1  1]); 

axis  tight; 

xlabel ( ' X  (in) ' ) ; 

ylabel ( ' Y  (in) ' ) ; 

title ( { RootName , DateTime } ) ; 


end 


Figure  A-  8  Matlab  Code  of  Function  MatlabExampleOOl 


The  Matlab  plot  results  are  shown  in  Figure  A-  9. 
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S-2  F400  ALL  06 
2007-04-26  15:02:26 


X  (in) 


Figure  A-  9  Matlab  Plot  Results  from  MatlabExampleOOl 


A.6  CONCLUDING  REMARKS 

The  Matlab  plot  shown  in  Figure  A-  9  is  perhaps  itself  not  useful;  however,  the  purpose 
of  describing  these  routines  and  showing  a  simple  example  is  to  demonstrate  the  utility  of 
post-processing  MR  sensor  data  in  a  graphical  interactive  environment.  As  was  mentioned 
earlier  at  the  end  of  Section  4.1.2,  the  interpretation  of  this  kind  of  data  is  complex  and 
requires  these  kinds  of  tools  to  at  times  programmatically  display  the  data  in  different  ways. 

Appendix  B  describes  the  development  of  such  an  interactive  graphical  user  interface 
environment  in  Mathematica  for  the  display  and  manipulation  of  MR  sensor  data. 
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APPENDIX  B 


Appendix  B  contains  a  write-up  that  Dr.  Ben  Koltenbah  (Boeing)  provided  with  the  data 
he  acquired  on  the  S-l,  S-2,  S-3  and  S-4  test  standards.  He  developed  a  post-processing  tool 
that  could  be  used  to  view  and  post-process  the  data  results.  This  write-up  has  been  slightly 
modified  from  his  original  to  better  fit  the  format  of  this  report. 


MR  Sensor  Data  Post-Processing  Tool 

PostProcess04B 
Dr.  Benjamin  E.  C.  Koltenbah 

Applied  Physics  -  Phantom  Works  -  The  Boeing  Company 
Seattle,  WA  USA 
January  21,  2008 


B.l  INTRODUCTION 

This  is  a  brief  manual  that  gives  instructions  for  use  of  PostProcess04B,  a 
Mathematica  notebook  that  I  developed  to  aid  in  viewing  MR  Sensor  data.  This  tool  is  still 
very  much  under  construction,  which  means  that  there  are  quite  a  number  of  bugs  and 
inconveniences  in  the  program.  However,  as  a  whole  the  program  has  proven  to  be  a  nice  tool 
for  fast  post-processing  analysis  and  should  serve  to  be  a  food  framework  for  further 
developments. 

By  no  means  do  I  propose  this  Mathematica  format  as  an  analysis  tool  for  incorporation 
into  data  acquisition  packages.  However,  I  do  think  that  it  serves  as  a  standalone,  rapid 
development  framework  for  developing  new  post-processing  algorithms. 

The  central  portion  of  the  notebook  is  the  Manipulate  function,  which  is  a  new  feature 
within  Mathematica  6.  This  constructs  a  graphical  user  interface  panel  in  which  “Control 
Objects,”  such  as  sliders  and  input  fields,  can  be  used  to  control  variable  values  and  thus 
manipulate  the  objects  within  the  panel.  I  built  this  tool  partially  as  a  means  to  learn  more 
about  the  Manipulate  function  and  test  its  utility  for  my  work.  I  am  mostly  satisfied  with  the 
outcome  thus  far  and  will  pursue  further,  more  advanced  applications  in  the  future.  As  a 
consequence  of  this  developmental  environment,  one  will  need  a  copy  of  Mathematica  6  in 
order  to  run  and  change  my  notebook.  I  will  explore  the  possibility  of  using  the  free 
Mathematica  reader  to  use  the  tool  in  the  future. 
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Finally,  note  that  this  notebook  is  being  provided  by  me  “as  is,”  which  means  that  I  make 
no  guarantee  as  to  its  accuracy  or  ability  to  run  on  all  platforms.  This  is  not  a  required 
deliverable  on  a  program;  therefore,  it  is  not  “certified  code.”  I  do  provide  it  to 
collaborators;  however,  in  that  I  wish  for  them  to  have  access  to  any  tools  I  have  developed, 
no  matter  how  primitive  or  unfinished  they  may  be.  My  time  and  ability  to  assist  in  using  the 
program  is  somewhat  limited;  however,  I  hope  that  this  manual  and  a  few  encouraging  words 
will  be  sufficient  to  get  one  started. 

B.2  RUNNING  THE  PROGRAM 

To  run  the  notebook,  one  must  load  it  into  Mathematica  6  and  use  the  “Evaluate 
Notebook”  command  under  “Evaluation”  on  the  top  menu  toolbar.  This  applies  the  various 
definitions  for  the  object  tools  and  creates  the  graphical  user  interface  pane.  All  interaction 
from  this  point  on  should  only  need  to  be  done  on  the  panel,  which  appears  at  the  very  bottom 
of  the  notebook.  If  one  wanted  to  change  some  of  the  tool  definitions  or  add  to  the  notebook, 
he  would  probably  want  to  save  the  file  and  then  evaluate  the  notebook  again.  I  have  taken 
steps  to  insure  that  re-evaluations  can  be  done  cleanly,  however  note  that  any  information  on 
the  panel  will  be  lost  and  the  session  will  start  with  a  blank  picture  again. 

The  remaining  sections  of  this  manual  cover  the  various  control  object  blocks  that  appear 
along  the  left  side  of  the  panel. 

B.2.1  Control  Selection 

I  provided  these  check  boxes  to  display  and  hide  the  various  control  groups  as  not  all  of 
them  need  be  on  the  screen  at  once. 

B.2.2  Load  Data  Button 

This  is  likely  the  starting  point  of  any  session.  A  fde  dialog  box  appears  that  asks  for 
selection  of  a  *.MAT  file.  This  is  a  Matlab  format  file  that  contains  sufficient  data  for 
Mathematica  to  construct  c-scans. 

The  original  scan  data  file  format  is  *.SCD,  which  is  created  by  Lab  VIEW.  Elsewhere 
(see  Appendix  A)  I  provide  Matlab  routines  that  read  and  parse  these  files  completely  and 
thus  load  the  data  into  the  Matlab  environment.  Unfortunately,  Mathematica  is  unable  to  read 
MAT  files  that  are  saved  with  the  whole  entire  data  set.  I  had  to  develop  other  Matlab 
routines  that  create  a  subset  of  the  data.  I  believe  that  these  subsets  contain  sufficient  data  for 
viewing  the  more  pertinent  information  in  this  Mathematica  tool.  The  full  original  data  set  is, 
of  course,  still  available,  and  I  may  be  able  to  include  more  and  more  of  the  original 
information  as  needed.  Suffice  it  to  say  that  calibrated  X  and  Y  component  data  is  available 
along  with  minimal  information  about  the  original  scan. 

Once  a  file  is  selected,  the  tool  will  adjust  itself  to  show  a  plot  of  the  X-component  (in- 
phase)  of  the  data.  One  annoying  quirk  is  that  the  extents  of  the  plot  may  not  be  scaled  to  fit 
the  plot.  In  order  to  adjust  this,  one  need  only  select  the  Y-Range  (in)  tab  and  possibly  the  X- 
Range  (in)  tab.  The  scaling  ought  to  be  done  automatically  at  this  point.  I  have  not  found  a 
way  to  fix  this  yet,  but  it  should  only  be  a  minor  annoyance. 
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The  caption  and  data  acquisition  date  and  time  should  appear  in  the  title  of  the  plot. 
Other  information  that  may  appear  is  the  data  functions  being  shown,  the  detection  angle 
value,  and  smoothing  settings.  I’ll  write  more  on  these  later. 

The  tool  should  now  be  ready  for  manipulation  of  this  data. 

B.2.3  Save  Image  Button 

At  any  time  an  image  of  the  plot,  as  shown,  can  be  saved  to  file.  This  will  save  in  the 
JPEG  format.  Note  that  some  adjustments  are  made  to  the  font  size  in  the  image.  One  may 
find  it  necessary  to  adjust  this  functionality  based  on  the  computer  system  being  used.  I 
suggest  that  a  few  trial  images  first  be  saved  followed  by  adjustments  as  needed. 

B.2.4  Function  Selection 

Here  the  function  to  plot  is  selected.  Options  are  as  follows: 

Mix  Comp  -  This  plots  the  mixed  components  as  Z  =  X  cos  8  +  Y  sin  8 ,  the  familiar 
mixing  formula  for  MRS  data. 

Magnitude  -  The  magnitude  of  the  components  Z  =  Va2  +  72  . 

Phase  -  The  phase  (in  degrees)  of  the  components  relative  to  the  excitation  signal 
6  =  arctan (7/x),  where  6  spans  -180°  to  +180  °. 

Abs  Mix  Comp  -  This  is  the  same  as  1  but  shows  the  absolute  value  instead  |z| . 

B.2.5  Gradient  Selection 

Here  the  gradient  of  the  function  can  be  chosen.  Options  are  as  follows: 

None  -  No  gradient. 

Grad-X  -  The  derivative  (difference)  along  the  x-direction. 

Grad-Y  -  The  derivative  (difference)  along  the  y-direction. 

I  suggest  that  the  data  be  sufficiently  smoothed  for  less  noisy  gradient  plots.  See  the 
smoothing  control  group  below. 

B.2.6  Detection  Angle 

For  “Mix  Comp”  and  “Abs  Mix  Comp”  functions,  the  detection  phase  angle  can  be 
adjusted.  If  the  “Data  Z-Range”  group’s  “auto”  checkbox  is  selected,  the  plot’s  contrast  will 
adjust  itself  to  match  the  minimum/maximum  limits  of  the  data  being  displayed. 

Note  that  most  sliders  include  “min”  and  “max”  input  fields,  where  the  extents  of  the 
slider  can  be  changed.  These  can  be  set,  for  example,  to  span  a  much  tighter  range  of 
detection  angles  for  better  fine-tuning  of  the  plot.  There  is  also  a  “stp”  input  field,  where  the 
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step  size  can  be  adjusted.  When  the  left  and  right  arrows  are  pressed,  the  value  will  change 
up  or  down  by  this  step  size.  Finally,  there  is  a  “val”  input  field  (yellow)  where  the  value  of 
the  slider  can  also  be  set. 

B.2.7  Color  Scheme 

This  tool  can  be  used  to  change  the  color  palette  of  the  plot.  The  colors  listed  are  those 
provided  by  Mathematica.  Presumably,  one  could  add  more  to  this  list;  however,  I  haven’t 
attempted  to  do  so  as  of  yet. 

B.2.8  Data  Z-Range 

Here  the  data  contrast  can  be  adjusted  with  two  values  on  this  slider.  Data  outside  of  this 
range  will  be  presented  as  white  (black)  when  above  (below)  the  range. 

In  addition  to  the  aforementioned  slider  controls,  there  is  also  a  “lock”  check  box  for 
locking  the  difference  between  the  minimum  and  maximum  values.  Now  when  adjusted,  this 
difference  will  remain  the  same  as  the  extents  move  in  tandem. 

There  is  also  an  “auto”  check-box  which  causes  the  extents  of  the  graph  to  self-adjust  as 
other  controls  change  the  plot.  I  suggest  that  this  be  checked  until  one  is  finished  adjusting, 
say,  the  detection  angle  or  extents  of  the  plot  and  is  then  ready  to  fine-tune  the  contrast  range. 

B.2.9  X-Range  (in)  and  Y-Range  (in) 

These  controls  set  the  x-  and  y-extents  of  the  plot  itself.  Note  that  these  do  not  adjust  the 
extent  of  the  data  being  displayed,  although  some  data  may  be  clipped  and  truncated.  Z-range 
contrast  is  auto-adjusted  based  on  all  the  data  selected,  not  on  these  plot  ranges. 

B.2.10  Data  X-Index  Range  and  Data  Y-Index  Range 

These  sliders  adjust  the  subsets  of  the  matrix  data  to  be  displayed.  Only  the  subset 
defined  here  will  be  shown  and  used  to  determine  the  auto-adjusted  contrast  range. 

B.2.11  X-Index  Smoothing  and  Y-Index  Smoothing 

These  controls  are  used  to  adjust  the  window  size  for  running  averages  along  both  the  x- 
and  y-directions.  From  my  experience  I  suggest  using  values  of  around  5  for  both  these 
controls.  By  setting  this  to  1,  one  can  see  the  unsmoothed  data. 

B.2.12  Lighting 

The  plotting  command  that  I  use  in  this  tool  isReliefPlot,  which  is  a  new  feature  in 
Mathematica  6.  This  function  can  be  used  to  plot  2D  data  and  provide  lighting  and  shading  to 
enhance  the  topology.  It  appears  to  have  been  originally  intended  for  visualizing  terrain 
contour  data  and  provides  realistic  color  and  shading.  Any  2D  matrix  data  can  be  plotted  in 
this  manner,  and  it  may  be  of  some  value  to  explore  the  lighting  and  shading  capabilities.  By 
default,  I  have  adjusted  these  parameters  to  essentially  turn  off  the  shading  and  just  show  fully 
lit  data. 
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Lighting  Angle  -  This  sets  the  angle  of  the  light  source. 

Lighting  Elevation  -  This  sets  the  elevation  angle  of  the  light  source. 

Lighting  Altitude  -  This  sets  the  relative  altitude  above  the  data  of  the  light  source. 

Image  Size  -  This  sets  the  image  size  in  pixels,  defaulted  to  900.  One  may  need  to 
adjust  this  to  get  a  better  fit  of  the  data  on  the  screen. 

B.2.13  Geometry 

I  developed  this  tool  primarily  for  viewing  particular  data  from  the  Boeing  S-l,  S-2,  S-3, 
and  S-4  test  standards.  This  tool  group  can  be  used  to  view  and  manipulate  the  basic 
geometry  of  these  standards  on  top  of  the  data.  More  geometries  may  be  added  in  the  future 
based  on  need. 

Geometry  Selection  -  The  choices  are  None,  S-l/S-3,  S-2/S-4.  One  can  also  show  or 
not  show  the  text  with  a  text  box. 

Color  -  This  selects  the  color  of  the  geometry  lines.  It  may  be  beneficial  to  use  a  light 
color  if  the  data  appears  dark. 

Line  Width  -  This  sets  the  line  width  in  points.  It  may  be  useful  to  increase  the  line 
width  just  before  saving  an  image. 

X-Shift  &  Y-Shift  -  There  is  no  automatic  correlation  of  data  position  with  the 

geometry.  These  tools  can  be  used  to  shift  the  geometry  accordingly.  One 
may  need  to  adjust  the  extents  of  the  plot  in  order  to  see  the  full  geometry. 

B.3  EXAMPLE 

Figure  B-l  shows  an  example  of  post-processing  of  the  S-l  standard  at  a  frequency  of 
350  Hz.  Here  I  have  adjusted  the  x-  and  y-extents  of  the  plot  to  show  the  entire  geometry. 

The  data  is  clipped  to  show  the  right  and  left  lobes  of  the  Ti  fasteners  with  cracks  from  the 
inside  and  outside  rows,  respectively.  I  have  adjusted  the  detection  angle  and  contrast  such 
that  the  four  inside  cracks  can  be  seen  well  (saturated  black)  and  at  least  two  of  the  outside 
row  cracks  can  be  seen  well  (saturated  white).  There  is  minimal  saturation  for  Tl  1,  and  no 
discemable  crack  on  T3.  Interestingly,  T7  shows  up  more  clearly  than  Til  even  though  its 
notch  crack  is  supposedly  smaller. 
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Figure  B-l  Example  of  S-l  Sample  Post-Processing 


Figure  B-2  shows  the  same  post-processing  parameters  but  with  zoomed  in  extents. 
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Figure  B-2  Zoomed  Re-plot  of  Figure  A-2 


Figure  B-3  shows  the  corresponding  saved  JPEG  image.  One  can  see  here  that  it  may  be 
useful  to  adjust  the  font  sizes  and  line  weights  to  make  higher  quality  plots.  This  can  be 
adjusted  as  needed  by  working  on  the  various  definitions  in  the  notebook. 
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S-l  F350  ALL  01  :  2007-04-26  22:09:35 
Mix  Comp  6  =  63.5°  :  Smoothing  =  (5,5) 


Figure  B-3  Saved  Image  from  Figure  A-3  Post-Processing  Session 


B.4  SUMMARY 


In  summary,  PostProccess04B  provides  a  nice  format  within  the  Mathematica 
environment  to  view  and  process  MR  Sensor  scan  data,  particular  during  post-processing 
algorithm  development.  It  is  my  hope  that  others  can  take  the  work  done  here  and  extend  it  to 
more  powerful  post-processing  tools  and  formats. 
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LIST  OF  ACRONYMS 


BEM 

boundary  element  method 

BIM 

boundary  integral  method 

ECT 

eddy  current  testing 

FDM 

finite  difference  method 

FEM 

finite  element  method 

GMR 

giant  magnetoresistive 

JSF 

Joint  Strike  Fighter 

MAG 

magnitude  (amplitude-based  features) 

MAUS 

Mobile  Automated  Scanner 

MR 

magnetoresistive 

MSU 

Michigan  State  University 

NDE 

Nondestructive  Evaluation 

ODA 

optimum  detection  angle 

PDE 

partial  differential  equations 

POD 

probability  of  detection 

SNR 

signal-to-noise  ratio 

TFQMR 

Transpose  Free  Quasi-Minimal-Residual 

VIM 

volume  integral  methods 
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